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Abstract 

In this report, instead to give comprehensive review on two important re- 
search fields during my first term postdoctoral working period: Higgs and FCNC 
physics, I will collect part of my recently works on it. Charged Higgs is the 
distinguished signature of new physics, in this report, I review my two works on 
charged Higgs associated production with top quark and W boson at hardron col- 
liders. Our researches show that these two charged Higgs production mechanisms 
are important channel not only in finding charged Higgs, but also in studying the 
quantum structure of new physics. Flavor Changing Neutral Current (FCNC) 
processes are forbidden at tree level in the Standard Model (SM), so they act as 
the ground to test quantum structure of the SM and also very important channel 
in finding new physics beyond the SM. In this report, I focus on the studies on 
FCNC processes on linear colliders and in B-factories. 
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Part I 

Higgs Physics 

At the moment, because of lack of imagination, one cannot do much 
more than try to calculate effects due to the Higgs system in order to 
make comparisons with experiments results. 

-adapted from M. Veltman 1997 "Reflections on the Higgs system". 



1 Preface 

The standard model (SM) gives an excellent theoretical description of the strong and 
electro-weak interaction. This theory which is based on an £77(3) x £77(2) x U(l) gauge 
group, has been proven extraordinarily robust. Albeit its success, the SM still has one 
part untested which is the mechanism of electroweak symmetry spontaneous breaking 
(EWSB), through which the gauge bosons and fermions gain their masses. In the SM, 
EWSB is realized through one fundamental scalar field - Higgs field. After EWSB, the 
physical world is left with one neutral Higgs boson. The mass of Higgs boson is not 
predicted by the theory and can only be determined by high energy experiments. 

Although the SM is robust, there are theoretical aspects of the SM, e.g. triviality |2| 
and naturalness |§ etc., which suggested the need for new physics. In addition, there 
are certain open questions within the SM, such as too many free parameters, origin 
of CP violation and flavor problem etc., whose answer can only be found by invoking 
physics beyond the SM. 

Among various new physics, supersymmetry (SUSY)||, especially minimal super- 
symmetrical standard model (MSSM), is the most elegant candidate. In order to 
preserve the SUSY and keep theory anomaly free, in the MSSM, there should be intro- 
duced two Higgs doublets to break the electroweak symmetry. After SUSY breaking 
and EWSB, there are five physical Higgs bosons: three neutral Higgs and two charged 
Higgs bosons. To find the Higgs bosons predicted by the SM and the MSSM and study 
theirs properties are the primary goal of present and next generation colliders for both 
theoretical and high energy experimental scientists. 

The Higgs masses are not predicted by the SM and MSSM (in the MSSM, there 
is a theoretical upper limit for lightest Higgs boson < 140 GeV), which can only be 
determined by experiments. The results coming from direct search for the Higgs in the 
process e + e~ — > ZH at LEP 200 are, for the SM Higgs boson || 

m H > 107.7GeU (95% C.L.) 

which is compatible with the result of the SM fit of all precision data || 

M H = (761^) GeV 
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or 

M H < 262 GeV (95% C.L.) ; 
for MSSM CP-even Higgs boson (tan/3 > 1) 

m H > 85.2 (95% C.L.). 

At upgraded Fermilab Tevatron (Run III), the mass of SM Higgs boson can be 
pushed up to ~ 180GeV combined subprocesses qq' — > WH and gg — > H |7|]. And, 
it is commonly thought that the combination of large hadron collider (LHC) and next 
linear collider (NLC) will cover the mass range of the Higgs boson up to 1 TeV or so. 

In the above, I have given a brief general review of this field. In the following, I 
shortly describe our works related to this topic: 

• The Higgs boson production in 77 collisions at NLC M 

High energy 77 collision is the collision mode realized at the NLC with almost 
the same center-of-mass energy and luminosity, and provide more clean place to 
study the properties of the Higgs bosons. In the framework of the MSSM, we 
studied the Higgs boson production in 77 collisions at NLC. Especially, the light 
neutral Higgs boson pair production involves many Feynman diagrams arising 
both from general two-Higgs-doublet model particles and the supersymmetrical 
virtual particles. We found the total cross section for the Higgs boson pair pro- 
duction is sensitive to the model parameters, such as tan/3, triple soft breaking 
terms A t , A b and the Higgs boson masses etc. 

• Higgs boson associated production with W at hadron colliders [[| 

Before the LHC comes to operate, to discuss the Higgs discovery potential at 
present collider Fermilab Tevatron is an urgent task. With the integrated lumi- 
nosity through the process PP — > qq 1 — > WH followed by H — > bb and 
W —>■ £u, Tevatron can find the mass of the Higgs boson up to 125GeV. In these 
works, we have studied the Yukawa corrections arising from the top loop as well 
as the leading electroweak corrections including the Higgs contributions besides 
the top contributions. And we found that in the SM, the corrections are small 
and at most few percent; however, in the MSSM, the corrections could reach tens 
of percent in the favorable parameters space. 



Higgs boson discovery potential through bg channel [[10 1 

For hadron colliders, especially LHC, the gluon distribution grows rapidly, it 
may play an important role in producing Higgs boson in particular for large 
tan p because the couplings of down-type quarks with Higgs can be enhanced in 
this case. We study the Higgs boson discovery potential through bg channel for 
both neutral and charged Higgs bosons. Indeed, we found that it is possible to 
find the SUSY neutral Higgs boson at Tevatron if tan/5 > 10. For charged Higgs 
production, we also calculated Yukawa correction, and found the magnitude of 
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the radiative corrections can exceed -20% and not sensitive to the mass of charged 
Higgs boson, these effects could be observable in the experiments. 

Higgs boson production at NLC |TTJ 

We have also studied the Higgs production at the NLC associated with W boson 
or heavy quarks. The WH production is the loop-induced process, and we found 
that the cross section can reach 1 fb, but decease rapidly with the increment of 
the Higgs mass. The tiH and bbH production have been considered by many 
groups, we re-study this process under the framework of M-theory. Our results 
show that the cross sections are sensitive to the model parameters. 



Radiative Higgs boson decay beyond the standard model |T2 



At LHC, the decay mode H — > 77 is used in searching Higgs boson for the 
intermediate mass Higgs boson. However, this searching strategy is suffered by 
the low decay rate of this mode. In this work, we study the possibility of using 
H — > //7 in searching intermediate mass Higgs boson where / represent light 
fermions. Our study shows that, at least, this channel can be used as discriminant 
between SM and MSSM for a wide range of parameter space. 

In the first section of this part, the supersymmetric electroweak correction for bg — > 
tH~ at hardron colliders will be presented in details; in second section, we will study 
the process of bb —>■ W~H + at CERN Large Hardron Collider in supersymmetrical 
model. 
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2 Supersymmetric Electroweak Corrections to Charged 
Higgs Boson Production in Association with a 
Top Quark at Hadron Colliders 

ABSTRACT 

We calculate the 0{a ew m^^ j /m\ v ) and 0(a ew m^ b ^/m^ v ) supersymmetric elec- 
troweak corrections to the cross section for the charged Higgs boson production 
in association with a top quark at the Tevatron and the LHC. These corrections 
arise from the quantum effects which are induced by potentially large Yukawa cou- 
plings from the Higgs sector and the chargino-top(bottom)-sbottom(stop) couplings, 
neutralino-top(bottom)-stop(sbottom) couplings and charged Higgs-stop-sbottom 
couplings. They can decrease or increase the cross section depending on tan/3 but 
are not very sensitive to the mass of the charged Higgs boson for high tan (3. At 
low tan /?(= 2) the corrections decrease the total cross sections significantly, which 
exceed — 12% for m#± below 300GeT^ at both the Tevatron and the LHC, but 
for m#± > 300GeT^ the corrections can become very small at the LHC. For high 
tan/3(= 10,30) these corrections can decrease or increase the total cross sections, 
and the magnitude of the corrections are at most a few percent at both the Tevatron 
and the LHC. 

2.1 Introduction 

There has been a great deal of interest in the charged Higgs bosons appearing in the two- 
Higgs-doublet models(THDM)[l], particularly the minimal supersymmetric standard 
model(MSSM) [2] , which predicts the existence of three neutral and two charged Higgs 
bosons h,H,A, and H ± . When the Higgs boson of the Standard Model(SM) has a 
mass below 130-140 Gev and the h boson of the MSSM is in the decoupling limit 
(which means that is too heavy anyway to be possibly produced), the lightest 
neutral Higgs boson may be difficult to distinguish from the neutral Higgs boson of 
the standard model(SM). But charged Higgs bosons carry a distinctive signature of 
the Higgs sector in the MSSM. Therefore, the search for charged Higgs bosons is very 
important for probing the Higgs sector of the MSSM and, therefore, will be one of the 
prime objectives of the CERN Large Hadron Collider (LHC). At the LHC the integrated 
luminosity is expected to reach L = 100/6 -1 per year in the second phase. Recently, 
several studies of charged Higgs boson production at hadron colliders have appeared in 
the literature[3,4,5]. For a relatively light charged Higgs boson, m H ± < m t — m&, the 
dominate production processes at the LHC are gg — > tt and qq — > tt followed by the 
decay sequence t — > bH + — > 6t + z/ t [6]. For a heavier charged Higgs boson the dominate 
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production process is gb — > tH~ [7,8,9]. Previous studies showed that the search for 
heavy charged Higgs bosons with m^± > m t + m b at a hadron collider is seriously 
complicated by QCD backgrounds due to processes such as gb — > tib, gb — > ttb, and 
gg — > ttbb, as well as others process[8]. However, recent analyses[10,ll] indicate that 
the decay mode H + — > r + z/ provides an excellent signature for a heavy charged Higgs 
boson in searches at the LHC. The discovery region for is far greater than had 
been thought for a large range of the (m^±,tan/3) parameter space, extending beyond 
m#± ~ ITeV and down to at least tan /3 ~ 3, and potentially to tan /? ~ 1.5, assuming 
the latest results for the SM parameters and parton distribution functions as well as 
using kinematic selection techniques and the tau polarization analysisfll]. Of course, 
it is just a theoretical analysis and no experimental simulation has been performed to 
make the statement very reliable so far. 

The one- loop radiative corrections to H~t associated production have not been 
calculated, although this production process has been studied extensively at tree- 
level[7,8,9]. In this paper we present the calculations of the 0{a ew m\^ j /m\ v ) super- 
symmetric(SUSY) electroweak corrections to this associated H~t production process at 
both the Fermilab Tevatron and the LHC in the MSSM. These corrections arise from 
the quantum effects which are induced by potentially large Yukawa couplings from 
the Higgs sector and the chargino-top(bottom)-sbottom(stop) couplings, neutralino- 
top(bottom)-stop(sbottom) couplings and charged Higgs-stop-sbottom couplings which 
will contribute at the 0(a ew mt b Jm^ v ) to the self-energy of the charged Higgs boson. 
In order to get a reliable estimate this process has to be merged with the related gluon 
splitting contribution gg — > H~tb. This leads to a suppression by about 50% at LO[12]. 
However, the complete one-loop QCD corrections are probably more important, but 
not yet available. 



2.2 Calculations and formulas 

The tree-level amplitude for gb — > tH~ is 

M = M ( S S) + M ( S* ) , (1) 

where and represent the amplitudes arising from diagrams in Fig. 1(a) and 
Fig. 1(6), respectively. Explicitly, 



M, 



and 



W9s 



V2m w (s — ml) 
-m 6 tan/3y 4 /;Pfl]w(p 6 )e At (A;)7?., 



u(p t ) [2m t cot Pp^Pl + 2m b tan (3p%P R - m t cot (3^ ]kP L 



(2) 



Mi 



(*) 



W9s 



ml 



\^2m w (t — 
-m b t a nPYfiPRMp b )e^k)T«, 



u(p t ) [2m t cot !3ptP L + 2m b tan/3rfP H - m t cot @y jkP L 



(3) 
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where T a are the SU(3) color matrices and s and t are the subprocess Mandelstam 
variables defined by 

s = (p b + kf = (pt+PH-) 2 , 

and 

t = (p t - kf = (p H - - p b ) 2 . 

Here the Cabbibo-Kobayashi-Maskawa matrix element V C KM[bt] has been taken to be 
unity. 

The SUSY electroweak corrections of order 0(a ew m 2 ^/m 2 v ) and 0(a ew m^/m^ v ) 
to the process gb — > H~t arise from the Feynman diagrams shown in Figs.l(c)-l(v) 
and Fig. 2. We carried out the calculation in the t'Hooft-Feynman gauge and used di- 
mensional reduction, which preserves supersymmetry, for regularization of the ultravi- 
olet divergences in the virtual loop corrections using the on-mass-shell renormalization 
scheme [13], in which the fine-structure constant a ew and physical masses are chosen to 
be the renormalized parameters, and finite parts of the counterterms are fixed by the 
renormalization conditions. The coupling constant g is related to the input parameters 
e,m w , and m z by g 2 = e 2 /s 2 u and = 1 — m^/m|. The parameter (5 in the MSSM 
we are considering must also be renormalized. Following the analysis of ref. [14], this 
renormalization constant was fixed by the requirement that the on-mass-shell H + lvi 
coupling remain the same form as in Eq.(2) of ref. [14] to all orders of perturbation the- 
ory. Taking into account the 0(a ew m 2 ^/m 2 v ) Yukawa corrections, the renormalized 
amplitude for the process gb — > tH~ can be written as 

M ren = M(S s) + Mjf ] + 5M Vl(s) + 5M Vl(t) + 5M s(s) + 5M s(t) + 5M v ^ s) 

+5M V2(t) + 5M b(s) + 5M m = M (s) + M (t) + ^ 5M l , (4) 

i 

where 5M V ^ S \ 6M Vl ®, 5M s( - s \ 5M s{ ^\ 5M V ^ S \ 6M V >®, 5M bi ~ s \ and 5M b ^ represent the 
corrections to the tree diagrams arising, respectively, from the gbb vertex diagram 
Fig. 1(c)- 1(d), the gtt vertex diagram Fig. 1(f)- 1(g), the bottom quark self-energy di- 
agram Fig.l(i), the top quark self-energy diagram Fig.l(k), the btH~ vertex dia- 
grams Figs.l(m)-l(n) and Figs.l(p)-l(q), including their corresponding counterterms 
Fig. 1(e), Fig. 1(h), Fig.l(j), Fig. 1(1), Fig.l(o), and Fig.l(r), and the box diagrams 
Figs.l(s) — 1(f). J2i $M l then represents the sum of the contributions to the Yukawa 
corrections from all the diagrams in Figs.l(c)-l(v). The explicit form of SM l can be 
expressed as 

SM l = - -^^— C^^fl-fPr, + fi>fP R + MP L + MP R + fbtP L 

+f[ 2 P^P R }u{p b )e^k), (5) 

where the C l are coefficients that depend on s, i, and the masses, and the f\ are form 
factors; both the coefficients C l and the form factors f\ are given explicitly in Appendix 
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A. The corresponding amplitude squared is 

EI M -| 2 = E|M« + M W | 2 + 2i?eB(E^)(MS S) + M^y\, (6) 



where 



ZlVC'll'W \ b ~ " l b! 

- m\p t ■ k + 2p b ■ kp b ■ p t - m 2 b p b ■ p t ) + 2m 2 b m 2 t (p b ■ k - m 2 b )] 

+ ~r ^ [( m t cot2 P + m t tan2 (3){pb ■ kp t ■ k + m\p b ■ k 

(t - mf) 2 



m t p b -Pt) + 2m b m t (p t • k - m t )\ + 



(s - m 2 b ){t - mf) 
x [(m 2 t cot 2 /3 + m 2 b tan 2 /3) (2p b • kp t ■ k + 2p b ■ kp b ■ p t - 2(p b ■ p t ) 2 
- m 2 b pfk + m 2 p b -k) + 2m 2 b m 2 (pfk-p b -k-2p b -p t )]}, (7) 



Y^5M\M^ = -—- J*f\ -C l f]h^f l t , (8) 

^ V ; 64N C x 167i 2 m 2 w {§ - m 2 ) fr[ 1 Jt ' KJ 

and 

E^(M«)t = - //% . r yT,hWfl (9) 

647V C x lQiT 2 mf v (t - mf) ~{ 

Here the color factor Nq = 3 and and /if^ are scalar functions whose explicit 
expressions are given in Appendix B. 

The cross section for the process gb — > tH~ is 

^sim^i 2 ^ (10) 

with 

imm = ^ 2V^ S_ (m t + m // -) 2 )(s - (m t -m H -) 2 ), 

and 

2 2 ^ 1 

tmax = ^ h 2 ^ ~ + m H-) )( S - ( m t ~ m H-Y)- 

The total hadronic cross section for pp — > gfe — > tiJ - can be obtained by folding the 
subprocess cross section <j with the parton luminosity: 

a{s) = ( dz^a^gb —> tH~ at s = z 2 s). (11) 

J {rat + m H - ) l\fs dZ 
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Here yfs and y/§ are the CM energies of the pp and gb states , respectively, and dL/dz 
is the parton luminosity, defined as 

Jz =2Z J z i ^rfo/piXi^fa/pi^/XiV)' ( 12 ) 

where fb/p(x, /-0 and f g / P (z 2 /x, jj) are the bottom quark and gluon parton distribution 
functions. 



2.3 Numerical results and conclusion 

In the following we present some numerical results for charged Higgs boson production 
in association with a top quark at both the Tevatron and the LHC. In our numerical 
calculations the SM parameters were taken to be mw = SOAlGeV, mz = 91.187GreV, 
m t = 176GeV, a s (mz) = 0.119, and a ew (mz) = i^g[15]. And we used the running b 
quark mass ~ 3GeV and the one-loop relations[16] from the MSSM between the Higgs 
boson masses m h H A H ± and the parameters a and f3, and chose m H ± and tan (5 as the 
two independent input parameters. And we used the CTEQ5M[17] parton distributions 
throughout the calculations. Other MSSM parameters were determined as follows: 

(i) For the parameters Mi,M 2 , and fi in the chargino and neutralino matrix, we 
put M 2 = 300GeV and then used the relation M l = (5/3)(g' 2 / g 2 )M 2 ~ 0.5M 2 [2] to 
determine Mi. We also put \x = — lOOGeV except the numerical calculations as shown 
in Fig.6(b), where /i is a variable. 

(ii) For the parameters m% ^ ^ and A t) b in squark mass matrices 



with 

m Il = m 2 Q + m 2 q + m 2 z cos 2(3{I 3 q L - e q sin 2 9 W ), 
M\ R = m 2 jj p + m 2 + m| cos 2f3e q sin 2 0w, 

M ^ M ^(t-7Z !!:?)• (14 » 

to simplify the calculation we assumed mj~ = m 2 ^ = m 2 ^ and A t = A^, and we put 
rriQ = 500GeV and A t = 200GeV. But in the numerical calculations of Fig.6(a) 
A t = A b are the variables. 

Some typical numerical calculations of the tree-level total cross sections and the 
0{a ew m 2 ^ ) /m 2 v ) SUSY electroweak corrections as the functions of the charged Higgs 
boson mass, A t = A b and /i, respectively, for three representative values of tan/5 are 
given in Figs. 3-6. 

Figures 3(a) and 4(a) show that the tree-level total cross sections as a function of 
the charged Higgs boson mass for three representative values of tan (3. For mn± = 
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200GeU the total cross sections at the Tevatron are at most only 0.7 fb and 0.1 fb 
for tan/5 = 2,30 and 10, respectively, and for m#± = 300GeU the total cross sections 
are smaller than 0.15 fb for all three values of tan/3. However, at the LHC the total 
cross sections are much larger: the order of thousands of fb for m#± in the range 100 
to 2A0GeV and tan/3 = 2 and 30; and they are hundreds of fb for the intermediate 
value tan/3 = 10. When the charged Higgs boson mass becomes heavy(< 500 GeV), 
the total cross sections still are larger than 100 fb and 10 fb for tan/3 = 2, 30 and 10, 
respectively. For low tan /3 the top quark contribution is enhanced while for high tan /3 
the bottom quark contribution becomes large. These results are smaller than ones 
given in ref. [8,9] because we used the running b quark mass ~ 3GeV in the numerical 
calculations. We have confirmed that if we chose m b = 4.5GeU, our results will agree 
with ref. [8,9]. 

In Figs. 3(b) and 4(b) we show the corrections to the total cross sections relative to 
the tree-level values as a function of m#± for tan/3 = 2, 10, and 30. For tan/3 = 2 the 
corrections decrease the total cross sections significantly, which exceed —13% for m H ± 
below 300GeU at the both Tevatron and the LHC. But the corrections decrease as m H ± 
increase. For example, as shown in Fig.4(b), the corrections range between —13% ~ 0% 
when niH± increase from 300GeU to ITeV at the LHC. For high tan/3(= 10, 30) these 
corrections become smaller, which can decrease or increase the total cross sections 
depending on tan/3, and the magnitude of the corrections are at most a few percent 
for a wide range of the charged Higgs boson mass at both the Tevatron and the LHC. 

In Fig. 5 we present the Yukawa correction from the Higgs sector and the genuine 
SUSY electroweak correction from the couplings involving the genuine SUSY parti- 
cles(the chargino, neutralino and squark) for tan/3 = 30 at the LHC, respectively. One 
can see that the Yukawa correction and the genuine SUSY electroweak correction have 
opposite signs, and thus cancel to some extent. The former decrease the total cross 
sections, which can range between —8% ~ —4% for m H ± below 300GeU, but the latter 
increase the total cross sections, which range between 10% ~ 7% for m#± in the same 
range. In such the combined effects just are about 2% ~ 3%. 

Figs. 6(a) and 6(b) give the corrections as the functions of A t = A b and \i for m H ± = 
300GeU at the LHC, respectively, assuming tan/3 = 2, 10 and 30. From Figs. 6(a) and 
6(b) one sees that the corrections increase or decrease slowly with increasing A t = A b 
and the magnitude of /i, respectively, for tan/3 = 30, 10, and the corrections are not 
very sensitive to both A t = A b and \i for tan/3 = 2, where the corrections are always 
about —12% and —13%, respectively. In general for large values of A t and small values 
of tan/3 or large values of fi and tan/3, one can get much larger corrections since 
the charged Higgs boson-stop-sbottom couplings become stronger. For tan/3 = 30, 
comparing Fig. 4(b) with Fig. 6(b), we can see that the corrections indeed become larger 
as the values of /i increase. But for tan/3 = 2 from Fig.4(a) and Fig. 6(a) we found 
that the corrections almost have no change when A t = A b become larger. Obviously 
the effects from the stronger couplings have been suppressed by the decoupling effects 
because with an increase of A t = A b all the squark masses are still heavy, which almost 
is same as discussed in Ref. [18]. 
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In conclusion, we have calculated the 0{a ew m 2 , b Jm 2 v ) and 0{a ew mt b Jrriy V ) SUSY 
electroweak corrections to the cross section for the charged Higgs boson production in 
association with a top quark at the Tevatron and the LHC. These corrections decrease 
or increase the cross section depending on tan f3 but are not very sensitive to the mass of 
the charged Higgs boson for high tan/5. At low tan/3(= 2) the corrections decrease the 
total cross sections significantly, which exceed —12% for m H ± below 300GeV at both 
the Tevatron and the LHC, but for m#± > 300GeV the corrections can become very 
small at the LHC. For high tan/3(= 10,30) these corrections can decrease or increase 
the total cross sections, and the magnitude of the corrections are at most a few percent 
at both the Tevatron and the LHC. 



2.4 Appendix A 

The coefficients C l and form factors f\ are the following: 



(jVi(s) = m b (jVi (t) = m t C s(s) _ m b 



m 2 w {s-m 2 b y mlyii-m 2 )' m 2 w {s-m 2 b ) 2 ' 

m t „v,(*\ m b mt „vi(t) rn b m t 



(js(t) = "H C V 2 (s) = ln b"H c v 2 (t) 



^{i-m 2 ) 2 ' m 2 w (s-m 2 b Y m ^(t - m t 2 )' 



C b(s) = C b(t) 



m t m b 

2 ' 

m w 



V ^[m b (g^ s) -9^ s) )-2 Pb .kg^% 

^[2( 5 r i(s) +^ i(s) )+- b (,r i(s) 



+.r i(s) ) 



2p b ■ kg 8 



Vi(sh 



V W [2( 5 r iW + 9? {S) ) + m b ( 9 ? {s) + ^ 1(S) ) + 2 Pb ■ kg™] 

. (i) [-( 5 r i(s) +^ (s) )+^(^ (s) +^ w )], 



Ms) 
J9 


= V { %^ S) + 2g^ + m b ^ is) - 








JlO 


= V { %^ S) + 2g^ + m b ^ is) - 








f VM 


= 2p b -kg^ s \ 


fV2(s) 
J2 — 


2 Pb -kg^ s \ 




fV2(s) 

J3 


= 2g\ 2(s) + 2m t cot 0(8 A^ + 5kf 


+ 8Kf) 


- 2m t g^ s) 


+ 2m b g^ s) 


fV2(s) 


= 2g^ 2(s) + 2m b tan 0(8 A^ + 8A% ] 




+ 2m b g^ 2(s) 


-2m t g V2{s) 


f V 2 (s) 

h 


1 ,V 2 (s) 

2 h > 


fV2(s) 

J8 


. 1 f V 2 (s) 

2 U ' 






= 2 Pt .kg^\ 


fV2(t) _ 

H — 


2 Pt ■ kg™ 




fV2(t) 


= 2g\ 2(t) + 2m t cot 0(8 A^ + dkf 


+ 8A^) 


- 2m tg ™ 


+ 2m b g^ t \ 


fV2(t) 

J6 


= 2gl 2(t) + 2m b tan 0(8 A^ + SA^ 


+ 5Af) + 2m b g^ t) 


-2m t g^\ 


f v 2 (t) 
h 


1 f V 2 (t) 
2 J5 > 


fV2(t) _ 

J8 - 


1 ,v 2 (t) 

2 /6 . 
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tf« = 2r ] ^ Pb .k[g^ + m b (g^+g^)], 

/ 3 S(S) = 2 V ^[m b (gf s) + g-< s) ) + 2(m 2 + p b ■ k)gf s) + (m 2 + 2p b • k)gf' ] + mgrf W ], 
£ W = 2r ? ( 1 )[m 6 ( 5l s(s) + ^ (s) ) + 2(m 2 + Pb • %* (s) + m 2 gf s) + (m 6 2 + 2p b • k)gf\ 

rs(s) 1 /•s(s) f s ( s ) _ ^ f s ( s ) 

~~ _ 2 ' ~~ ~~ 2 ' 

/i (s) = E ^ (3) {^ (2) [2m 6 (-3D 3 i2 + (1 - Ci)£>27) + mj|(A> + D 12 -D 22 

i=H ,h°,G°,A 

-D 32 ) - m 2 m b (D 23 + 2L> 39 ) - 2m fe p fe • fc(2£> 36 + D 24 + Q(D + D 12 )) 
+2m bPt ■ k(D 25 + L> 310 ) + 2m bPb ■ p t (D 26 + 2D 38 )] + r/ 1 ) [2m t (-3D 313 + (1 
+d)D 27 ) - mf(D 33 + (1 + Ci)D 23 ) + m 2 b m t (D l3 - 2D 38 + (1 + &)(A) 
-D22)) + 2m tPb ■ k(D 13 - D 3W - (1 + Ci)(-Di2 + D 24 )) + 2m tPt ■ k(2D 37 
+ (1 + d)D 25 ) + 2m tPb ■ Pt (2D 39 + (1 + C;)£>2 6 )]} 
(-k, -pb,Pt, m b, m b , mi, m t ) 
%\[2vtiyj 



sin 2(3 



E N kiN^ 3 Ri (b)Rj (t)aijD 27 (-k, - Pb , Pt , m~ b . , m~ b . , m^o , m~ tj ) , 
f Hs) = f ^ {r] W^^), Ll ^ RhNkl ^ N * kll 

f b 3 {s) = E V? ) {V W 1-4D 27 + 2m 2 b (D 22 -D - (1 - Q)(D 12 + D 22 )) 

i=H ,h°,G°,A 

+2m 2 {D 23 - (1 + d)D 26 ) + ^Pt ■ k(D 2e - D 25 )] + V ^2m t m b (l + C*)(£>22 
-D 12 - D 26 )}(-k,-p b ,p t ,m b ,m b ,mi,m t ) 

_ gVW ^2 aijl-mtN^N^Ri^Rj^D^ + m b N^ 4 Nk 3 Li (b)Lj (t) (D\ 2 
sin z. ;■ ) . . 

+D 22 ) + m^o N^ 4 N^ 3 Ri (b)Lj (t)D 12 ](-k, -p b , p t ,m~ b . , m~ b . , m^o , ), 
f Hs) = »r,W,L l ~R i ,N kl ~Nj> l ), 

/5 (S) = E V? ) {v (1) [^D 313 + 2m 2 (2D 38 -D 13 + (l-Q)(D 13 

i=H ,h°,G°,A 

+D 2& )) + 2m 2 (D 33 + (1 + d)D 23 ) + 4p fc • k(D 25 + D 310 ) - 4p t • fc(D 2 3 
+2L> 37 ) - 4p t • p 6 (£>23 + 2D39)] + n {2) 2m t m b {\ + C;)(£>i3 + £>23 
-D 26 )}(-k, -p b ,pt,m b ,m b ,mi,m t ) 

S V 2„ M r Y,<Ti j [-m t N k 4Nt 3 R i (b)R j (t)D 23 + m^N^L^L^D^ 



sin 2/? 



+L>2e) + m x^ N^N^Ri (b)Lj (t)D 13 ](-k, -p b , p t ,m~ b . , m~ b . , m^o , m^. ), 



^JW = fl {s \^~^), Ll ~R h N kl ~Ni 



klJi 



f 7 (s) = E ^ (3) {^ (1) [6(^27- £>31l) +mg(Dn- 2^2-2^22 

-2Z) 36 + (1 + Ci)(£>o + £12)) - m 2 (2L>23 + 2Z) 37 + (1 + 0)^13) - 2 Pb ■ k(D 12 
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+2D 24 + 2D U ) + 2 Pt • fc(£> 13 + 2D 25 + 2D 35 ) + 2 Pt • Pb {D 13 + 2D 26 
+D310)] +r] {2) m t m b (l + &){D 12 - D 13 - D )}(-k, -p b ,p u m b ,m bl m h m t ), 

f b s is) = 7?% (1 W 2) ), 

f b 9 {s) = E V? ) {r] {1) 2m t [-D 13 -D 26 + (l + C t )(D 12 + D 24 )] 

i=H ,h°,G°,A 

+7? (2) 2m fe [- D22 + D 24 + d(D + 2D 12 + D 24 )]}(-k, -p b ,p t ,m b ,m b ,mi,m t ) 

E ^N k4 N* k3 R t (b) Rl (t)(Di2 

i,j,k 

+D 24 )(-k,-p b ,p t ,m~ b .,m~ bi , m^o , m^), 
fi ( o S) = /{ W fa (1) ^V (2 \U^RuN kl ^Nt l ), 

/if = E V? ) W 1) 2m t [D 23 -(l + d)D 25 ] - 7] ^2m b [-D 26 + D 25 

i=H ,h°,G°,A 

+(i(Di 3 + D 25 )]}(-k, -pb,Pt, mb, "m b , m*, m*) 

+ £ ^ N k4N* k3 R t (b)R 3 (t)(Di3 

+L>25 ) ( - k , -p b , p t , m~ b . , m- b . , m^o , m fj . ) , 
f b(s) = f ^) {v (i)^ v (2) jLl ^ Rl:Nkl ^ N * i)} 

where -Do, Dij,Dij k are the four-point Feynman integrals [19]. The explicit forms of 5M Vl ^\ 5M S ^ , 5M b ^ 

can be respectively obtained from 5M Vl ( s \ 5M S ( S \ 5M b ^ by the transformation U defined 

as 

PO^Pt, S-f, k^-k, t£^T&*\ 

m t ^m b , <-> r/ (2) , X b <-> A t , m t ~ <-> mg., V^, N i3 ^ N* 4 , 

Li(b)~Li(t), Ri(b)~Ri(t), P^P L (R)^ptP R (L), -f frPL^-f fiP R . 

All other form factors /J not listed above vanish. In the above expressions we have used 
the following definitions: 

??^=m b tan/3, rj ( >=m t cotp, A 6 = ^= -, A* = —= — 

V 2m\v cos p v^H' sin p 

Li(g) = cos6> g , L 2 (g) = -sin^, i?i(g) = sin0 ? , i2 2 (?) = cos0 ? , 

(1) _ cos 2 a (1) _ sin 2 a m _ 2 (1) _ 
^° " cos 2 /?' Vh ° - cos 2 /?' Va ° - ta ^ Vg ° - ' 

(2) _ Sin 2 ^ ( 2 ) _ COS 2 ^ (2) _ 2 g (2) _ j 

(3) __ _ (3) _ sin a cos a (3) _ (3) _ 1 
n n sin p cos p u A 

Cff- = 9 ^' = 2' £ff- = 2 ^ an ^' = 2' 

n m b ^ m b n mt ^ mf 

4^ = ^n 2 (3, 4 3 1 = 1, = cot 2 0, = 1, 
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(h° = Ch° = Ch- = ~(a° = ~(g° = ~Cg- = 1, 



m w ( . 0R "ifetan (5 + m 2 cot p 

— =-(sin2/3 b - — 2 )L i (b)L j (t) 

V 2 m w 



+ m * mb (tan + cot P)Ri(b)Rj(t) - ™ b (fi - A b tan P)Ri (b)Lj (t) 
—Jp—fji-At catp)Li(b)Rj(t), 
gi l(s) = E r,f\[\-2C 2A + m 2 b {-2C ll + C l2 -C 21 + C 23 )-s{C 12 

i=H°,h°,G°,A° 

+C 23 )](-Pb, -k,mi,m b ,m b ) + [-F + F 1 + 2m 2 Gi 
-(1 +Ci)2m^G ](m^,mj,m f) )}, 

= E 2{ef ) [^-2C f 2 4 + m2C + m2(-Co-2C 11 + C 12 -C 21 + C 23 ) 

i=H~,G- 

-s(C 12 + C 23 )](- Pb ,-k,m i ,m t ,m t ) + [4 (1) (-F + F x ) - 2m 2 t (iG 



+rnl(£ ) +$ ) )(G 1 -CiG )](mlm i ,m t )} 

E{Af [^(6)1^3 1 2 (-F + Fi) + m 2 |Ay 2 (-G + d) - 2m fe m^o 



4m^ 



b 1,3 



xL J (b)R,(b)N* 3 2 G ](mlm b . ,m#>) + [-2m 6 m + \ b \ t L j (t)R j (t)VgUgG 
+X 2 R 2 {t)\V l2 \ 2 {-F Q + Fx) + m2(A 2 i? 2 (t)|^ 2 | 2 + A^L 2 (t)|?7 l2 | 2 )(-G 
+Gi)](m 2 ,m f - 2A 2 i? 2 (6)| 2 C 24 (-pb, -k,m^o,m b ,m~ b ) 



-2A 2 -R 2 (t)|Fi 2 | 2 C 2 4(-p b , -fe, m-+,m^,m£.)}, 

^Vi( s ) = ^ r ? f ) 2m 6 [(7o + 2(7ii + C' 2 i + Ci(C'o + C'ii)](-p 6 ,-fc,m i ,m 6 ,m 6 ) 

^0 ) + 2C 11 + C 21 ) + " ;r 

m 



E 4m 6 [^ 3) (C + 2C n + C 21 ) + ^|Ci(C + C u )](-Pb, -Km u m u m t ) 



i=H-,G~ 



+C2i)](-P6, -A;, m^o,m b . , m~ b . ) 



+ [m-+A 6 A t L j (t)i? j (t)y; 2 [/* 2 (C , o + C n ) - m 6 A 2 ^ 2 (0l^2| 2 (Cn 
+C 2 i)](-P6, ~k, m-+ , m t -. , m fj )}, 
= ^W (C (D ^ C (3) )F . 2 ^ ^ )iVj3 ^ iv* )jL .( 6 ) <_> A 6 L^) <- XtRjit)), 

S^ l(s) = - E r/f ) m 6 (C + Cii + CiC )(-p6,-A;,mi,m 6 ,m 6 ) 

i=ff°,fc°,G°,A° 

- E 2m fe [^ 3) (C + Cn) H ^CiCo}(-p b ,-k,mi,mt,m t ), 

i=H-,G~ m b 
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9s l{S) = E 2r lt\Ci2 + C 23 ){-p b , -k,mi,m b ,m b ) 

i=H°,h°,G°,A° 

+ E ^l 1 \c 12 + C24)(-Pb,-k,m i ,m t ,m t ) 

i=H-,G~ 

_^^{\iR 2 (b)\N i3 \\C 12 + C 23 )(-p b ,-k, m^,m- b ,m~ b .) 

ITljy ■ ■ % 3 3 



,2 

h3 



+XtRj{t) I V i2 \ z (C12 + C 23 ) (-P6, -k, m-+ , m f . , m f . )}, 
= ~ if ] ,V l2 ~ U* 2 ,N i3 ~ N^Ljib) ~ i^(&), A b L^) ~ X t R 3 (t)), 

9i 2(s) = E V? ) W 1) [-l + ^C 24 + m 2 t (C + 2Cn + UCo + C ll ) 

i=H°,h°,G°,A° 

+C 21 - C12 - C 23 ) + m\-(C 22 - C 23 ) + s(Ci 2 + C 23 )] + r/( 2 )m fe m t [CiCn 

+(1 + Ci)C ]}(-pt,-p i? -,mi,m t ,m 6 ) + Y^[ m tRi(b)Rj(t)Nk 3 N M 

sin ■ - j 

x (-Cn + Ci 2 ) + m^Lj^RiQ^N^N^Colaiji-pt, -p H - , m^o , , m t ~ . ) , 

g Ms) = g v^) {71 m^ v (2)^ Ll ^ RhNkl ^ N * i): 

9s 2(s) = E ^ (3) {^ (1) ^[Co + C 11 + C,(Co + C 12 )]+ ?? ( 2 )c l m fe C 12 } 
(-pt, -p H -,mi,m t ,m b ) 



E Ri{b)Rj(t)Nk 3 N k 4aijC 12 (-pt, -p H -,m^m~ b ., 



sin 2/3 



E V^i^i— + 4C 24 + m 2 (C + 2C U + Q(C + C n ) 



9 

i=H ,h°,G°,A 



+C 2l - C12 - C 23 ) + m^-(C 22 - C 23 ) + i(Ci 2 + C 23 )) 
+7]^m b m t [C + Ci(C + Cn)]}(-p b ,p H -,mi,m b ,m t ) 

+^f^i m t>L i (b)L j (t)N* k3 N k4 (-C 11 + Ci 2 ) 
+m i oLj(*)i2i(6)JVj 3 JV^Co]a-ij(-p6, pjy- , m^o , m~ b . , m^), 

g w) = g ^) {r] m^ v m :Ll ^ RhNkl ^ N ^ 

g 3 2(t) = E nT\n {1) rn b [C G + C 11 + aC Q + C 12 )}+rt 2 \ i m t C 12 \ 

i=H°,h°,G°,A° 
(-Pb,PH-, m i,m b ,m t ) 

+ 4 s^2/T S ^Wi t )W^ c i2(-ft»Pir-. m x t . m i i » m f j ). 
^(*) = ^%V>» r ,WL l ~R i ,N u ~Nj> l ), 
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5i S) = E m b r ] \ 1) {-QF ( Pb + k,m l ,m b ) + [C i Fo-2m 2 b (l+Ci)G 

i=H ,h°,G°,A 

+2m b G 1 ](m b ,mi,m b )} + E 2m bi |Ci-Fo(p& + k,mi,m t ) 

i=H-,G~ m b 

2 

+[-2m 2 t QG + m 2 {^ ) + ^ ) )(G 1 - CiG ) + ^FoKml^m)} 

b i,j 

-m~+\ b \ t L j (b)R j (b)V£U* 2 F (p b + k,m i .,m-+) + [ml\ 2 \N t3 \ 2 (-G + d) 
-m^X 2 b L j (b)R j (b)N* 2 (2m 2 b G - F )](m 2 ,m s .,m^o) + [mf(A 2 L 2 (t)|C/ i2 | 2 
+A 2 ^)|^ 2 |2 )( _ Go + Gi) - m-+A 6 A t L J (t)i? J (t)^*^(2m 2 G 
-F )](m 2 ,m f .,m~+)}, 
6'2 (S) = E ryf^-Fo + FO^b + Zc^^mb), 

53 (S) = E ^ (1) [Fo-i ? i-2mfG 1 + 2(l + C i )m 6 2 G ](m 2 ,m J ,m 6 ) 

i=H°,h°,G°,A° 

+ E a^C-Fo + FiJCpb + fc.mi.mtJ-^C-Fo + Fi) 

i=H-,G~ 

-2Qm 2 G + m 2 b (£ ) + # ) )(G 1 - QG )](m 2 b , mi ,m t )} 

_4m|, ^ {A 2 [jR 2 (6)!iv . 3 |2 ( _^ o + Fi) + \ N . 3 \2 m 2 { _ Go + Gl) 

m b ij 

-2m b m^oLj(b)Rj(b)N* 3 G ](rnl,m- b ^ , m^o) 

+ [A 2 i? 2 (t)|l/ l2 | 2 (-F + Fi) + m^A 2 ^ 2 ^)!^! 2 + A 2 L 2 (t)|C/ l2 | 2 )(G 1 - G ) 
-2m 6 m-+ (t)^- (i) A 6 A t V^^Go] (m 2 , m t ~ . , m -+ ) 

-X 2 b R 2 (b)\N l3 \ 2 (-F + Fi)(p 6 + fc,m^.,m^o) 

-A 2 i? 2 (t)|y i2 | 2 (-^o + Fi)(p 6 + fc, m t ~. , m-+)}, 
= ^(f W <_> ^f), y i2 - l/f 2 , N t3 <- TV*, Lj(6) <- A 6 L,(t) <- A t i?,(t)), 

= g< a \N^^N a ,V^^V a ,Ut 2 ^U i2 ), 

SA^ = (1 - cot 2 W ) [2m 2 (In ^ - 1) + m\ + m 2 - ^mfV + m 2 b F Q 

om w /i o 

+(m 2 h - m 2 - 2m 2 v )F 1 ](m 2 v ,m b ,m t ) + cot 2 6 w {-hg b v ) 2 + {g b A ) 2 

+ (sl0 2 + + [((^) 2 + (^) 2 )(2m t 2 In 4 + m 2 F Q - 2m 2 z F 1 ) 

-{{9v? ~ (^) 2 )3m 2 F ](m| ,m t ,m t ) + [{(g b v ) 2 + (^) 2 )(2m 2 ln4 
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+m 2 b F - 2m 2 z F l ) - {{g b v f - (g b A ) 2 )3m 2 b F ](m 2 z ,m b ,m b )} + -^£[(cot 2 /? 

-l)m 2 F + (m 2 - m 2 - 2m 2 cot 2 (3)F 1 + (m 2 cot 2 (3 + m 2 tan 2 
+2ml)m 2 Go — (mf cot 2 (3 + m 2 tan 2 f3)mj I -Gi](m 2 H - , mt, m^) 

+ E ^{m 2 b 4 1) [F 1 -F -2m 2 b (G + QGo-G 1 )](m 2 b ,m l ,m b ) 
i=H°,h°,G°,A° m W 

-mUP [-(1 + 2Ci)F + F 1 + 2m 2 (l + C;)G - 2m 2 t G 1 }(m 2 t , m h m t )} 

+ E ^{^[4 {1) (-Fo + F 1 )-2m 2 C,G + m 2 (4 (1) +4 (3) ) 
i=H-,G~ m W 

x(Gi - CiGoJKmg.mi.mt) - m?[-^f + CfV^o + Fi) + 2mgCiG 

m ( 

-m?(^ 2) +^ 4) )(Gi - CiG )]{r^,m i ,m b )}-2NcY l {2<T ij <T ij G 



1,3 



1 1 

mi m 



+ _ i((6)L . (t)|ii(6)i . (()(; ^ + _l_ )cos2/J 

+Ri(b)Rj(t)mtm b (tan 2 (3 — cot 2 , , rn b ), 

5kf = -2Y^{^[-^-Lj{t)Rj(t)N^(F Q - m 2 G ) + \N i4 \ 2 (R 2 (t)(-F + Fx) 

2m-+ 

-m 2 (-G + G 1 ))](m 2 t ,m i .,m i o) + [ -^-\ b X t Lj (b)Rj (b) U* 2 V* 2 (F 

-m 2 G ) + X 2 b R 2 (b)\U i2 \ 2 (-F + F ± ) - m 2 (\ 2 L 2 {b)\V i2 \ 2 + \ 2 b R 2 (b)\U l2 \ 2 ) 



x(-G + Gi)](m 2 ,m^,m-+)}, 



<5A£ j = 2j2{$l\Ni3\ 2 (R 2 (b)(-F + Fi) + m 2 b (-G + Gi)) - 2m b m^L j {b)R 3 {b) 

i,3 

xN* 2 G ](m 2 b ,m b .,m^) + [-2m-+m b A 6 A t L J (t) J R J (i)C/*y i *G 
+A 6 2 L 2 (6)|^ 2 | 2 (-F + Fi) + m 2 b (X 2 R 2 (t)\V l2 \ 2 + \ 2 b L 2 {t)\U l2 \ 2 ) 
x (-G + Gi )] (ml, m~ t . , m ~+ ) }, 

*Ag } = E 7 V{m 2 7 ? f ) [-F + F 1 -2m 2 (Go + OGo-G 1 )](m 2 ,m l ,m t ) 

-mg^hFo + F 1 - 2QF + 2m 2 (l + C*)G - 2m 2 G 1 ](m 2 m l; m b )} 
+ E -Vw 2 Kf ) (-Fo + F 1 )-2m 2 C,Go + m 2 (^+^)(G 1 



i=H-,G- m W 



-( i Go)](m 2 ,m i ,m b )-mll-^QFo + & ) (-F + F 1 ) + 2m 2 C t G 

m b 

-mli^ +ef } )(Gi - CiGoWmg.rm.mt)}, 
5A« = 5Ai 2) (f/), <5Ag } = Mi 3) ([/). 
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Here Co, C\j are the three-point Feynman integrals[19] and C24 = — 5 A + C24, while 

-q 2 y(l - y) + m\ (1 - y) + m\y 



and 



F n (q,m 1 ,m 2 ) 
G n (q,m 1 ,m 2 ) 

1 4 



f dyy n \n[- 
L dy -q 2 y(l - u) + m{( 1. - ;/) + mitj' 



y n+1 (i-y) 



9v= 2~3 SiU ° W ' 9a = 2' 



9 b v 



which are the SM couplings of the top and bottom quarks to the Z boson. The definitions of 
Qq, Uij, Vij, Nij, /j, A q can be found in ref.[2]. 

2.5 Appendix B 



no 



12 



Am t r/ (2) (2p b -k- • p 6 ) - 4m b r? (1) (p W -p t +p t -k), 

^) (r? (l) ^ ?? (2) )) 

2r? (2) (2p fe • /cp b • p 4 - m 2 ]^ • fc - 2p w • • p t ) + 2m b m t r]^ (p b • fe 
4*) ( ^(i) ^(2)^ 



- h (0 



2r/ 2 ^ (m 2 pb 
,0 



fc - 2p W • • pt) + 2m b m t m t r]^ (p t ■ k - 2p W ■ p t ), 



4^(2) ^,(1) . Wt . _ . ^ . pt _ pb . fc p (t) . pt _ 2 p 6 . fcp t . k) 

—4m b mtr)^ p0 ■ k, 

Am t rj^p b ■ k(p b ■ k - p w • p b ) - 4m b ?7 (1) p (l) ■ PbPt ■ k, 

4m t r)Wp b -k(p t -k- p W • p t ) - 4m b rj^p t • kp {l) ■ p t , 
*i?(»7 (1) ~»7 (2) ), 

(t) 



where the index i represents the two channels s and t, and p^ = p b , p K "> = pt 
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Fi gure 1: Feynman diagrams contributing to 0{oL ew 

m \b)/ m w) Yukawa corrections to 
gb — > tH~: (a) and (b) are tree level diagrams; (c) — (i>) are one-loop diagrams. The 
dashed lines represent H, h, A, ,G° and for diagrams (c) and (/); H,h,A and G° for 
diagrams (m),(p),(t) and (u); i, 6, -ff, /i, A^H^ ,G° and for (i) and (fc), where the solid 
lines represent charginos and neutralinos if the dashed lines represent squarks. For diagrams 
(d) and (g), the solid lines in the loop represent x° an d X + an d the dashed lines represent 
squarks. 
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Figure 2: Self-energy Feynman diagrams contributing to renormalization constants: The 
dashed lines represent t,b, H,h, A, ,G° and for diagram (a), where the solid lines 
represent charginos and neutralinos if the dashed lines represent squarks. 
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Figure 3: The tree- level total cross sections (a) and relative one-loop corrections (b) 
versus m H ± at the Tevatron with y/s = 2 TeV. The solid, dashed and dotted lines 
correspond to tan/? = 2, 10 and 30, respectively. 
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Figure 4: The tree-level total cross sections (a) and relative one-loop corrections 
(b) versus m H ± at the LHC with y/s = 14 TeV. The solid, dashed and dotted lines 
correspond to tan (3 = 2, 10 and 30, respectively. 
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Figure 5: The radiative correction from top, bottom quarks (dashed line) and genuine 
SUSY particles (dotted line), as well as total contributions (solid line) when tan (3 = 30 
at the LHC with y/s = U TeV. 
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Figure 6: Relative one-loop corrections versus A t , (a) as well as /i (b) at the LHC 
with y/s = 14 TeV, where m H ± = SOOGeV" and the solid, dashed and dotted lines 
correspond to tan (3 = 2, 10 and 30, respectively. For (a), \x = — lOOGeV, and for (b), 
A t = A b = 200GeV. 
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3 Supersymmetric Electroweak Corrections to W ± H T 
Associated Production at the CERN Large Hadron 
Collider 

ABSTRACT 

The 0(a ew m^ b ^/m'^ v ) and 0(Q eBJ mtn/m^) supersymmetric electroweak correc- 
tions to the cross section for W^H^ associated production at the LHC are cal- 
culated in the minimal supersymmetric standard model. Those corrections arise 
from the quantum effects which are induced by the Yukawa couplings from the 
Higgs sector and the chargino-top(bottom)-sbottom(stop) couplings, ncutralino- 
top(bottom)-stop(sbottom) couplings and charged Higgs-stop-sbottom couplings. 
The numerical results show that the Yukawa corrections arising from the Higgs 
sector can decrease the total cross sections significantly for low tan/?(= 1.5 and 2) 
when m H +(< 300)GeV, which exceed —12%. For high tan/3 the Yukawa correc- 
tions become negligibly small. The genuine supersymmetric electroweak corrections 
can increase or decrease the total cross sections depending on the supersymmetric 
parameters, which can exceed —25% for the favorable supersymmetric parameter 
values. We also show that the genuine supersymmetric electroweak corrections de- 
pend strongly on the choice of tan/3, A t , Mq and /i. For large values of A tl or large 
values of fi and tan /3, one can get much larger corrections. The corrections can 
become very small, in contrast, for larger values of Mq. 

3.1 Introduction 

One of the most important objectives of the CERN Large Hadron Collider (LHC) is 
the search for Higgs boson. In various extensions of the Higgs sector of the standard 
model(SM), for example, in the two-Higgs-doublet models(THDM) [1], particularly the 
minimal supersymmetric standard model(MSSM)[2], there are physical charged Higgs 
bosons, which do not belong to the spectrum of the SM and therefore their discovery 
would be instant evidence of new physics. In much of the parameter space preferred by 
the MSSM, namely m H ± > m w and 1 < tan/3 < m t /mb[3,A], the LHC will provide the 
greatest opportunity for the discovery of charged Higgs boson. Previous studies have 
shown that for a relatively light charged Higgs boson, m#± < m t — rrib, the dominate 
production processes at the LHC are gg — > ti and qq ti followed by the decay 
sequence t — > bH + — > 6t + z/ t [5], and for a heavier charged Higgs boson the dominate 
production process is gb — > tH~ [6,7,8]. Besides the processes mentioned above, in 
Ref.[9] Dicus et al. also studied the production of a charged Higgs boson in association 
with a W boson via bb annihilation at the tree level and gg fusion at one loop at 



3 Supersymmetric Electroweak Corrections to W ± H T Associated Production at the 
CERN Large Hadron Collider 32 



hadron colliders. Since the leptonic decays of W boson would serve as a spectacular 
trigger for the charged Higgs boson search, these processes seem attractive. But the 
authors of Ref . [9] only considered the case where the value of tan f3 to be in the range 
0.3 — 2.3. Recently Barrientos Bendezu and Kniehl[10] further studied these processes 
and presented theoretical predictions for the W ± H Zf production cross section at the 
LHC and Tevatron's Run II, where they generalize the analysis of Ref. [9] for arbitrary 
values of tan/3 and to update it. They found that the W ± H Zf production would have 
a sizeable cross section and its signal should have a significant rate at the LHC unless 
m#=F is very large. 

As analyzed in Ref. [7, 11], the search for heavy charged Higgs bosons with m H + > 
m t + rrib at a hadron collider is seriously complicated by QCD backgrounds. For 
example, the processes suggested in Ref. [10] suffer from the irreducible background due 
to top quark pair production, qq — > tt and gg — > tt with subsequent decay through the 
intermediate state bbW + W~ , and heavy charged Higgs boson produced in association 
with W ± gauge bosons cannot be resolved at the LHC, via semileptonic W + W~ decays, 
for charged Higgs boson masses in the range between 2m t and 600GeV at neither low 
nor high tan/3[ll]. However, recent analyses[12,13] have shown that the decay mode 
H + — > t + u, indeed dominant for light charged Higgs bosons below the top threshold 
for any accessible tan /3 [14], provides an excellent signature for a heavy charged Higgs 
boson in searches at the LHC. The discover region for H ± is far greater than had 
been thought for a large range of the (m H ±,tan (3) parameter space, extending beyond 
m H ± ~ ITeV and down to at least tan/3 ~ 3, and potentially to tan/3 ~ 1.5, assuming 
the latest results for the SM parameters and parton distribution functions as well as 
using kinematic selection techniques and the tau polarization analysis [13]. Of course, 
it is just a theoretical analysis and no experimental simulation has been performed to 
make the statement very reliable so far. 

Since the contributions to the W ± H T production cross section due to bb annihila- 
tion at the tree level are greater than ones due to gg fusion which proceeds at one-loop, 
it is important to calculate the one-loop radiative corrections to the W ± H^ production 
via bb annihilation for more accurate theoretical predictions for the cross sections. In 
this paper we present the calculations of the O^e^m^/m^) and O^a^m^^jm^) 
supersymmetric(SUSY) electroweak(EW) corrections to this W ± H T associated pro- 
duction process at the LHC in the MSSM. These corrections arise from the quantum 
effects which are induced by potentially large Yukawa couplings from the Higgs sec- 
tor and the chargino-top(bottom)-sbottom(stop) couplings, neutralino- top(bottom)- 
stop(sbottom) couplings and charged Higgs-stop-sbottom couplings which will con- 
tribute at the 0{a ew m\(^ ) lm\ v ) to the self-energy of the charged Higgs boson. The 
relevant QCD corrections are expected to be larger, but not yet available. 

The arrangement of this paper is as follows. In Sec. II we give the analytic results. 
In Sec. Ill we present some numerical examples and discuss the implications of our 
results. Some notations used in this paper and the lengthy expressions of the form 
factors are summarized in Appendix A, B. 
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3.2 Calculations and formulas 

The Feynman diagrams for the charged Higgs boson production via b(pi)b(p 2 ) —* 
W ± {k)H T (p 3 ), which include the SUSY EW corrections to the process, are shown 
in Fig.l and Fig.2. We carried out the calculation in the t'Hooft-Feynman gauge and 
used dimensional reduction, which preserves supersymmetry, for regularization of the 
ultraviolet divergences in the virtual loop corrections using the on-mass-shell renor- 
malization scheme[15], in which the fine-structure constant a ew and physical masses 
are chosen to be the renormalized parameters, and finite parts of the counterterms 
are fixed by the renormalization conditions. The coupling constant g is related to the 
input parameters e, mw, and mz via g 2 = e 2 /s 2 u and = 1 — m w /m 2 z . As far as the 
parameters (3 and a, for the MSSM we are considering, they have to be renormalized, 
too. In the MSSM they are not independent. Nevertheless, we follow the approach of 
Mendez and Pomarol[16] in which they consider them as independent renormalized pa- 
rameters and fixed the corresponding renormalization constants by a renormalization 
condition that the on-mass-shell H + lui and hll couplings keep the forms of Eq.(3) of 
Ref. [16] to all order of perturbation theory. 
We define the Mandelstam variables as 

s = (P1+P2) 2 = (k + Ps) 2 , 
i= (pi - kf = (P2-P3) 2 , 

u = (pi - P3) 2 = (P2 - k) 2 . (1) 
The relevant renormalization constants are defined as 

rn^o = m w + $ m W, m zo = m z + S m z, 
tan p — (1 + 8Zp) tan f3, 

sinao = (1 + SZ a ) sin a, 

Wt = (1 + 5Z w f' 2 W^ + iZ]it w± d^H^, 

H^il + dZH^ 1 ' 2 ^, 

Z% = (1 + 5Z z f 2 Z» + iZ^A, 

A Q = {l + 5Z A fl 2 A, 

H = {l + 5Z H f' 2 H + z]l 2 h h, 

h = (l + 5Z h y/ 2 h + Zl /2 H. (2) 

Taking into account the 0{a ew m 2 ^ j /m 2 v ) and 0(a ew m^/m^ v ) SUSY EW correc- 
tions, the renormalized amplitude for bb — * W~H + can be written as 

M ren = M (s) + M (t) + [5M Vl{s) + 5M s(s) + 5M V2{s) ](Hi) + [5M Vl{s) 

+5M s(s) + 5M V2{s) ](A) + 5M Vl{t) + 5M s{t) + 5M V2{t) + 5M box , (3) 

where M (s) and M (t) are the tree-level amplitudes arising from Fig. 1(a) and Fig. 1(6), 
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respectively, which are given by 

MS S) = -*E ff Myi 2 \ I>i+ ^ 12 2N (M 1 -M 2 + M 3 -M 4 ) (4) 



and 



M (t) = -=-^ — -(2h b p 12 M 2 - h b m b [3 12 M 5 + h t m t p u M 6 - h b (3 l2 M l2 ). (5) 
V2(t - ) 

Here h b = gmj, / v^^w/ cos /3 and /i t = gm t / \^2mw sin /3 are the Yukawa couplings from 
the bottom and top quarks, pi and p 2 denote the momentum of incoming quarks b and 
b, respectively, while k and p^ are used for the outgoing W~ Boson and H + Boson, 
respectively. The notations q;^-, and ipij used in the above expressions are defined 
in Appendix A, and Hi stands for Higgs Bosons h with % — 1 and H with % = 2. M« 
are the standard matrix elements, which are defined by 

M 1 =v(p 2 )P R u(p 1 )p 1 -e(k), 
M 2 =v(p2)PLu(pi)pi-e(k), 
M 3 = u(p2)-Pflw(pi)p2 • £(*;), 
M 4 = v(p 2 )P L u(pi)p 2 ■ e(k), 
M 5 =v{p 2 ) 4{k)P R u{ Pl ), 
M 6 =v{p 2 ) 4(k)P L u( Pl ), 
M 7 = v{p 2 ) Mk)P R u(pi)pi ■ e(k), 
M 8 = v(p 2 ) M k ) p Lu(pi)pi ■ e(h), 
M 9 = v(p 2 ) Mk)P R u(pi)p 2 ■ e(k), 
M w = v{p 2 ) jt{k)P L u{ Pl )p 2 ■ e(k), 
M u = v(p 2 ) jk Mk)P R u( Pl ), 

M 12 = v(p 2 ) jk /t{k)P L u{ Pl ), (6) 

where Pl,r = (l=F75)/2. The vertex and self-energy corrections to the tree-level process 
are included in 5M V ' S , which are given by 



sin((3-a)sma 1/2 cos(/3 - a) 

H ^ ( tan aoZ a + Z hH ) 5 — (smadZ a 

s — m H " s — ml 

-cosaZ 1 I { 2 h )}j2M j + 6M v ^\H), 



V2(s — mjj ^6 2 2 2 
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+ ^s^ Z ^ " ^ + ^ - M 4 ) + 6M™{A), 
5M^m = ^ Ej^y^< - (s - mlJ5Z Hi - (s 

-m 2 H )Z]i 2 h — (s — m\)Z)£\ £ Mj + 5M S ^(H), 

j'=i 

5 MS(s ){A) = igh^P [5m 2 A m 2 A)5 z A ]( Ml -M 2 + M 3 - M 4 ) 

V 2(s — m A ) 

+5M s(s \A), 

cos a cos(/9 — a) „ cr7 cr7 i/ 2 
(sm p cos pdZp — tanaoZ a — Z h ' H 



s — m 



H 

i/2 v sin a sin(/? - a) . 

+m H // f j H ,) H : g [sm (J cos pdZfs 

s — m h 

- ttma5Z a + ZH + m w Z]&)} ]T M j + 5M v ^ s \H), 

SM^\A) = -^^[^ + W + hz B+ 
V2(s - m A ) g 2 2 

+ \dZ w -](M l - M 2 + M 3 - M 4 ) + 5M y2(s) (A), 

5M V ^ = ! 9 (2h b (3 12 M 2 - h b m b (3 12 M 5 + h t m t f3 u M, 
\/2{t-mi) 

-h b (3 12 M 12 ){^ + l -5Z{ + \bZ h L + l -5Z w -) + 5M^\ 
= 19 2 J (2mf-^ + mlSZl - t5Z{){2h b p 12 M 2 

\ Z[t — 777- ^ J 

1 1 (5 

-/i 6 m 6/ 3 12 M 5 - h t (3 12 M l2 + -h t m t (3 u M 6 ) + -(2t— 

2 2 m t 

+m 2 t 5Z t R - tSZ^htrritPuMe] + 5M s{t \ 
6 &™ = ^f^yK cot^ - cos^Z, + 1^ + \sz R 

+\SZ H+ + ^Lz^)M 6 + m 6 tan/3(^ + sin 2 /WZ, + \bZ\ 

+ l -8Z h R + ^Z ff+ - ^-Z)l 2 w ){2M 2 - M 12 - m b M 5 )] + 5M^\ (7) 
2 2 tan p 
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with 



5g 5e 1 Sm 2 z 1 5m 2 z — Sm^ 
g e 2 m| 2 m|- 



5h b 5g 5m b 1 5m 2 w 2 

— = — + + COS f 362 ^ 

h b g m b 2 r?% 



5m t 1 5m 2 w . 2 

— = — + -— sin (35Z/3, 

h t g m t 2 



A7 5g lSntfy 1 rrny„ 



rl/2 

(? ' 2 m 2 ^ 2""" ' tan/5 
<SZ a = + 1^ - ±*Z fc - cot aZ]jl - sin 2 (38Z p . (8) 

Z TTljy Z 

The <5e/e appearing in Eq.(8) does not contain the 0(a ew m 2 ^/m 2 v ) corrections and 
needs not be considered in our calculations. And 5M Vl ^(Hi), 5M Vl(s \A), 5M s ^(Hi), 
5M s(s \A), 5M v ^(Hi), 6M V *W(A), 6M Vl ®, 5M S ^, 5M V ^ and 5M box represent the 
irreducible corrections arising, respectively, from the bbH(h) vertex diagrams shown in 
Fig. 1(c) — 1(d), the bbA vertex diagrams shown in Fig. 1(c) — 1(d), the H and h boson 
self-energy diagrams in Fig.l(i) — l(k), the A boson self-energy diagrams shown in 
Fig.l(i) — l(fc), the H(h)W~H + vertex diagrams shown in Fig.l(/) — l(/i), the AW~H + 
vertex diagrams shown in Fig.l(/) — 1(h), the btW~ vertex diagrams Fig.l(Z) — l(o), 
the top quark self-energy diagrams Fig.l(r), the tbH + vertex diagrams Fig.l(p) — l(q), 
and the box diagrams Fig.l(s) — l(x). All above 5M V ' S and 5M box can be written in 
the form 

5M v,s,box = • ^ fY^Mi, (9) 
i=i 

where the fY' s,box are form factors, which are given explicitly in Appendix B. 

Calculating the self-energy diagrams in Fig.2, we can get the explicit expressions 
of all the renormalization constants as following: 

- E 32^^ 2m *^ + h l m ^ii)Bf Ht + 2h b h t (3 ll (32 l B t bH >] 

+ E a^-^l^l 2 ^ + ^) + m^eUN% + n;!)b^} 
+E 5^W«i) 2 l^il 2 + ^(^ 2 ) 2 I^I 2 ](^ + + 
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+ E T^M^i' + h 2 b P 2 2i )B tbH ? + 2h b h t m b [3 li p 2l B , tbH >} 



-16tt 2 



+h b h t m^9 b il et 2 (u j2 v j2 + u; 2 v; 2 )b'^} : 

w = 16^2**-' 6 ~~ *^ + 2mJ ° ' ~ * 

VtK 2 - m 2 ) 2 + (m 2 b + m 2 )m 2 w }BD, 



2m w 



S7 ^ f ( m fc m t) 2 , 7 v;/./ pU7</\ r/,,,2 ,,,2\ 

5Zw ~ 32^<; { (5 ° ) + [^- m *) 



■06t D W6t\ , iy2 _2\2 



+(m 2 + m 2 )m^]^ 6t }, 

,2„2 ™2 



5m z = T^^ 3 - 2^)(^ 66 + iC) - m?(3 - 4f w )(B*« - B™)] 

LOC w 7I Z 

+ o^KW* " 25$*) " + 2B§*)], 

-^/^m^n/^i^ + E (^) 2 (^) 2 (^yi + ^yi) 2 ^^], 

= i|pK-2/M fc [m?(l + T + 2< fc ") + m^Sf* tt ] - 2Z£a|K(l 
+^ + 25j?*») + m^Sf + £ [(M^V ^O 2 < fe ^ 
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Hh^Ql^fB^] + £ hlmjal k (l + B^] 

i 



5Z Hk = J— 2 {2h 2 t a\ k {B^ tt + 2m 2 t B' Hktt + m^B'***) + 2h 2 b a 2 2k (B^ bb 
+2m 2 b B'^ bb + m 2 Hk B'^ bb ) + £ [M^e],^^ 

+{h b 9 b ii ,9 b jj , Q 2 lj/k ) 2 B Hkb ' b3 ] } , 
5m 2 Ak = ^{2h 2 t p 2 lk [m 2 (l + B™) + m\B^\ + 2/#&[mg(l + 

+m^»] - E [(^^e? YJfc ) 2 s ^ + (M^^e^ fc ) 2 s ^] 
+ E ^K^( 1 + ^) + E ^H/^ 1 + 

i i 

8Z Ak = ^mtPUB^ + m\ k B^) + 2h 2 b (3 2 2k (B?* bb + m\ B^ bb ) 
- E \( h tO\ i :d] j ,Ql rk ) 2 B^ Lktlt: ' + (h b 9 b iil 9 b jj/ Q^ j , k ) 2 B Akb ' bj ]}, 

Z H+W = - -[(h t m t p u + h b m b [3 12 )((m 2 b - m 2 )(B° bt - B» +bt ) - m 2 H+ Bf bt ) 



16V2n 2 m 2 H 



.m 



w 



+ E 9 b l 9 b i ,9 t jl 9) j ,{h b el j , l + ^e? Y1 )(m|. - m|)(S oS< ^ - Bf^')], 



+ 32^4 4 E. {M^^'©Jvi[(3 - 2^)6^ - 2^^ 2 ^ 2 ]K - 
-B ^) - M^e* vi [(3 - 4^)^ - 4 S 2 ^. 2 ](m| - m|)(B °^ - 

z - = i 6 ^K-U) [2mg(1 + + 25r) " 2m?(1 + + 25r) 

-mUB™ - B™)\ 

E [(hb9 b i/ 9 b jj ,) 2 Q 2 , j , 1 Q 2 , j , 2 B^ btb] + (/i t 6 l - i ,6'* j ,) 2 6-, j , 1 0- Y2J B ( f* l<J ] 



3an«i2 



167T 2 (m 2 - .j^., 



167r 2 (m 2 — r?%) 



+ ^) + h 2 t m\{\ + B ^)], 



^Hfc - Z hH\h~H, (10) 
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with 

^ = (_ir 1 d Vy(y-lH^(l-y)+mIy- (12) 

The notations 0|- and 0*- used in above expressions are defined in Appendix A. Ai 
stands for A with % — 1 and G° with i = 2. Hf stands for H + with % = 1 and G + 
with i = 2. 5Z b L , 5Z b R can be obtained, respectively, from <5Z£, <5Z*j by the 
transformation: 

h b <-> h t ,m b <-> m t , rag. <-> ra t ~, <-> a 2 i, ft* <-> /?2i, ^ <-> 0^-, iVi4 -> iV i3 , C/ i2 -> V^. 
The corresponding amplitude squared is 

£|M ren | 2 = ]T|M (s) + M<*>| 2 + 2i?e£[(]T5M)(M (s) + M (t) )t]. (13) 

The cross section for the process bb — > W ± H T is 

1 



167TS 2 



S|M ren |^t (14) 



with 



4 = ^ + ™ g ~ S ± \^{s - {m w + m H -f){s - {m w - m H -)*). (15) 

The total hadronic cross section for pp — > 66 — > W ± H T can be obtained by folding the 
subprocess cross section <j with the parton luminosity: 

a(s) = f 1 dz^a(bb -> W^if* at s = ^ 2 s). (16) 

Here and v^S are the CM energies of the pp and bb states , respectively, and dL/dz 
is the parton luminosity, defined as 



2z —f b/P (x,fj,)f- b/P (z 2 /x,fj), (17) 

J z X 



dL „ f 1 dx 

dz 



where f b /p(x,ii) and f b / P (z 2 /x, /i) are the bottom and anti-bottom quark parton dis- 
tribution functions, respectively. 
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3.3 Numerical results and conclusion 

We now present some numerical results for the SUSY EW corrections to W ± H T as- 
sociated production at the LHC. The SM input parameters in our calculations were 
taken to be a ew (m z ) = 1/128.8, m w = 80.375GeV and m z = 91.1867GeV[17], and 
m t = 175.6GeV and m& = 4.7GeV, which were taken according to Ref.[10] for compar- 
ison. We used the CTEQ5M parton distributions throughout the calculations [18]. The 
one-loop relations[19] between the Higgs boson masses M^h^h^ an d the parameters 
a and (5 in the MSSM were used, and m H + and (5 were chosen as the two independent 
input parameters. Other MSSM parameters were determined as follows: 

(i) For the parameters Mi, M2 and \i in the chargino and neutralino matrix, we take 
M 2 and \i as the input parameters, and then used the relation M\ = (5/3)(</ 2 /g 2 )M 2 — 
0.5M 2 [2] to determine M x . 

(ii) For the parameters m|~ q ^ and A t b in squark mass matrices 

with 

M\ L =m 2 Q + m 2 q + m 2 z cos 2(3( l\ L - e q sin 2 6 W ) , 
M RR = rrtjj jj + m + m z cos2pe g sin 0w, 

{T-m «=§)■ (19) 

to simplify the calculation we assumed Mq = Mjj = M^ and A t = A b , and we used Mq 
and A t as the input parameters except the numerical calculations as shown in Fig.6, 
where we took , m~ b and A t = A^ as the input parameters. 

Some typical numerical calculations of the Yukawa corrections and the genuine 
SUSY EW corrections are given in Fig.3-4 and Fig. 5-9, respectively. 

In Fig. 3 we present the Yukawa corrections to the total cross sections relative to 
the tree-level values as a function of for tan (3 = 1.5, 2, 6 and 30. For tan /3 — 1.5 
and 2 the corrections decrease the total cross sections significantly, which exceed —6% 
for m H + < 500GeV and —12% for m H + < 300GeV. For tan/3(= 6) these corrections 
also decrease the total cross sections, although relatively smaller, which exceed —2.5% 
for m H + < 500GeV and exceed -5% for m H + < 250GeV. But for high tan /3(= 30) 
these corrections become positive, which increase the total cross sections slightly. Note 
that there are the peaks at m#+ = 180.3GeV, which arise from the singularity of 
the charged Higgs boson wavefunction renormalization constant at the threshold point 
m H + —m t + m b . 

In Fig. 4 we show the Yukawa corrections as a function of tan f3 ioim H + = 100, 150, 200 
and 300GeV. For tan (3 < 4 the corrections reduce the total cross sections by more than 
10% with m H + = 100, 150 and 200GeV. With m H + = 300GeV the corrections are only 
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significant for 1 < tan/3 < 5. For high tan/3(> 10) the corrections become negligibly 
small for all above m#+ values. 

Fig. 5 gives the genuine SUSY EW corrections as a function of for tan/5 = 
1.5,2,6 and 30, respectively, assuming M 2 = 300GeV, \x = — lOOGeV, A t = A b = 
200GeV, and Mq = Mfj = Mfj = 500GeV. From this figure one sees that the correc- 
tions are very small and negligible, which is reasonable because the squark masses are 
now very large and also the couplings of the charged Higgs boson-squarks are small for 
the values of A ttb , Mqfj jj and /i used in those numerical calculations. In contrast, in 
Fig. 6 when we take the lighter sqarks masses: = lOOGeV and m~ bi = 150GeV, and 
put A t = A b = ITeV, which are relatively larger, assuming M 2 = 200GeV, ji = lOOGeV 
and Mq = Mfj, the genuine SUSY EW corrections are enhanced significantly, espe- 
cially for low tan/3(= 1.5) and m H + below 250GeV, which can exceed —30%. But 
when rriH+ > 250GeV the corrections increase the cross sections, which can exceed 
10%. For tan /3 — 6 and 30 the corrections are at most 10% and become small with an 
increase of m H + . The sharp dips at m H + = 250GeV are again due to the singularity of 
the charged Higgs boson wavefunction renormalization constant at the threshold point 
m H + = + m~ hi = 250GeV. 

Fig. 7, Fig.8 and Fig.9 give the genuine SUSY EW corrections versus A t = Ab, 
Mq = Mfj = M^ and fx, respectively, for tan/5 = 1.5 and 30. In each figure we fixed 
m H+ = 200GeV and M 2 = 300GeV. 

Fig. 7 shows that the corrections are negative for tan (3 = 1.5 and positive for tan (3 = 
30, assuming Mq = Mfj = M b = 400GeV and fi = lOOGeV. For both tan/3 = 1.5 
and 30 the magnitude of the corrections increases with increasing A t = A^. When 
A t = Ab = ITeV the corrections can reach —6% and 7.5% for tan/3 = 1.5 and 30, 
respectively. Otherwise, when A t = A b decrease to lOOGeV, the corrections become 
negligibly small. This result is due to the fact that large values of A t = A b not only 
enhance the couplings, but also give a large splitting between the masses of ti(bi) and 
£2(^2), and in consequence lighter t\ and b\. 

Fig.8 also show that the corrections are negative for tan/3 = 1.5 and positive for 
tan/3 = 30, assuming A t = A b = 500GeV and fi = lOOGeV. When = 250GeV 

the corrections can reach —3.6% for tan/3 = 1.5 and 7.3% for tan/3 = 30. But the 
magnitude of the corrections drops below one percent when Mq ^ jj increase to 750GeV. 
This is because for larger values of Mq ^ ^ the squarks have larger masses and their 
virtual effects decrease due to the decoupling effects. 

In Fig.9 we present the genuine SUSY EW corrections as a function of fi, assuming 
A t = A b = 500GeV and M Q = M = M b = 400GeV. For tan (3 = 30 the magnitude of 
the corrections increase with an increase of which varies from 0% to 5% when \fj,\ 
ranges between ~ 500GeV. For tan/3 = 1.5 the corrections are relatively small and 
increase slowly from about 0% to 3.5% when fi ranges between — 500GeV~ 500GeV. 
This result indicates that large values of ji and tan/3 can enhance the corrections 
significantly since the couplings become stronger. 

In conclusion, we have calculated the 0(a ew m^/m^ v ) and 0{a ew m^/m\ ¥ ) SUSY 
EW corrections to the cross sections for W ± H Zf associated production at the LHC in 
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the MSSM. The numerical results show that the Yukawa corrections arising from the 
Higgs sector can decrease the total cross sections significantly for low tan/?(= 1.5 and 
2) when m#+(< 300)GeV, which exceed —12%. For high tan/5 the Yukawa corrections 
become negligibly small. The genuine SUSY EW corrections can increase or decrease 
the total cross sections depending on the SUSY parameters, which can exceed —25% 
for the favorable SUSY parameter values. We also show that the genuine SUSY EW 
corrections depend strongly on the choice of tan/3, A t , Mq and /i. For large values of 
A t , or large values of /i and tan/3, one can get much larger corrections. The correcan 
become very small, in contrast, for larger values of Mq. 

3.4 Appendix A 

We present some notations used in this paper here. We introduce an angle ip = (3 — a, and 
for each angle a, (3, tp, 9 l or 9 b , we define 

(cos a sin a \ „ I cos (3 sin (3 \ I cos ip sin ip \ 

— sin a cos a J ' y I — sin cos f3 J ' ^ l— sin <p cos ip J ' 

at _ ( cos #* sin61 * \ R b _ ( cos9 b sm6 b \ 
lJ \ — sin 9 l cos 

We define six matrix 9* fci ,i = 1 — 6 for the couplings between squarks and Higgses: 



2m t cos a A t cos a + [i sin a 
A t cos a + \x sin a 2mt cos a 



el 



2mt sin a A t sin a — // cos a 



u2 \/2 \ A t sin a — ,u cos a 2mt sin a 

q2 -1 I 2rrib sin a Af, sin a + cos a 

ul v 7 ^ V ^4fe sin a + cos a 2m5 sin a 



e 2 . 




2m;, cos a A?, cos a — u, sin a 



lj2 \/2 \ Af, cos a — sin a 2mb cos a 



1 / A t cos /3 + ^ sin /3 

s/2 \ —M cos /3 — [i sin /? 

3 I I A t sin [3 — u, cos /3 

ij2 ~ 72 V -^*sin/3 + /icos/? 

1 / A;, sin /3 + cos /3 



6 



1 



^2 V -A 6 sin/3-^cos/3 

4 1 / — Aj, cos /3 + ^ sin (3 

ij2 ~ 7^ V A bCosp - fism(3 



a 5 , 



nib sin /? 
A;, sin (3 + u, cos (3 nit sin /3 

— nibcosfl \ 



ij2 ~ I -A b cos (3 + [i sin /3 
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„ 6 _ ( m t cos (3 A t cos (3 + n sin (3 \ 
\ m b cos/3 J 

6 / m t sin (3 A t sin (3 — n cos (3 \ 

3.5 Appendix B 

The form factors defined in Eq.(9) are the following: 

+(4m 2 . + s)Ci](s,mg,mg,mg,mg,m^.)} 

, ~9^ l b^2i a 2j ( Pjl rr>bbAi i , 2 /- 2 2 2 2 2 m 

+ L QO /» V 2 J ^o - (4m 6 -s)Ci(s,m 6 ,m 6 ,m 6 ,m 6 ,m A .)] 

+2h b h t m 2 t (3ii[32i + hlm b m t (3 2i )Co + (2h 2 m b m t p' 2 i + h b h t s(3iif3 2 i 
+2hlm b m t pl)C 1 ](s, m 2 , m 2 , m 2 , mf, m^)} 

+ E E ^^4? foW^xl^^o + d + C 2 ) 
.^-i .^T 1 . 16it z (s — mir ) J 

i,3,k,li',]' v ■ H i / 

+^+^^^2^it / fe2V fc2 Co - h 2 t m b 0\ 2 0) 2 |y i2 | 2 Ci](s, m 2 ,m 2 ,m|,m|,m~+) 

, \- \- 9hl<Piie b jf e b kk/ N l3 e 2 fk/i 

hi h l^ 2 {s-m 2 H ) + C i + ^) 

i,3,k,lj',k' y H i> 

-m b e b j2 0k2 N k c i + m 5i o6 b jl e b k2 N lz CQ]{s, m 2 b ,m 2 b ,m\^m\ k ,m^), 

f 2 Vl{s) (H) = /^(^(Mnl <- Mn2,#nl ^ ^2,^2 - K*2,^n3 - K 3 ), 

fi l{s) (A)=f^\A) a + ^ {s \A) b , 



where 



A™ W« = E E ^7^y [-^^il^l 2 ( g o + d + d) 

i,j,ki',j' 

-m^h b h t Q\ 2 e) x \J k2 V k2 CQ + /ii 2 m 6 ^ 1 ^ 2 |l/ j 2| 2 C 2 ](s, m 2 , m 2 , m|, mf., m~+} 

+ £ E ^1r^tr [-^^(Ci + c 2 + c 3 ) 

i,j,k%',j' v A/ 

/ 2 Vl(s) (A) a = /^(^(/l^l - M<2,0nl - ^2^2 - F n * 2 , iV n3 ~ AT* 3 ), 
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fz l{s \A) a = f? (s) (A) a , 
fI l{s) (A) a = f^\A) a , 

h {A)b -\z2V2^{s-m\) {B ° - [AmbC ° 
+(4m 2 - s)C'i](s,mg,mg,mg,m^m^.)} 

+ 1^ oo m JT^o - \ Am b - ■VC 1 (s,m b ,m b ,m b ,m b ,m A .)} 

i 32v27r 2 (s — m A ) 

+ [(h 2 m b m t (3f i - 2h b h t ml(3iiP2i + h 2 rn b mtPli)Co 
-h b h t {^m 2 b - s)PiiP2iCi](s,m b ,m b ,mf,m'f,m 2 H+ )}, 

i 

f2 Vl(s \A) b = -fI l(s) (A) b = f?<'\A) h = -f? {s) (A) b ; 

ft (S \H) = 2 M Ul + O + (2mW + IB? 6 )] 

i,j V ^ S m H x )V s m Hj) 

y- y- 3^6 Q2fc (gn (6>^-> ) 2 (g|Q^Q|yfcQ|vj j^. 

ilkife 16v^(S-m^)(5-m^) 

v v ^6^«2^ fc i(^0 2 (^V) 2 e^, / 9),, fc j?gf - if - j 



8\/27r 2 (s-m 



8\/27r 2 (s - '/// 



i,j i',y 



. M2 K(i + C) + ^f] 

y- y- 9 h bfcl( d jj>) 2 ( e ii>) 2 (®i>j>l) 2 gsbjbj 

16V2TT 2 (s-m 2 A ) 2 



9hbh 2 t f32 1 (9y 2 (et l ,) 2 (ef fl r sut 



y^ y^ a"-o"-t^n" J3 v W \~i'j'U ^ 



2 \2 
A) 



* 3 



/ 2 * ( *>(A) = -/f >(A) = /f>(.4) = -^<*>(A); 
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+h b m 2 (3 2 i)C + h b mwP2iCi - (h b m b (3 2 i + h b m 2 (3 21 
-2h t m b m t Pu)C 2 - 2h b f3 21 C o - h b /3 21 (t + u- 2m b )C 12 
-2h b m 2 H +(3 2 iC 22 ] (s, m 2 H+ , m 2 ^, m 2 , m 2 bl m 2 t )} 
gh b h t aua 2 i f 3 m 2 

+ ^ 16V2, 2 (s- m 2 H J { 2 ktf3nBo + ^ m ^ 

-h b m b m t (3 21 )C + htrri^Pu^ - {h t m 2 b [3u + h t m 2 t f3u 
-2h b m b m t p 2 i)C 2 - 2h t (3i 1 C o - h t f3u(t + u- 2ml)C 12 
-2h tm 2 H +{3 xx C 22 \(s, m 2 H+ , m 2 w , m 2 , m 2 , m 2 )} 

+ J2 16vr 2 f? - m 2 ) a2 ^^') 2 ^i^i'^i^' ^(^ i'fc'i 

i,j,k,li',j',k' ^ Hi) 

+ht&i'k'i)C 2 (s, m 2 H+ , mly^m 2 ^ , m~. , m| ) 

+ Yj H 16vr 2 fs -m 2 T a2 '^«') 2e ii e iJ' e *i^ e i'< , '^ 6e *'*'i 
i,j,k,li',j',k' \ Hi' 

+h t @t, i , 1 )C 2 (s, m 2 H+ ,m 2 w ,ml , m|, mfj, 

/ 2 y2(s) (i7) = / 3 ^W(fl) = jf (s) (JO = A y2(s) (#); 

A %W (^) = 1fl T^P 1 2 A &iiff 6 " [(M*m*0ii - h b m 2 (3 21 )Co 
16V2ir 2 (s - m z A ) 2 

+/i fe m 2 y /3 2 iCi + h b (3 21 (ml - m 2 )C 2 - 2h b (3 2 iC 00 - h b (3 21 (i + u 
-2ml)C 12 - 2h b m 2 H+ P 2 iC 22 ](s,m 2 H+ ,m 2 v ,m 2 ,ml,m 2 )} 
-ghhPnfa A ht(3uB stt + [{htm 2 pil _ hb m b m t (3 21 )C 



16V2ir 2 (s-m 2 A ) L 2 
+/i t m 2 v /3i 1 Ci - h t (3u{m b - mf)C 2 - 2h t fi\\C m - h t (3n{i + u 
-2m 2 )C 12 - 2h t m 2 H+ PiiC 22 ](s,m 2 H+ ,m 2 v ,m 2 ,m 2 ,ml)} 

i,j,ki',j',k' \ A) 

+ 1! E Wir^-m 2 )^ 21 (^') 2 ^i^V^i 6l ^' i'*'^ /lfe0 ^'i 

i,j,ki' ,j',k' ^ M> 

+h t @i> i > 1 )C 2 (s, m 2 H+ ,mw > ^> m |' m S fc )> 

f Vi(t) _ v- -ghbhfaua^Pu Wbt 2 2 

~~ 16v27r z (£ — m^J 
+(m 2 , + t)Ci2 + iC 22 ](ml, m 2 w ,t, m 2 Ht ,ml,m 2 )} 

\ - -ghh 2 (3ii(3 2i Pii nWbt.r 2 r , or , 
+ \l^ 2 (i-m 2 ) { ~ B " +[-^Co + 2C o 

+m 2 C n + (m 2 + £)Ci2 + iC22](mJ,m^,f,m^,mj,rai)} 
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+ E 1 * gh f^ a ^i lW \ [(htm b mtPij - h b m 2 t (3 2j )C + h b f3 2j (m 2 b - mf)Ci 

+h b (3 2j (i - mf)C 2 + 2h b (3 2j C 00 + (-h t m b m t (3ij + h b mlf3 2j + h b i(3 2j )C 12 

+(h b t(3 2j - h t m b m t (3ij)C 22 ]{m b , m^J, m 2 ,m 2 H+ ,m 2 H J 

j 

+ E ^f n/?1 * [htm.mMCo + 2d + 2C 2 + C n + 2C 12 + C 22 ) 
i l6\/2TT z (t — mf) 

+2h b f3 2i C 00 + h b m 2 p 2i (C + C\ + C 2 ) + /i 6 mg/? 2i (Ci + Cn + Ci 2 ) 
+hP 2i i(C 2 + C i2 + C 22 )](mb, mvj/,i, m 2 ,m 2 + ,m 2 A J 

i 

+ E f^ a 2 f 11 ^ [^^(^o + C 2 - C n - C 12 ) 

-2/ii/3 H C 00 + /if>m 6 mi/3 2i (-C + Cn + 2C i2 + C 22 ) - h t (3 u i{C 2 
+C 12 + C 22 )\(m 2 b ,m 2 w ,i, ml,m 2 Hj ,m 2 H+ ) 

+ E ?l b 1f n ^ 1 2 j 1 {^3^4^i^i[^g(C 1 + C n + C 12 ) + t(C 2 + C 12 
i,j,fc ov^vr [t-m t ) 

+C 22 ) + 2C 00 ] - (iV^iV^^^m^o + 7V fe3 ^4^ 1 ^ 2 m t m ;e o)(Co 

+Ci + C 2 ) + N&Numbmt^uid + C 2 + C n + 2C i2 

+C 22 )}K,^,t,m|o,m|,m|) 

+ E y^%{-MLO,t^ 2 ^ + ^ + MHi^Cm^msoC^Co 

+0,t(2C o + m 2 b C u + m 2 b C 12 + tC 12 + tC 22 - m|C )) + h t m b Of*V j2 9{ 2 (C 
+Ci + C 2 )](m^,m^,t, m|,m^+,m|o)} 



+ E o iTT 2 ; [^m-^LlCo + Cx + C 2 ) + ^m b 4 2 ^ 2 (C 

+2Ci + Cn + 2Ci 2 + 2C 2 + C 22 )](m 2 b , m 2 w , t, m|, m|+, m|o) 

+m 2 b (C u + C 12 ) - mf fc C + t(C 12 + C 22 )) + m~+m^O^C Q ) 

+0 b k A 2 O^(h b m t m^N j3 U i2 + ^m 6 m-+ AT; 3 y i2 )(C + d + C 2 ) 

+/ l6 m b m t (^ 2 ) 2 ^; 3 ^2O^(C + 2d + C n + 2C i2 
+2C 2 + C 22 )](mf,m^,i, m~ fc ,m|o,m|+)}, 

v l(t) = E ^W^fe (4Co + 4Cl + 4C2 + Cn + 2Cu 
^ 16V2ir 2 (t - mf) 

+C 22 )(ml,m 2 v , i, m\. , m\ , m 2 ) 
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\^ gh b htm b mtPiiP2iP2i (r or r w 2 2 12 m 2 2, 

+ Z_^ 72 AT ill 2~7~^ u + 2Cl2 + c 22){m b ,m w ,t,m Ai ,m b ,m t ) 

^ 16v 2ir z (t — mf) 

+ E if^^J^ 21 ^ [( h bm b m t P2j - htmfpi^Co + h t (3ij{m 2 b - m 2 )d 
^ 16V2vr 2 (t - mf) 

+/i t /?ij(t- mf)C 2 + 2h t f3ijC Q Q + {-h b m b m t (3 2j + h t m\p Xj + M7?ij)Ci 2 

+ (h t iPij - h b m b m t P2j)C22](rn b ,m^J,m^,m 2 H+ ,rn 2 H .) 

j 

+ E ~ 9 £ bk 9 f. lii3 \ [h b m b m t (3 2i {C + 2d + 2C 2 + C n + 2C 12 + C 22 ) 
i \§\J2-K l {t - mf) 

+2h t (3 li C O0 + h t m 2 t (3 u {Co + d + C 2 ) + hm^d + C n + C 12 ) 
+h t (3iii(C 2 + C 12 + C 22 )](mf,m^,t, mf,m 2 H+ ,m 2 A J 

i 

+ + C 2 - C n - C 12 ) 

-2/i b /3 2i Coo + h t m b m t (3ii{-Co + Cn + 2C i2 + C 22 ) - h b f3 2i i{C 2 
+C 12 + C 22 )](m^,m^,t, ml,rn 2 H] ,m 2 H+ ) 

+ E fl^^l 1 . {^iV^^^mfem^d + C 2 + C n + 2C 12 + C 22 ) 

-(iV^V^u™*™^ + ^3iV fc4 ^ 1 ^ 2 m 6 m ;e o)(Co + d + C 2 ) 
+7V* 3 7V fe4 ^ 2 4[m^(Ci + Cn + C 12 ) + i(C 2 + d 2 + C 22 ) 
+2C 00 ] } {m b , , t , m| , m~ , m~. ) 

+ 2^ a 2 /r 2 -—\- h tV k 2 (J ji v j2 b o + [nhmhm^oV^Oji U j2 (C 

i,j,k 87r l t_m C 

+Ci + C 2 ) + h t 0i 2 Of*V j2 {2C 00 + mgCn + mgCi 2 + tC 12 + tC 22 

- m 2 k C + m^o m -+ C )] (m^ , m^ , i , m| , m|+ , m|o ) } 

„ gh b h t m t f3i 2 eLOf*N* 
+ E 2T 2 M^iM^t) + 2d + C n + 2C 12 + 2C 2 

+C 22 ) + h t m^+e{ 2 Vj 2 {C + Ci + C 2 )](m^,m^,t,m|,m^+,m|o) 

- E 8 J^ 2) {^(^L) 2 ^3^ 2 o^i? ^^ - [/i b (^ 2 ) 2 ^;3^2(o^( 2C, oo 

i ,j,k ^ t I 

+m 2 b {C u + C12) - mf fe C + i{C 12 + C 22 )) + m~+m^oO£*C ) 

+0 b ki e b k2 Ofj*{h t m t m^>N* 3 V t2 + h b m b m~+N j3 U i2 ){C + Ci + C 2 ) 

+^m 6 m t (^i) 2 ^3^20^(Co + 2d + C n + 2C 12 
+2C 2 + C 22 )] {m 2 b ,m 2 w , i, m^ , m|o , m|+ )}, 

fVljt) _ ST 9 h l h t m t a ^ a 2ih\ m Wbt , \( Arn 2 1 _2 , im 2 r 

/ 5 "E 32 ^(t-m?) {B ° +^ + ™H l )C» + Am b C l 



3 Supersymmetric Electroweak Corrections to W ± H T Associated Production at the 
CERN Large Hadron Collider 48 



+2(m 2 h + i)C 2 -2C 00 ](m 2 b ,m 2 v ,i,m 2 Hi ,m 2 b ,m 2 t )} 

+ E ~t}X 2 t-<^ [r +(mi ^° - 2Cm){ml m ^ 1 m% ' mlml)] 

+ E + W^ooK, t, m?, < + , <) 

^ 16V2n 2 (t - mf) "j 

+ 2^ ^ at o/r ^{htiribPii - h b m t P2i)Coo{m b ,m w ,t,m t ,m H+ ,m A ) 

^ l6\ / 2Ti z (t — mf) 1 

+ E -. f h ^ 23 ^ 2lLP3 1 {h t m t /3ii + hmh^Cooim^mwJ, m 2 b ,m 2 H m 2 +) 
~~ lQ\J2n z (t — m z ) J « 

i,j,k v ^ \ t ) 

1 2 2 1 2 2 2 \ 

(m b , m w , t, m-o , , m t ~ ) 

i,j,k imr l t_m *i J 
+£C + tC 2 ) + h h m^e\ x Of*XJj< 1 (m 2 b Cz + m 2 Ci + tC 2 ) 
-^ t m 6 4 2 Oj;*^ 2 (-2Coo - m^Ci + t(C + Ci + 2C 2 ) + m|C 
- "^0 m xj C oM m h m w > *> m | > m ~+ , m|o ) } 

^ gh b htmtPi20ki N iA r h at n R* T r R w *l^ , r h _. gt n L* TI r 

+h b e t kl Of*U j2 {-2C m + m 2 C + 2mgCi + mfC 2 + tC 2 + m\C Q ) 
+h t m b e t k2 V j2 {m ji +Of* - rn^oOjf )(C + C\ 

2 m? - 2 - 2 

,2 , 



+C 2 )](m 2 ,m 2 y ,f,m|,m 2 + ,m 2 o)} 



+ E 16 / 2 (|! 21 m2) {(-^(^i) 2 ^3^or + h(el 2 ) 2 N* 3 u t2 o%*)B^ 

-\h t m t (9 h kl ) 2 N j3 V l2 (Of*{-2C m + m^(Co + 2d + C 2 ) + m? fe C + tC 2 ) 



-m-+m^O^C ) + ^i^(^^3[/ i2 (m.+O^(m^ + mid + tC 2 ) 
-m^pf*{m 2 b d + iC + tC 2 )) + h t m b m t N* 3 V i2 (m^Of-* (m b m t C 
+m b m t d + m 6 m t C 2 ) - m-+O%*(C + d + C 2 ))) 
+h b m b (d b k2 ) 2 N* 3 U l2 (m^m ji0 OffCo + O%*(2C 00 - m 2 d - m\C Q 
-i(C + Ci + C 2 )))](m6,m^,f,m? ,m| ,m?+)}, 

«k A.j Xj 



m 



+(2mg + m 2 + i)Ci + (m 2 + 3i)C 2 - 2C 00 ](m 2 , m 2 ^, t, mi , mg, m 2 )} 
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Y> -gh b h^m b pii(3 2 iPii , pW bt , r _ 2 r , 2 rw r , r x 



-2C 00 ] (m 6 , m w , t, m Ai ,m b ,m t )} 

+ E -, f 1 ^ 1 ^!^^ i. h t m tPij + h b m b p 2 j)C 00 (ml,mw,i, m\ , m 2 + , m# ) 
16V27r^(t - m£) j 

+ StTTJ-^ on ( h bm b (3 2 i - htrntPii)Coo{ml,rnw,i, m 2 ,m] l+,m 2 A ) 

i 16^2ir z (t — mf) i 

+ E f^t^off 11 ^ {htmbPu + h b m t f32i)C 00 (ml,m 2 v ,t, m 2 ,m 2 H ,m] 2 +) 
16V27r 2 (t - mf) 



<;3 

+ E 'tltta^ll ( m tN* k3 N k4 e b j2 8t 2 - mtNnNh^ODCoo 
~ 8V27T 2 (t-mf) 

12212 2 2 \ 
(m b , m w , t, m-o , m g , m f . ) 



+^m 6 4 1 0,f ^ 2 (-2C o + m^Ci + tC + tCi + 2tC 2 + m\C Q ) 
-htm^OjaPftV^mlCo + m 2 b d + tC 2 ) + htm-^O^V^mid + tC 
+iC 2 )} {ml , m 2 w , i, m| , m|+ , m|o ) } 

+ 2^ 167r 2 (t-mf) j * J + [ h brribVki u j2K m x° u ji 

-m~+Of*)(C + Ci + C 2 ) - ^^ 2 0^y j2 (2C o - m 2 C - 2m 2 Ci - m^C 2 

- E I ^r^y {[-^(^) 2 ^3^og* + ^(^J^^oflBo^ 

-[/ l6 m t (^ 2 ) 2 ^; 3 ^ 2 (0^(-2Coo + m^(Co + 2d + C 2 ) + m? fc C + fC 2 ) 
-m- t m^Oij*C ) + e b ki e h k2 (h t N* z V a (m^O^{m 2 C Q + m 2 d + tC 2 ) 
-m^oOfj \m\C\ + tC + tC 2 )) + h b m b m t Nj 3 Ui 2 (m^oOfj* (m b m t C 
+m b m t C l + m 6 m t C 2 ) - m~+Of*{C Q + d + C 2 ))) 



+h t m b (9 b kl ) 2 N 33 V l2 (m jlt m^O%*C + O£*(2C 00 - m 2 d - m? fc C 
-i(C + Ci + C 2 )))](m6,m^,f,m| fc ,m?p,m|+)} 



f Vi(t) _ \- -9h 2 htrn t aua2ifa\ (or , ^ r w 2 2 ,~ 2 2 ^ 
J7 - 2^ 32^/2tt 2 (£- m 2 ) 2 L22 >^ m b^ m W^^ m H^ m b^ m t) 

x ^ -ghlh t m t (3 li p2iP2i (r r w 2 2 f 2 __2 _.2<i 

+ 1^ /o o/x 2V + L '22){m b ,m w ,t,m Ai ,m b ,m t ) 

i l6\/2TT z (t — mf) 
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+ gh b h tai p 2m + + + hbmt(3 {Cq + Ci 

^ 16V2vr 2 (t -m|) 

-C12 - C22)}(ml,m^,i,m^,m 2 H+ ,m 2 H .) 

j 

+ E 9h " h f' ihl 2 ^ [^m^HCCi + Cn + C 12 ) + h b m t (3 2i (C + d + 2C 2 
~^ 16v27r 2 (i — m|) 

+C12 + C 2 2)}(ml,m 2 v ,i, m 2 ,m 2 + ,m 2 A .) 

i 

+ + + C 22 ) + h b m b p 2i (C + C 2 - C u 

-C 12 )](m b , m 2 w , i, ml , m 2 H , m 2 + ) 

J i 

+ E Z^a 1 ^ i m bN* k3 N k4 9 b j2 9UC2 + c n + c l2 ) 
-NkzN^eynv (C + Ci + C 2 ) + m t iV fc3 Ar* 4 ^i0'i(C 2 + C 12 

~ 2^ a 2 /r ^[hbm^+O^Oji U j2 Ci + h b m^d kx O^ t/ i2 C 2 

+h t m b 6 t k2 Of*V j2 (C 1 + Cn + C12)] (mg, m^, i, m| , m|+ , m?o) 
- E ^f^^^t f/ j2 ^(C 12 + C 2 + C 22 )K, £, m | , m 2 m| ) 



- E a Al 2 J W^) 2 ^ 2 o£*(c 12 + c 2 + c 22 ) 

+h b e b kl 9 b k2 N j3 U i2 (m-+Ofj*C 2 + mqOffd) 

+h b m b (9 b k2 ) 2 N* 3 U i2 0**(C n + C 1 + C 12 )](m 2 , t, m?, , m|o, m|+), 

f Vi(t) _ \- -g^^"ib«H«2i/3ii , lfm 2 m 2 r 2 TO 2 2^ 

is - 2^ 32^/27r 2 (£ - m 2 ) Ci2){rn b ,rn w ,t,m H .,rn b ,m t ) 

x ^ gh b h 2 m b (3ijf32i(3n (r r u 2 2 ? 2 2 2 , 
i 16v 2vr 2 (i — m^J 

+ E^fet%[« 3 (Ci + C n + C 12 ) + htmtlhjiCo + d 
^ 16V27r 2 (t - m£) 

-C12 - C 2 2)](ml,m 2 v J,m 2 ,m 2 + ,m 2 H .) 

j 

+ [hniMC! + C u + C 12 ) + h t m t (3 u {C + d + 2C 2 

i 16v27r z (t — m£) 

+C i2 + C 22 )]{ml,m 2 v ,i, m 2 ,m 2 H+ ,m A J 

i 

+ j2Zf^0mi [hbmt(32ii c 2 + C n + C22) + h t m b /3 u (C + C 2 - C n 
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-C12 )] {ml , m 2 w , i, ml , m 2 H ,m 2 H+ ) 

J i 

+ E tf^f^t [m t N* k3 N k4 e b j2 eUC2 + C 12 + C 22 ) 

i,j,k °v Z7r V' ~ m t ) 

-N^N^elm^iCo + Ci + C 2 ) + mtN^N^elid + C n + C 12 )} 

( 2 2 1 2 2 2 \ 

{m b ,m w ,t,m- ,m~ .,m t ~.) 

+h b m b 6 t kl Of*U j2 {C 1 + Cn + Ci 2 )](mg, m^,£, m| fc ,m|+,m|o) 

Egh^mtl3i\{9 k2 ) 2 L 2 2 x 2 2 2 \ 

——j 2^ j* ^2^ 4 (C 2 + C12 + C22)(m b ,m w ,t,m i m +,m~o) 

+ E ^^ [^(^2) 2 ^3^20r(gl2 + C 2 + C 22 ) 

+h t e b kl e b k2 N^v l2 (m^o^C2 + m^p\;c x ) 

3 f> \2 A7 T/ r\L*(r< \ r< \ r* \U™2 „,2 £ _2 _2 _2 



+^m 6 (^ 1 ) 2 7Y ? - 3 y i 20^(Cii + Ci + Ci 2 )](mg,m^,t,m? fc ,m|o,m|+), 

/s w = e ^fCT-% {Br + [rr4Co + 2m ' Ci + (m? + £)C2 

-2C 00 ](^, m^,£, m 2 H .,ml,m 2 )} 
^ -gh b h 2 m b pii(3 2 iPii r pW bt , r _2 r / 2 ^ 

+ E ^MM) {i?0 +[-^o-K-t)c 2 

- 2C 00 ] (m 2 , mvi/ , i , . , m 2 , m\ ) } 

^ gh b h t aiif3n(32j<Pij n ( 2 2 X 2 2 2 \ 

—ghbh 2 f3liP2iPll „ , 2 2T 2 2 2 \ 

+ E 1C at 271 h~:Coo(m b ,m w ,t,m t ,m +,m A .) 

i lG\/27r z (t — mf) « 

^ gh b h 2 a 2 j(3ii(3uipji ,22x22 2 \ 
E 16^7r 2 (t-m f 2 ) °° ' ™ W ' ' m& ' H ' ' "V } 

^-^ gh b h 2 /3ii(9 t ji) 2 (Oj^) 2 1 2 2 ~ 2 2 2\ 

+ E o at 2^ ^NmN^Coo (m b ,m w ,t, m~ ,m~ m- ) 

+MfciOji*^2(-2Coo + m 2 Ci + tC 2 + m|Co) + h t m b 9 t k2 V j2 {m-+Ofi 
- m^o Of* ) Ci] (m 2 , m 2 ^ , t , m| , m 2 + , m|o ) } 



^ S^ ^^^^2(Co + C 2 ) - h b m^ kl Of*U 32 C 2 

i,j,k v * ' 
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+h t m b ei 2 Of* V j2 {C + d + C 2 )] (m 2 b , m 2 w , t, m\ , m|+ , m|o ) 



+O£*(-2C 00 + mid + m? fc C + tC 2 )) ~ h kl 9 h k2 (h b m t N^U l2 (m^O^C 2 

-Og*(m^oC + ^+C2)) + h t m b N* 3 V i2 {m^oOfj*d - O%*(m~+C + m^oCi))) 

+/i 6 m 6 m t (^ 2 ) 2 Af; 3 C/ i2 0^*(Co + Ci + C 2 )](mg,m^,*,m? fc ,m|o,m|+)}, 

j?Vi(t) —ghlh t aiia 2 il3 2 i . , , 2 2 r 2 2 2\ 

/12 = E 32y feir 2 (i _ m 2) (2C + Ci + C 2 )(ro 6 , ro w , t, m 6 , m t ) 

^ gh b h^aiiP 2 iPij<fij ,22122 2 \ 

+ E 1fi at uz 2V Coo(m b ,m w ,t,m t , m + , m H ) 

^ 16V 2ir 2 (t - mf) i 



+ E ,a AT 271 2T C oo(m 6 , m w , t, m t ,m +,m Ai ) 

x ^ gh\a 2j fi 2X fi 2 npji 2 2 r 2 2 2 \ 

+ E 1C at 271 2 -zCoo{m b ,m w ,t,m b ,m H ,m +) 

^ -ghlhtfae^e^el^ 2 2 . 2 2 2 

+ E 8^2^ _ m 2) — ^fcs^Coo K.%.<. ^ ' m s, > m * i ) 

+ E 167r 2 (t-m t 2 ) { fcl Ji j L-"& m # m *+ e fci°ji ^' 260 
+MlliOjf^2(-2Coo + r^Ci + tC 2 + m|C ) + htm^V^m^+Of* 
- m^o Of* )Ci](ml, m 2 w , i, m| , m 2 . + , m?o ) } 

i,j,k v * ' 

+/i 6 m 6 ^ 1 O|*C/ j2 (C + Ci + C 2 )](mf,m 2 y ,t, m|,m 2 + ,m 2 o) 



+O**(-2C 00 + mfd + m? fe C + tC 2 )) - 9 b kl 9 b k2 (h t m t N* 3 V i2 (m~+0*j*C 2 
-Off(m^oC + m~+C 2 )) + h b m b N j3 U i2 (m^Of*d - Off(m~+C + m^d))) 
+h t m b m t (e b kl ) 2 N j3 Vi2O^(C + d + C 2 )](mg, , t, , m|o, m|+)}; 
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t"f 8\/27r 2 (i - rr?|) 2 



,2//it \2 



/ 6 s(t) = e "^y^ 2 [-(^ 2 - ^ + 

Z ^32V27r 2 (t-- 2 ^ LV * ; 



m 



+ 2 , 2 M^fe(m 2 + t)**^ + (« + ^KtB?^] 

+ E 1fi ^tf 1 iJ ^MifWd 2 ^ - h b h t etA2{m 2 t u j2 v 32 
+tu; 2 v* 2 )B^~ bi + h^m^iv^B^], 

,s(t) _ 1 ,s(t). 
J12 — 2 

/J* W = E I^fr^^ 1 ^" + Mn«0) + 2mgd + m 2 C 2 



m 



+tC 2 ) - h b m b m t f32i(4Co + 2C\ + 2C2)](m b ,m H+ ,t,m H .,m b ,m t )} 

+ E E 8 "^tf 'l^ (^ e 'vi + ^9^, 2 )(m t iV fe *3iV fe4 ^ 1 ^ 1 C 2 
,j,ki',j' v ^ v 



m 



m^o N^N^e^Co + mbNksN^e^d ) (m 2 , m 2 ^ , f, m|o , m 2 . , mf. ) , 



f V 2 (t) _ _m b y 2 ( t ) 
/5 ~~ 2 ^ 2 ' 

,y 2 (t) \ - -gh b h t a li a2i , DJ H-+ta . r , « / 2 ^ 

/e =2^ 32^ 2 (£ _ m 2 ) \ h b m thiB + [h b m t (32i{m H .C 

+2tC 2 + 2m 2 b C\) - h t m b p n {(m 2 + i)(2C + Ci) 
+2tC 2 )] (ml,m 2 H+ , i, m 2 H .,m 2 ,m 2 )} 
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-gh b h t p li p 2 i f , „ jjH+bt . r. Ian 
^ 32V^it^m^ 21 [hbm t m A .p2iCo 

+h t m b f3n (m? - t)d] (mg , , i, m% , m 2 , m 2 )} 

+ E E 1fi V ^'i + ^©'v2)(m 6 m t iV fe *3iV fc4 ^ 1 ^ 1 C 1 

ij.fei'j' ib V^7T [t m t ) 

-mtm^N^Nk^e^Co + iN k3 N^e) 2 0\ 2 C 2 )(m 2 b ,rn 2 H+ X m|o , m?., m| ), 



f v 2 (t) _ _! ,v 2 (t). 
J12 — 2 1 ' 2 ' 



gh 2 cx 2 - 

ff ] = E 16V2tt-2 { h bP 2lC °( m H+' m w, s,rn b ,m 2 ,m 2 ) + \-h t m b m t P\\ {D 13 

+D 23 ) - h b p 21 2D 00 + h b m 2 b f3 2l {2D 3 - D n - D 12 + D 13 + D 23 ) 
+h b m 2 H+ (3 21 {D 13 + D 23 ) - h b tf3 21 (D 12 + D 13 + D 22 + ^23) 
+h b m 2 H .p 21 D ]{m b ,rn 2 H + , m^, m 2 , i, s, m 2 Hi ,ml,rn 2 ,m 2 )} 
gh 2 (3 2 - 

+ E 16% /2^2 {- h bP2iC (m 2 H+ , m 2 w , s, ml,m 2 ,ml) + [-h t m b m t (3ii{Di 3 

+D 23 ) + 2h b p 21 D 00 + hmlp^Du + £> 12 + £> 13 + D 23 ) 
-h b m 2 H+ (3 21 {D 13 + L> 23 ) + h b tp 21 {D 12 + L> 13 + D 22 + L> 23 ) 
-h b m\(3 21 D }(ml,m 2 H+ , m 2 w , m b , i, s, m 2 A .,ml, m 2 ,m b )} 

+ E gh 8 b ^ 2 i [h 2 rn b p 11 p li (D 1 + D u + D 12 + D 13 ) - h b h t m t f5 lx f3 2i {p x 
+D 12 + D 13 ) + h 2 b m b p 21 p 2i (D 12 + D 13 )} 

(22 2 2 - - 2 2 2 2\ 

(m b , m w , m H+ ,m b ,u,s, m + ,m t ,m b ,m t ) 

i 

i,j,k,l V ,k' 

+D 23 ) - m^N l3 e^ 2 e b n (D + D t + D 2 ) + m b 0^ 2 d b j2 N^(D 1 + D 2 + D u 
+2D 12 + D 22 )](ml,m 2 H+ ,m 2 v ,mlJ,s,m 2 l o,ml^m 2 k ,m 2 b ) 

+ E E 9 ° ll ^2 ij ' ( h » e h>i + ^ef^vOt/igm^g^^CDia + D 23 ) 

-hhtm.+e^Uk^Ds + /^m^^Ps + A*)] 
(m 2 , m 2 ^ , , m 2 , i, m 2 + , m| , m| , mf. ) 

+ E ^16^^ {- fe 6fti2ff2jC2(mg, "4+ , *, rn 2 H .,ml,m 2 ) 

i,3 * 7F 

+[/t t 2 mg/3ii/3i j (D 23 + 2D3 + 2D33) - h b h t m b mtpiif3 2j (D 23 + 2D 3 ) 
+h b h t m b m t (3i 2 (3ijD 33 - hl(3i 2 (3 2j (ml(D 23 + D 33 ) + m^D^ + 6Z> 23 

+m^+ L> 3 )] (m 2 ^ , m 2 , ,m 2 ,u,t, m 2 H+ , mf? . m 2 , m 2 ) } 

j 3 
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+ E iG^n 2 ^ h ^ 21 ^ C ° + Cl + C2 )( m J?+' m b'*' m h m b, m % 
+ [h 2 b (3 21 p 2i (m 2 b (D 12 - D u ) + m 2 w D 13 - uD 12 - m 2 +D 1 ) 

i 

+h b h t m b m t (l3 2 i(3iiDn - /3 n /3 2i D 12 ) + h 2 ml(3iif3iiD 12 } 
(ml,m 2 H+ , ml, m 2 w ,u, i, m 2 H+ , mf , m 2 , m 2 A .)} 

i 



+ E E g fe 8 g' ^ (^e^ 1 + fe t e^ 1 )M 1 ^ 1 Lr J - 2 (m 6 JV5Qg > i> 2 3 

% -jk ^1 k^ . 

/ 2 22 2 ~ r 2 2 2 2 \ 

(m w , m b , m H+ ,m b ,u,t, m _+ , m ~o , , m t ~ ) , 
,(6) _ 



= E ftr^ i-htmmtPupD! + 2D 2 + D n + 2D 12 + D 22 ) 

+h b m 2 b p 21 (4D + 6Z?i + 2£>n + 4D 12 + D 13 + 6D 2 + 2D 22 
+D 23 + 2D,)] (mlm 2 H+ , m 2 w , m 2 ,t, s, m 2 Hi , mg, m 2 , m 2 ) 
-nh 2 B 2 

+ T / -^0i h t m bm t Pu(D 11 + 2D 12 + D22) 

+h b m b p 2 i (D 13 + D 23 )] (ml , m 2 H+ , m^, ml,i, s, m\. , m 2 b ,m 2 u m 2 b ) 

+ E g^^2 ^ h ^ nf3li ^ C ° + Cl + C ^)( m W^ m2 H+^,m 2 t ,m 2 b ,m 2 t ) 

+[-h b htm b m t (Pufc\D 1 + /3n/3 2 ^ii) + ^m^ 2 i/%(.Di + £>n) 
+h 2 {3 11 p li (-mlD 11 + m^Dxz - uD 12 - uD 13 

2 r-> M/ 2 2 2 2 ~ ~ 2 22 2\i 

-m^+L>i)J(m & , m w ,m H+ ,m b ,u,s, m H+ ,m t ,m b ,m t )} 



2 



iij->k,l i f ,k' 

+D U + 2Z?i2 + £22) - m^N? 3 e b a 8 b j2 (D +D ± + D 2 ) + m b 9 b 2 9 b 2 N l3 (D 13 
+D 23 + D 3 )] (ml , m 2 H+ , m 2 w ,ml,i,s, m|o , m~ . , m| , m| ) 

+ E E 96f f^f"' (h b ef, fl + ^efyO^m^^-^AB 

-hhtm^ele^U^Da + h 2 m b e\ 2 e) 2 V k \(D X3 + Das)] 

,-22 2 2 ~ ~ 2 2 2 2 \ 

(m 6 , m H+ ,m b ,u, s, m + , m f . , mr , m f . ) 



+ E ^ig^l^ ^ Ai AjC2(mg, , t, m 2 H .,ml,m 2 ) + [-hlml(3 12 /3 2j (D 23 

+2£> 3 + 2D 33 ) - h b h t m b m t (3ii(3 2 jD 33 + h b h t m b m t (3i 2 (3ij(D 23 + 2£> 3 ) 

+hlf3 11 p lj {ml(D 23 + D 33 ) + m^Dis + fiD 23 

+m^+ D 3 )] (m^ , , ,m 2 ,u,i, m 2 H+ , m# . ml,m 2 )} 
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~E iQ^-2 ^ h ^ ll(3li ^ Co + Cl + C 2 ){m 2 H+ ,m 2 b ,i, m 2 ,rn b ,m 2 A .) 
+ [h 2 t PnPii(m 2 (D 12 - D u ) + m 2 w D 13 - uD 12 - m 2 + A) 

i 

+h b h t m b m t ((3 ll p 2i D u - f3 21 f3 u D 12 ) + hl/3 21 /3 2i D 12 ] 
(ml,m 2 H+ , ml,mw,u, i, m 2 H+ , m 2 t , m 2 , m A .)} 

i 

+ E E 9hb ^f w (h b el, n + h t ei /n )[h t e b k2 ej 2 v j2 (m b N i3 o^D 23 

+m jc+ N* 3 0?;D 3 ) + h b rn b b kl O t ll NZ i U j2 Og'D a3 ] 



(m w ,m b ,m H+ ,m b ,u,t,m-+,m^pm~ bk ,m {i ), 

= E 16 ^^ 2 {-^ m &/ 3 2iC 2 (m^ + ,m^,g,mb,mt,mfc) + [h t m b m t (3n(2D 3 

+D 33 ) - 2h b m 2 b (i 2l D 33 + h b m 2 w (i 2l D l3 - h b t/3 21 (D 13 + D 23 ) 
-h b m 2 H .(3 21 D 3 ](rn 2 ,m 2 H+ , m^,mlX s, m 2 Hi ,ml,m 2 t ,m 2 )} 

+ E { h bP2iC 2 (m 2 H+ , m 2 w , s, m 2 b ,m 2 ,m 2 b ) 

+ [h t m b m t (3 u D 33 - h b rn 2 w fi 21 D 13 + h b t(3 21 (D 13 + D 23 ) 
+h b m\(3 21 D 3 }(ml,m 2 H+ , m 2 w , m b ,t, s, m A .,m b , m 2 ,ml)} 

+ Y,-0l^l h t m b m tPiMDo + D l + D 2 + D 3 ) - h b h 2 m 2 b (3 u Pu(Do + A 

+D 12 + D 13 + 2D 2 + D 22 + 2D 23 + 2D 3 + D 33 ) - h b h 2 m 2 pnPii(D + D 2 
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I 2 2 2 2 ~ ~ 2 2 2 2\ 
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(mw, ml,m 2 H +,ml,u, i, m 2 H + , m^.m 2 , m 2 )} 

i 

+ E iq^2 { h bfoi&2i c i (. m H+ > m b » i m ti m l m %) + [hlP2ip2i(.ml(D 12 
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+ E E - 167] . 2 — N l3^ii'^jl^kl^kk' 

(h b Q% kll + hte^ k/1 )[m b 9 b t 29 b j2 N l3 D 33 

i,j,k,l i',k' 
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+h b h t m-+e t ii e t j2 U k2 V k2 (D + D X + D 2 ) - h 2 m b B^e) 2 V k \(D x + D 2 + D n 
+2D 12 + D22)]{ml,m 2 v ,m 2 H+ ,ml,u,s,m 2 ^+,m^,ml [ ,m 2 j ) 

+ E g ^g^|^2 3 {-tffaiPijCi (ml, m 2 H + , i, m 2 Ht ,ml,mf) + [h 2 b m 2 b l3 12 (3 2j (2D 2 

+D22 + 2D 23 ) + h b htm b mt(Pii02jD2a + 0i2/3ij(D 22 + 2D 2 )) 
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~ E { h t(3iif3iiCi(m 2 H+ ,ml,i, m 2 , ml , m A .) + [h 2 t f3uf3ii{m 2 b (Di 2 

-D22) + m 2 v D 23 + uD 22 + m 2 H+ D 2 ) - h h h t m b m t ((3 2 i(3iiDi 2 - (3n(3 2i D 22 ) 
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+/ l ^ t mfm t /3 li /3 21 /3 2i ( J Di + i} 2 + £> 3 ) - hlm b p 21 pl(-2D 00 + mg(Z>i + £> 2 
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8 ^2 ■ "iw,.,!^, 

(m| , , ,m 2 b ,u,s, m 2 + , m\ ,m\,m\) 
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+ E 9 ( h t m bPnPij - hhtPufoj ~ 2hlm b f3 12 (3 2j )D 00 



1,3 

J2 



I 2 2 2 2 " 2 2 2 2 2\ 

(m w ,m b ,m H+ ,m b ,u,t,m H+ ,m Hi m b ,m t ) 
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(ml,m 2 H+ , m 2 w ,m 2 b ,i, s, m%,ml,m 2 ,ml)} 

+ E ^^ {^PiMhbmMd + C 2 ) - /itm^uCo] 

(m^,m 2 H+ ,s,m 2 ,ml,m 2 ) + [h b h 2 m b m 2 (3 2 i(3 2 iD - h 2 h t m b m t /3ii(3 2 iP 2 i(2D 
+D 1 + D 2 + D 3 ) + h 3 b m 3 b p 21 pl(D + D 1 + D 2 + L> 3 ) - /itWii/^K^i 
-m^L> 3 + £(#2 + #3) + m 2 H+ D ) + h b h 2 m b pnPiiP2i(4Doo + m 2 (£>i 2 

+D 13 + D 23 + D33) - m 2 H+ D 13 + m t 2 r»i - m 2 w {D 23 + D 3 - D 33 ) 
+u(D 12 + D 13 + D 2 + 2D 22 + 3D 23 + D 3 + L> 33 ) + m 2 H+ {D 2 + D 3 ))\ 

1 2 2 2 2 - - 2 22 2\i 
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+ E E ^16^-2 6 \ N l3\ 2 ^ii'^i2 d jl^j2^kl^kk'( h b&i'k'l + frf®i'fc'l)A)0 
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_ ^9h 2 t e b e b Jel 2 ) 2 e\ l e t .. l v 2 / , . 
+ E E sX 2 - (^ef yi + ^ef yi )^oo 
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+ E E ^(^(Mfc + ^e^{-ML4^ 2 o£* 
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+h t m b d b kl 9 t l2 N* 3 V j2 (-m-+0^ + m^O%)D 2 + h b m b b k2 9 t ll N i3 U j2 (m-+O^ 

/-,fi*-, n 1/2 2 2 2 ~ X 2 2 2 2 \i 
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(m b ,m H+ ,m w ,m b ,t, s, m A .,m b ,m t , m b ) 
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+ X! \hlm b (32if32iDi2 - h b h t m t {Pii(32iDi + (kijhiDw) 

+h 2 m b f3 1 if3 li (D 1 + D u + D12 + D^^m^m^+^^mly^J^ 2 +,m 2 ,ml,m 2 A J 
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+ E E §^fci^fc'^*i^V(M3^ 
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/ 2 22 2-X2 2 2 2\ 

(m w ,m b ,m H +,m b ,u,t,m H +,m H .m b ,m t ) 
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~ E i h t m bPiiPuDi2 - h b h t m t ((3 2 if3uD 1 + PufciDu) 
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1,2 2 
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+ E |^ih^Wll/?Hpl2 + £>13 + D 2 + D 2 2 + D 23 ) 

+h b h t mtf3 11 (3 2i {D 2 + D 22 + D 23 ) - h 2 m b PuP 2 i(D 22 

,7-1 M/ 2 2 2 2 ~ ~ 2 2 2 2\ 

+.D23 ) J (m 6 , m w , ,m b ,u,s, m H + ,m t ,m b ,m t ) 

+ E E ^16^-2 6 I 1 2 1 (^1 ) 2 <9fcl ^fcfc^ ( ^6 Qf" fc' 1 + h t @i>k'l) D 23 
I 2 2 2 2 ; - 2 2 2 2 \ 
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- E E^^^iWM 2 *^^ 

+ E ^§^[-^11^(^2 + ^2 + ^22 + ^23) 

+h b h t m t PiiP 2j (D 12 + £>2 + ^22) + hlm b p 12 (3 2j {2D 12 + D 22 )] 
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+ E j^fja [-h 2 b m b lhilhiD22 + h b htmtPiiP2i{D 2 + D 22 + ^23) 
-h 2 m b f3iiPii(Di2 + D 2 + D22 + D23)](m 2 ,m 2 H+ ,ml,m 2 v ,u,i, m 2 + ,m 2 ,m 2 ,m 2 A .) 
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~ E E |^|^fel^feifc'^U^V(^&®l'/'l + h t Q% V1 )N* z Uj 2 0^* (D 12 + £> 2 

+D22 ) (m^ , ™& , ™h+ , m b,u,i, m 2 + , m 2 -o , m~ fc , m| ) , 

-(b) \- -g^^«2^21 , 2 2 2 2 i - 2 2 2 2n 

/10 = iGv^tt 2 Dl ^ m b^ m H+^ m w^ m b^^,m HiJ m b ,m t ,m b ) 

Egh^TYifofSoifi^ . 2 2 2 2 ^ ^ 2 2 2 2\ 
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+ E ^ 8^2" ^ htmt(3M + h b m bP2i(D 23 + D3 + D 33 )} 
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+ E 9 i£% V 2 i-hfabPnPijiDu + D22) ~ hhtmtPuPijiDu + D 2 + D 22 ) 
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ECjhu /?21 /~^2? 7—. / 9 9 2 2 ^ ^ 2 2 2 2\ 

16a/2tt 2 ' ' m& ' U ' * ' ' m * ' m& ' ^ 
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+ E ^^^ {-hhtPuPiifoiiCo - Ci)(m^ ,m 2 H+ J,m 2 ,m 2 ,m 2 ) 
+ \h\m b m t $ lx 0{ i D 3 - h 2 b h t m b m t (3 2 i{(3iif3 2l D 2 - fhifeDi) 
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+h 3 b m 2 b (3 21 (3lD 2 - h b h 2 t (3 n f3iM-4D 00 + mg(£>i - D 12 - D 23 

+D 3 ) + m 2 H+ D 12 + m 2 w D 23 - u(D 12 + D 22 + D 23 ) + m 2 + (D - D 2 ))\ 

i 

2 ™2 2 2 ~ s 2 2 2 2\ 



(m 6 , m w , m H+ ,m b ,u,s, m H+ ,m t ,m b ,m t )} 

i 

9h%Ol 2 ift\\fi\ji*pij . 2 22 2 - r 2 2 2 2\ 

+ 1^ ,r /o 9 - D o{m w ,m b ,m H+ ,m b ,u,t,m +,m H m b ,m t ) 

gh b h 2 f3u f3\i(3 2 i . 2 2 2 2 - 1 2 2 2 2 \ 

+ 16 ^ 7r2 ^00 (m b , m H+ ,m b ,m w ,u, t, m H + , m t , m bJ m A% ) 

+ E E ^^fi *f + ^TiK^™^^^ 

i,j,k,lk',l' 

+h t 9 b ki e t l2 N* 3 V j2 [-m^O^(D + D 1 + D 2 ) + m^Of*{D 1 + £> 2 )] 
+h b m b el 1 9j 1 N* 3 U j2 0^*D 3 }(m 2 v ,ml,m 2 H+ ,ml,u, i, m|+,m|o,m| fc ,m|); 

All other form factors /j not listed above vanish. 

Here Ao, Cj, and Z)jj are the one-, three- and four-point Feynman integrals [20]. The 
definitions of Uij, Vij, N^, Ofj and Ofj can be found in Ref.[2]. 
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Figure 1: Feynman diagrams contributing to supersymmetric electroweak corrections to bb — > 
W~H + : (a) and (b) are tree level diagrams; (c) — (x) are one-loop corrections. The dashed 
line 1 represents H, h, A; the dashed line 2 represents H, h, A, G°; the dashed line 3 represents 
H + , G + . For diagram (r), the dashed line in the loop represents H, h, A, G°,H + , G + , t, b. 
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Figure 2: Feynman diagrams contributing to renormalization constants: The dashed line 
represents H, h, A, G° , H + , G + , t, b for diagram (a), and Hi in diagrams (d) — (/) represents 
H,h,A. 
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Figure 3: The Yukawa corrections versus mjj+ for tan/3 = 1.5, 2, 6 and 30, respectively. 
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Figure 4: The Yukawa corrections versus tan/? for m H + = 100, 150, 200 and 300GeV, 
respectively. 
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Figure 5: The genuine SUSY EW corrections versus m H + for tan/3 = 1.5, 2, 6 and 30, 
respectively, assuming M 2 = 300GeV, fi = -lOOGeV, A t = A b = 200GeV and M^ = M {] = 
M„ = 500GeV. 




Figure 6: The genuine SUSY EW corrections versus m H + for tan/3 = 1.5, 6 and 30, respec- 
tively, assuming M 2 = 200GeV, n = lOOGeV, A t = A b = ITeV, Mq = M„, m~ ti = lOOGeV 
and mr = 150GeV. 
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Figure 7: The genuine SUSY EW corrections versus At = A b for tan/3 = 1.5 and 30, 
respectively, assuming m H + = 200GeV, M 2 = 300GeV, fj, = lOOGeV, and Mq = A% = 
M„ = 400GeV. 




Figure 8: The genuine SUSY EW corrections versus Mq = Mjy = for tan/3 = 1.5 and 
30, respectively, assuming m H + = 200GeV, M 2 = 300GeV, fi = lOOGeV, and A t = A b = 
500GeV. 
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Figure 9: The genuine SUSY EW corrections versus [i for tan/3 = 1.5 and 30, respectively, 
assuming m H + = 200GeV, M 2 = 300GeV, A t = A b = 500GeV and Mq = M = M b = 
400GeV. 



4 Preface 

In this part, I will present some works on flavor changing neutral current (FCNC) 
processes. As already emphasized, the FCNC processes are forbidden at tree level in 
the SM, so it acts as the ground to test the quantum structure of the SM; at the same 
time, it is also the ideal place in searching new physics beyond the SM. 

In the first section of this part, tc associated production in the SM at linear colliders 
is presented, which corrects some mistakes in literature; the second section is the tc 
production in supersymmetrical models, the FCNC process is mediated by gluino; the 
third section is bs production in the SM at linear colliders; the fourth section is about 
inclusive process of B meson decay in a CP-softly broken two-Higgs-doublet model; 
the last section is given to B exclusive decay to K and lepton pair in supersymmetrical 
models. 



Part II 



FCNC 
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5 Top-Charm Associated Production at High En- 
ergy e + e~ Colliders in Standard Model 

ABSTRACT 

The flavor changing neutral current tcV(V=7,Z) couplings in the production vertex 
for the process e + e~ — ► tc or ic in the standard model are investigated. The precise 
calculations keeping all quark masses non-zero are carried out. The total production 
cross section is found to be 1.84 x 1CT 9 fb at ^/s=200 Gev and 0.572 x 10" 9 fb at 
y / s=500 Gev respectively. The result is much smaller than that given in ref. || by 
a factor of 10~ 5 . 

Top quark physics has been extensively investigated |L). The advantage of exam- 
ining top quark physics than other quark physics is that one can directly determine 
the properties of top quark itself and does not need to worry about non-perturbative 
QCD effects which are difficult to attack because there exist no top-flavored hadron 
states at all. The properties of top quark could reveal information on flavor physics, 
electroweak symmetry breaking as well new physics beyond the standard model(SM). 

One of important fields in top physics is to study flavor changing neutral current 
(FCNC) coupings. There are no flavor changing neutral currents at tree-level in the 
SM. FCNC appear at loop-levels and consequently offer a good place to test quantum 
effects of the fundamental quantum field theory on which SM based. Furthermore, they 
are very small at one loop-level due to the unitary of Cabbibo-Kobayashi-Maskawa 
(CKM) matrix. In models beyond SM new particles beyond the particles in SM may 
appear in the loop and have significant contributions to flavor changing transitions. 
Therefore, FCNC interactions give an ideal place to search for new physics. Any 
positive observation of FCNC couplings deviated from that in SM would unambiguously 
signal the presence of new physics. Searching for FCNC is clearly one of important 
goals of high energy colliders, in particular, e + e~ colliders 0. 

The flavor changing transitions involving external up-type quarks which are due to 
FCNC couplings are much more suppressed than those involving external down-type 
quarks in SM. The effects for external up-type quarks are derived by virtual exchanges 
of down-type quarks in a loop for which GIM mechanism H is much more effective 
because the mass splittings between down-type quarks are much less than those between 
up-type quarks. Therefore, the tc transition which is studied in the latter opens a good 
window to search for new physics. 

The FCNC vertices tcV(V=7, Z) can be probed either in rare decays of t quark 
or via top-charm associated production. A lot of works have been done in the former 
case [[|. And a number of papers on the latter case have also appeared || [£], [7|. In 
this letter we shall investigate the latter case in the process 

e + e~ — > tc or ic. (1) 

Comparing t quark rare decays where the momentum transfer q 2 is limited, i. e., it 
should be less or equal to mass square of t quark m 2 , the production process (1) allows 
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the large (time-like) momentum transfer, which is actually determined by the energies 
available at e + e~ colliders. The reaction (1) has some advantages because of the ability 
to probe higher dimension operators at large momenta and striking kinematic signa- 
tures which are straightforward to detect in the clean environment of e + e~ collisions. 
In particular, in some extensions of SM which induce FCNC there are large underly- 
ing mass scales and large momentum transfer so that these models are more naturally 
probed via tc associated production than t quark rare decays. 

The production cross sections of the process (1) in SM have been calculated in 
refs. H [7] . In the early references |?J a top quark mass m t < m z is assumed and the on- 
shell Z boson dominance is adopted. The reference || considered a large top quark mass 
and abandoned the on-shell Z boson dominance. However, the "self energy" diagrams 
have been omitted in ref. ||. This is not legal because the one- loop contribution for 
FC transitions is of the leading term of the FC transitions and must be finite, i.e., 
although there are some divergences for some diagrams they should cancel each other 
in the sum of contributions of all diagrams. Furthermore, the order of values of cross 
sections given in ref. |J is not correct. 

The order of values of cross sections for the process (1) in SM can easily be esti- 
mated. The differential cross section can be written as 

dcosd 32tts s A s ^ ns 

Where N c is the color factor, 6 is the the angle between incoming electron e~ and 
outgoing top quark t and M is the amplitude of the process. In eq. (2) the charm quark 
mass in kenetic factors has been omitted. Due to the GIM mechanism, one has 



E l M l 2 = e 8 l E v; t v ]C f(x 3 , y ,)\ 

spins j=d,s,b 

ml — m? s df 

m\j dxj 



_8|TT* V " b b - '"-8 "J I , |2 ,o\ 

6 \ v tbVcb — — 2 — aT7k.VJ=0 + -I > W 



where Xj = m 2 -jm 2 wl yj = m 2 -/ 's, and denote the less important terms for yfs > 200 
Gev. Assuming J^- |ay,iy=o = 0(1), one obtains from eqs. (2), (3) 

a ~ lO" 8 - 10 9 /6 
at y/s = 200 Gev. However, the results given in ref. || are 

a = 0.71 x 10" 2 /fr 

for m t =165 Gev and 

a = 4.1 x IQ^fb 

for m t =190 Gev, which are much larger than the above estimation by a factor of 10 5 . 
In order to test SM and search for new physics from observations of some process one 
needs to know what are the precise results for the relevant observables of the process 
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in SM. Therefore, it is necessary to calculate precisely the cross sections in the SM. 
In this letter we calculate the differential and total cross sections of the process (1) in 
SM. 

In SM for the process (1) there are three kinds of Feynman diagram at one loop, 
"self enengy" (actually it is a FC transition, not a usual self energy diagram), triangle 
and box diagram, which are shown in Fig.l. We carry out calculations in the Feynman- 
t'Hooft gauge. The contributions of the neutral Higgs H and Goldstone bosons G 0,± 
which couple to electrons are neglected since they are proportional to the electron mass 
and we have put the mass of electron to zero. 

We do the reduction using FeynCalc [[| and keep all masses non-zero except for the 
mass of electron. To control the ultraviolet divergence, the dimensional regularization 
is used. As a consistent check, we found that all divergences are canceled in the sum. 
The calculations are carried out in the frame of the centre of mass system (CMS) and 
Mandelstam variables have been employed: 

s = (Pi + P2) 2 = (h + k 2 ) 2 t=( Pl -h) 2 u = ( Pl -k 2 )\ (4) 

where pi,p 2 are the momenta of electron and positron respectively, and fci, k 2 are the 
momenta of top quark t and anti-charm quark c respectively. 
The amplitude of process e + e~ — > tc can be expressed as 

M = ^ 2 a 2 V* j V t j[giUa^PLVcV e J^PRU e + g 2 U t J^ P L V c V e ^^P L U e + g 3 U t PLV c V e h P R u e + 

j=d,s,b 

giUtPhVcVe h P LU e +g5Ut A P LV c V e ftlP L U e + g§V e ^ P L U e Ua ^ i>\PhV c + 
g^Utl^PRVcVel^PRUe + g^Utl^ ' PRV c V e ^^P L U e + g 9 U t PRV c V e faPRU e + 

9wU t P R V c V e h P LU e + gilVe^PLUeUafi 0iPrV c ] (5) 

where a is fine structure constant, V^- is CKM matrix element, P L is defined as 
(1 — 7 5 )/2, and Pr is defined as (1 + 7 5 )/2. The exact expressions of the coefficients 
gj(j = 1,2, ...11) are too long to be given. Instead, in order to show the essential 
points, we give them in the limit of mj/m (i=d,s,c, m=m w , m t , s) approach to zero. In 
the limit gj(j = 7, 8, 9, 10, 11) is zero, and the others are given as follows. 

gi = a 3 m 2 — 2a<±m 2 s\ ] + ^a^C^rri-m^s^ + Qm 2 u (2Cfim 2 + 2C 22 s + 2Cf 2 (m 2 + s))(a 3 + 2a^c 2 1) s 2 l) ) + 
12m 2 w C^(a 3 (m 2 + s) + a±(2sc 2 w s 2 w + m 2 c 2 w s 2 w - mfa* )) - B%(m1j - m 2 m 2 t + m 2 m 2 w + 
2m 2 m 2 w - 2m^,)(ai + 2a 2 s^,(3 - 4s^)) + 12CQ (a 3 (m 2 + Qm 2 w ) + a^i^m 2 + I2c 2 w m 2 w - 
m 2 s 2 J) + 6Cfm 2 w (2a 3 (m 2 + s) + 2a A s 2 w (c 2 w m 2 + 2c 2 w s - m 2 s 2 J) - 6C$m 2 w (2a 3 (m 2 - s) - 
2a 4: s 2 v (2c 2 u 'm 2 + 2c 2 w s + 2m 2 s 2 w - m 2 s 2 w - 3>m 2 c 2 w )) + 2CQm 2 (a 3 (m 2 + 2m 2 w - m 2 ) + 
a 4 s 2 w (3m 2 - 2m 2 s 2 w - 4m 2 w s 2 w + 2m 2 s 2 w )) - 2(2C C + C^m 2 + C c 22 s + C c 2 s + 
C i2( m t + s)){a 3 (m 2 + 2m 2 w ) + 2a A s 2 UJ (3m 2 J - m 2 s 2 w - 2m 2 w s 2 w )) + 
B^{a\(m 2 — m 2 v )(m 2 + 2m 2 w ) — 2aim 2 m 2 + 2a 2 s 2 l] {2>m A j — 6m 2 m 2 + 3m 2 m 2 w — ^m\, - 
Aml-sl - \m 2 m 2 w s 2 w + 8m*s^, + Sm 2 m 2 s 2 w )) - 2Cl{a 3 {2m 2 w s + m 2 m 2 ) - 

a4S 2 u (3m 2 m 2 + 2m 2 m 2 s 2 u — 6m^s + 4m^ss^,)) (6) 
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g 2 = a z m) + o 4 mJ(l - 24)(4 - 3C 2 c m 2 ) + + <D{ + Df + 2L>f + 2L>? 2 + 4£>f 3 + 

2L>f 3 + 2Z>f 3 ) + (2m?L>f + 2sDf 3 + 2L>f 3 m?)) + 6m 2 (2C 1 d 1 m f 2 + 2C 2 d 2 s + 
24(m t 2 + s))(a3 - a 4 c 2 (1 - 24)) + 6m 2 C^(2a 3 (m 2 + s) - a 4 (l - 2s 2 )(mfc 2 „ - m?4 + 
2sc 2 )) - -Bo( m i ~~ m 2 m 2 + m 2 m 2 w + 2m 2 m 2 w - 2m^)(ai - a 2 (3 — 4s 2 )(1 - 2s 2 )) + 
6C$ (2a 3 {m 2 + 6m 2 ) - a 4 (l - 2s 2 )(c 2 m 2 + 12c 2 m 2 - m 2 s 2 w )) + 6Cf m 2 (2a 3 (m 2 + s) - 
a 4 (l - 2s 2 w )(c 2 w m 2 t + 2c 2 s - m 2 s 2 )) - 6C^m 2 (2a 3 (m 2 - s) + a 4 (l - 2s 2 w )(-c 2 w m 2 t + 2c 2 s + 
2m 2 4 - m 2 s 2 w )) + C c m 2 (2a 3 (m 2 + 2m 2 w - m 2 ) - a 4 (l - 2s 2 )(3m 2 - 2m 2 s 2 w - 4m 2 w s 2 w + 
2m 2 4)) - 2(2C C + C{ x m\ + C 2 c 2 s + C c 2 s + C c 12 {m 2 + S ))(a 3 (m 2 + 2m 2 ) - 
a 4 (l - 24)(3m 2 - m 2 s 2 w - 2m 2 w s 2 J) + B%{ ai (m 2 - m 2 J(m 2 + 2m 2 ) - 2aim 2 m 2 - 
a 2 (l — 2s 2 u )(3mj — 6m 2 m 2 + 3m 2 m 2 w — 6m^ — 4777^4 — 4m 2 m 2 u s 2 + 87744 + 8m 2 m 2 s 2 u )) — 
C±(2a 3 (m 2 m 2 + 2sm 2 J + a 4 (l — 2s 2 u )(3m 2 m 2 — 6st?4 + Asm 2 u s 2 u + 2m 2 m 2 s 2 u )) (7) 

g 3 = l2a 4 s 2 w m t {2C^m 2 w - C^m 2 ) + Am t C^m 2 {a 3 - 2a 4 4) + 24m t m 2 (2Cf + C^)(a 3 + 2a 4 44) + 

8Cim t (a 3 m 2 - 2a 4 s^,m 2 ) + YlmtiC^ + Cf 2 )(a 3 (m 2 + 2m 2 ) + a 4 s 2 (4m 2 + 4c 2 w m 2 w - m 2 s 2 w )) + 
477 H (C 1 C 1 + C 1 c 2 )(a 3 (m 2 + 2777^) + 2a 4 4(3m 2 - m 2 4 - 2m 2 w s 2 w )) (8) 

g 4 = -2a 5 m t (2D e 23 + D e 2 + 2D e 33 + 2D e 3 ) + 6a 4 m t (C 2 c m 2 - 2C 2 d m 2 )(l - 24) + 4C c m t m 2 (a 3 + 
o 4 4(l - 24)) + 2477^774,(2^ + C d )(a 3 - a 4 c 2 (1 - 24)) + SCfm^m 2 + 
a 4 4m 2 (l - 24)) + 6(4 + C 1 d 2 )m t (2a 3 (777 2 + 2777^) - a 4 (l - 2s 2 J(c 2 w m 2 + 4c 2 w m 2 w - 
m 2 s 2 J) + 4m t (C 1 c 1 + C 1 c 2 )(a 3 (m 2 + 2777^) - a 4 (l - 24)(3m 2 - m 2 s 2 w - 2m 2 w s 2 w )) (9) 



55 = -4a 5 (Df 2 + Dt 3 ) (10) 

56 = a 5 m t (2^ 2 + 2£»f 3 + 2 J D| 3 + Z)| + 2 J D| 3 + 2D|) (11) 

with m 2 = ml (since m s , rrid have been omitted in the above expressions of g's), 
where a«(i = 1, 2, 5) are defined by 

1 1 1 

1 96s7T 2 s 2 777 2 777 2 ' 2 1 7 687r 2 c^s\ ] m\m 2 w {m 2 z — im Z I y z — s) ' 3 192s7r 2 s 2 m 2 

1 1 

4 3847T 2 44 m2 y( m2 — im^rz — s) ' 5 327T 2 S^ 



with c w = cos6 w and s w = sin6 w . In the presentation of g,j above, we have used the 
definition of scalar integrals Bs, Cs,and -Ds|§, and these functions, Bs, Cs, and Ds, 
with superscripts a,b,...,e have the arguments 

(0,m 2 ,m 2 ), (m 2 ,m 2 ,ml), (m 2 , 0, s, m 2 , m 2 , m 2 ), (m 2 , 0, s, m 2 , m 2 , m 2 ) 
(0, s, m 2 , it, 0, 0, 0, ml, m 2 , m 2 ) 

respectively. Here mj denotes the mass of down- type quark b. 

In the numerical calculations the following values of the parameters have been used 

& 

m e = 0, m c = lAGev, m t = 175Gev, mj = 0.005Get> , m s = 0.17Gev, 
m b = 4AGev, m w = 80AlGev, m z = 91.187Gef , T z = 2.5Gev, a = 
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Figure 1: Feynman diagrams of prosess e + e — ► tc 

In order to keep the unitary condition of CKM matrix exactly, we employ the 
standard parametrization and take the values |9|, 10 1 

s n = 0.220, s 23 = 0.039, s 13 = 0.0031, 5 13 = 70°. 

Numerical results are shown in Figs. 2, 3. In Fig.2, we show the total cross section 
a to t of the process e + e~ — ► tc as a function of the centre of mass energy <Js. One 
can see from the figure that the total cross section is the order of 10~ 10 ~ 10~ 9 fb, 
as expected, and decreases when center-of-mass energy increases and is large enough 
(> 250 Gev ). We fixed the centre of mass energy yfs at 200Get> . Differential cross 
section of the process at the energy as a function of cos 8 is shown in Fig.3. 

To summarize, we have calculated the production cross sections of the process 
e + e~ — > tc in SM. We found that the total cross section is 1.84 x 10~ 9 /6 at y/s = 
200 Gev and 0.572 x 10~ 9 fb at y/s = 500 Gev. It is too small to be of experimental 
relevance. Therefore, this is a remarkable situation that allows for a precise test of 
the SM and, in particular, of the GIM mechanism in SM. Even a small number of tc 
events, detected at LEP II or a NLC running with a yearly integrated luminosity of 
£ > 10 2 [/6] _1 , will unambiguously indicate new FCNC dynamics beyond SM. 
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Figure 3: Differential cross section of the process e + e — > tc, where = 200 GeV. 
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6 SUSY-QCD Effect on Top-Charm Associated Pro- 
duction at Linear Collider 

ABSTRACT 

We evaluate the contribution of SUSY-QCD to top-charm associated production at 
next generation linear colliders. Our results show that the production cross section 
of the process e + e~ — > tc or ic could be as large as 0.1 fb, which is larger than the 
prediction of the SM by a factor of 10 s . 

One of the most important physics in top quark sector is to probe anomalous flavor 
changing neutral current (FCNC) couplings. In the Standard Model (SM), FCNC 
couplings are forbidden at the tree level and much suppressed in loops by the GIM 
mechanism. Any signals on FCNC couplings in the processes of top quark decay and 
productions or indirectly in loops will indicate the existence of new physics beyond 
the SM. Recently in the framework of effective lagrangian, Han and Hewettfll] have 
examined carefully the possibility of exploring the FCNC couplings tcZ/tcy in the 
production vertex for the reaction e + e~ — > tc + ic and concluded that at higher energy 
colliders with 0.5 — 1 TeV center-of-mass energy, the resulting sensitivity to FCNC 
couplings will be better than the present constraints |J. In this paper, in the minimal 
supersymmetric standard model (MSSM) we study the process e + e~ — > tc + ic and 
perform an detail calculation of the contribution from the FCNC couplings in the 
vertex of gluino-squark-quark to the production cross section. We will point out that 
at higher energy e + e~ colliders the cross section could be as large as 0.1 fb which is at 
least eight order of magnitude larger than the prediction of the SM ~ 10~ 10 — 10~ 9 fb 

I- 

The MSSM is arguably the most promising candidate for physics beyond the SM. 
Beside many attractive features of supersymmetry in understanding the mass hierar- 
chy, gauge coupling unification, the weak scale SUSY models in generally lead to a 
rich flavor physics. In fact, SUSY models often have arbitrary flavor mixings and mass 
parameters in the squark and slepton sectors and these mass matrices after diagonal- 
ization induce FCNC couplings at tree level in the vertex of gluino-squark-quark etc. 
Phenomenologically one would have to assume certain symmetries or dynamical mech- 
anisms to prevent large FCNC among the first and second generations. On the other 
hand the flavor structure, especially among the second and third generations in the 
SUSY sector motivates us to seek for new physics and any experimental observation 
on the FCNC processes beyond the SM would undoubtedly shed light on our under- 
standing for flavor physics. In this paper we take model of Ref. 0, [|] where the FCNC 
couplings relevant to our calculation is given by: 

C rc = -V2g s T a K^P L qq L + h.c. (1) 

In (1), K is the supersymmetric version of the Kobayashi-Maskawa matrix, which is 
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explicitly expressed as: 



( I e e 2 
-e 1 e | (2) 

\-E 2 -E 1 



where e parameterizes the strength of flavor mixing and is shown to be as large as 1/2 
without contradicting with the low energy experimental data 0. 

In Fig.(l) we give the Feynman diagrams for the process e + (pi)e~(p 2 ) — > t(ki)c(k 2 ). 
In calculations, we have neglected the scalar u-quark contribution since it is highly sup- 
pressed by Ki 2 Ki 3 ; and we use the dimensional regularization to control the ultraviolet 
divergence. We have checked that all divergences cancel out in the final result with 
the summing up of all of the diagrams. The calculations are carried out in the frame 
of the center of mass system (CMS) and Mandelstam variables have been employed: 

s = (Pi + Vif = (h + k 2 f t={ Pl -k l ) 2 u = { Pl -k 2 ) 2 . (3) 

After a straightforward calculations, one obtains for the amplitudes 

M = ^v{p 1 ) 1 ^u{p 2 )u{k l )V^{tc 1 )v{k 2 ) 

HPihM - g e A ls)u(p 2 )u(h)V>*{tcZ)v(k 2 ) (4) 



2 cos 9 W {S- Mf) 



where, g e v = 1/2 - 2sm 2 6 w , g\ = 1/2, and V^tcrf) and V{tcZ) are the on-shell 
quarks effective vertices given by 1 

v»(tcr,z) = fT z i,P R + P2 Z i,Pl + ff z k 1 ,p R + fJ' z k^p L 

+fI' Z k 2 ,P R + W Z k 2f ,P L . (5) 

The form factors, are 

r-v (iljee^ cos(^-) sin(^)ms . n , n 2 2\ n / 2 2 2 m 

fi = E 12m 7T 2 [ ° ( ' ^ ^ ~ ° (m " ^ ^ + 

q=c,t 1 

H = E 2Am 2 n 2 K m 9- m 9~ 2 ) 5 o(0,m 5 ,m g - 2 )-(m 5 -m g - 2 + m t )B {m t ,m § , 



9=5,* 



+4m 2 C o] + R-R 
f , = (Tl)ee^ sin(0,) cos^m, ^ + ^ + R R 

q=c,t 

= (±l)eeg 2 sin(^) cos(6> g -)m g + + ^ 

q=c,t 

1 For simplicity, we only give the results in the limit of m c = 0. However in our numerical calcula- 
tions, we use the full formulas. 
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fl = E (±1)fe ^f Wm ' [C + 2 Clll+ R.R. 

q=c,t 
q=c,i 

ff = E ^Sff Mt,) W. < O ~ *° « ™% "4)1 + ** 
^ S. ( 9eSff {( - 3 + - "4 Wo. "4 "4) 

-(mf - m? 2 + m t 2 )S (m t 2 ,m|,m|)] + 4m 2 (-3 sin 2 (^~) + 4 sin 2 (^))C 00 



-f 2m 2 cos 2 Coo} + A-R 

,2 , 



/s - L 48cos(^)vr2 [(4sm (^) - 3sm (^))(C + 2^ 



q=c,t 

,2 



+3sin 2 (^)(C + 2C 1 )] + J R. J R. 
h ~ E 48cos(^) 7 r2 [(4sm(^)-3sm(^))(C + 2C 2y 

q=c,t v ' 



+3sin 2 (^)((% + 2C , 2 )] + J R. J R. 



.2 - 2/ 



/? = E (=Fl lg^ W2 )mt [(^in 2 (^) - 3sin 2 (%))(C + 2C U ) 

q=c,t ^ 

-3cos 2 (^)(C , + 2C'i 1 )] + J R. J R. 
/« - L 48cos( ^ )7r 2 t( 4sin - 3sin (WO, + 2C 12 ) 

q=c,t x ' 

-3cos 2 (^-)(C'o + 2(712)] + R.R. (7) 

where represents the replacement of 9q — > 7r/2 + 9g and <->• m^. The vari- 
ables of three point functions (7$, (7^- and (7j, (7^ are (m 2 , S 1 , 0, m~, m| 2 , m~ 2 ) and 
(m 2 , S 1 , 0, ra|, m| 2 , rn^J, respectively. 

In the MSSM the mass eigenstates of the squarks qi and g 2 are related to the weak 
eigenstates qi and qn by [|7| 

( f 1 U # f g ~ L 1 with # = f C ° S ^ Sin ^ "I . (8) 

For the squarks, the mixing angle 9g and the masses 2 can be calculated by diago- 
nalizing the following mass matrices 

2 / M 2 L m g M LR \ 
* U,% M 2 fl i' 
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M 2 LL = m| + m\ + m 2 z cos 2/3(J g 3L - e q sin 2 9 W ), 
M RR = m }j y r) + m l + m l cos 2f3e q sin 2 9 W , 

MLR - MRL -\A b -^ a nP (g = 6), (9) 

where mj~, Jn 2 ^ are soft SUSY breaking mass terms of the left- and right-handed 
squark, respectively; fi is the coefficient of the H1H2 term in the superpotential; A t 
and Ab are the coefficient of the dimension-three tri-linear soft SUSY-breaking terms; 
Ig L , e q are the weak isospin and electric charge of the squark q. From Eqs. |3| and we 
have 



mj 

ti.a 



ml -M£ L 



M 2 LL + M 2 RR =F J(M 2 LL - M 2 RR f + 4m t 2 M 2 H 



tan% = fl - — . (10) 

Now we present the numerical results. For the SM parameters, we take 

m z = 91.187GeV, m w = 80.33Ge\/ , m t = 176.0GeV, m c = lAGeV 

a = 1/128 , a5 = 0.118 (11) 

For the MSSM parameters, we choose fi = —lOOGeV and e 2 = 1/4. To simplify the 
calculation we have taken that m^ = mg, = mg = A t = m$ (global SUSY). In Figs. 
2-5, we show the cross sections of the process e + e~ tc as functions of m s , rrig, y/s 
and tan (3. One can see that the production cross section increases as squarks and 
gluino masses decrease, and it could reach 0.1 fb for favorable parameters. This is an 
enhancement by a factor of 10 8 relative to the SM prediction. Such enhancement could 
be easily understood as following: 



vsusy ( a s Am : 



a SM V am b 



:i2) 



where Am? represents the possible mass square difference among squarks. If Am? 
varies from 100 2 - 200 2 (Ge1/) 2 , = 10 7 ~ 10 8 . At the same time, this kind of 

aSM loll 

enhancement could also be observed in FCNC decay process of top quark Q. Due to 
the rather clean experimental environment and well-constrained kinematics, the signal 
of tc or tc would be spectacular |1[. We expect the SUSY-QCD effects studied in this 
paper be observed at higher energy e + e~ colliders. 
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Figure 1: The Feynmann diagrams for the process e + e — > tc. 




m s (GeV) 



Figure 2: The cross section for the process e + e — > tc as a function of ms, where 
v^s = 500GeV, tan /3 = 2, e 2 = 1/4 and /i = -lOOGeV. The solid and dashed lines 
represent rrig = lOOGeV and 500GeV, respectively. 
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Figure 3: The cross section for the process e + e — > tc as a function of rrig, where 
= 500GeV, tan/5 = 2, e 2 = 1/4 and /i = -lOOGeV. The solid and dashed lines 
represent m s = 300GeV and lOOGeV^, respectively. 
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Figure 4: The cross section for the process e + e~ — > tc as a function of v^5, where 
m § = lOOGeV, tan/3 = 2, e 2 = 1/4 and /i = -lOOGeV. The solid and dashed lines 
represent ms = 300GeV and lOOGeV^, respectively. 
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Figure 5: The cross section for the process e + e — > tc as a function of tan/3, where 
v^S = 500GeV, m~ g = lOOGeV, e 2 = 1/4 and n = -lOOGeV. The solid and dashed 
lines represent m s = 300GeV and lOOGeV, respectively. 
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7 Bottom- Strange Associated Production at High 
Energy e + e~ Colliders in Standard Model 

Abstract 

We investigate the flavor changing neutral current bsV(V=7,Z) couplings in the 
production vertex for the process e + e~ — ► bs or bs in the standard model. The 
precise calculations keeping all quark masses non-zero are carried out. Production 
cross sections are found to be the order of 10~ 3 fb at LEP II and the order of 10 _1 
fb when center-of-mass energy is near the mass of neutral gauge boson Z. 

7.1 Introduction 

There are no flavor changing neutral currents (FCNC) at tree-level in the stan- 
dard model (SM). FCNC appear at loop-levels and consequently offer a good place 
to test quantum effects of the fundamental quantum field theory on which SM based. 
Furthermore, they are very small at one loop-level due to the unitarity of Cabbibo- 
Kobayashi-Maskawa (CKM) matrix. In models beyond SM new particles beyond the 
particles in SM may appear in the loop and have significant contributions to flavor 
changing transitions. Therefore, FCNC interactions give an ideal place to search for 
new physics. Any positive observation of FCNC couplings deviated from that in SM 
would unambiguously signal the presence of new physics. Searching for FCNC is clearly 
one of important goals of the next generation of high energy colliders JTJ. 

The flavor changing transitions involving external up-type quarks which are due to 
FCNC couplings are much more suppressed than those involving external down-type 
quarks in SM. The effects for external up-type quarks are derived by virtual exchanges 
of down-type quarks in a loop for which GIM mechanism || is much more effective 
because the mass splitting between down-type quarks are much less than those between 
up-type quarks. Therefore, for example, the bs transition which is studied in the paper 
has larger probability to be observed than that for the tc transition. 

The b-hadron system promises to give a fertile ground to test the SM and probe 
new physics. The FCNC vertices bsV(V=7, Z) have been extensively examined in rare 
decays of b-hadron system || ^, [5] . The observation of FCNC processes in both the 
exclusive B — > K*^ and inclusive B — > X s ~f channels has placed the rare B decays on a 
new footing and has put a stringent constraint on classes of models M. Analyses of the 
inclusive decay B — > X s l + l~ show that in the minimal supergravity model(SUGRA) 
there are regions in the parameter space where the branching ratio of b — > sl + l~(l = 
e, /i) is enhanced by about 50% compared to the SM |7| and the first distinct signals 
of SUSY could come from the observation of B — > X s fi + fi~ if tan/3 is large ( > 30 ) 
and the mass of the lightest neutral Higgs boson m^ is not too large (say, less than 150 
Gev) The B factories presently under construction will collect some 10 7 — 10 8 B 
mesons per year which can be used to obtain good precision on low branching fraction 
modes. 
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The FCNC vertices bsV(V=7, Z) can also be investigated via bottom-strange as- 
sociated production. In the paper we shall investigate the process 

e + e~ — ► bs or bs. (1) 

Comparing b quark rare decays where the momentum transfer q 2 is limited, i. e., 
it should be less or equal to mass square of b quark m 2 , the production process (1) 
allows the large (time-like) momentum transfer, which is actually determined by the 
energies available at e + e~ colliders. The reaction (1) has some advantages because of 
the ability to probe higher dimension operators at large momenta and striking kine- 
matic signatures which are straightforward to detect in the clean environment of e + e~ 
collisions. In particular, in some extensions of SM which induce FCNC there are large 
underlying mass scales and large momentum transfer so that these models are more 
naturally probed via bs associated production than b quark rare decays. 

It has been shown that the cross sections of e + e~ — > tc in SM are too small to be 
observed at LEP or NLC |§. As pointed above, in SM the cross sections of e + e~ — ► bs 
should be much larger than those of tc final states. Are they large enough to be seen at 
LEP or NLC? In the paper we would like to address the problem by calculating cross 
sections and backward-forward asymmetry of the process (1) in SM. 



7.2 Analytic calculations 

In SM for the process (1) there are three kinds of Feynman diagram at one loop, 
self energy-type, triangle and box diagram, which are shown in Fig.l. We carry out 
calculations in the Feynman-t'Hooft gauge. The contributions of the neutral Higgs 
H and Goldstone bosons G 0,± which couple to electrons are neglected since they are 
proportional to the electron mass and we have put the mass of electron to zero. 

We do the reduction using FeynCalc |TJj and keep all masses non-zero except for the 



mass of electron. To control the ultraviolet divergence, the dimensional regularization 
is used. As a consistent check, we found that all divergences are canceled in the sum of 
contributions of all Feynman diagrams. The calculations are carried out in the frame 
of the center of mass system (CMS) and Mandelstam variables have been employed: 

s = (Pi + P2) 2 = (h + k 2 f t = (p 1 - h) 2 u = (j>i k 2 ) 2 , (2) 

where pi,p2 are the momentum of electron and positron respectively, and ki, hi are the 
momentum of bottom quark b, and anti-strange quark s respectively. 
The amplitude of process e + e~ — ► bs can be expressed as 

M = 1 6n 2 a 2 Vf f y js \giUb^PRV s v e 'y fl P R u e + g2u b ^PLV s v e j^P R u e + g 3 u b ^P R v s v e j^P L u e + 

j=u,c,t 

giUbl^PLVsVel^PLUe + g 5 U b P R V s V e %P R U e + g 6 U b P L V s V e ftlPRUe + g 7 U b P R V s V e h P LUe + 

g 8 u b P L v s v e fi\P L u e + g 9 u b i>\P L v s v e $\P L u e + 9ioV e ^ 'P L u e u b ^ ^ i>\P R v s + 
gnv e ^P L u e u b 'y fl tfiP L v s ] (3) 
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where a is fine structure constant, Vij is CKM matrix element, Pl is defined as 
(1 — 7 5 )/2, and Pr is defined as (1 + 7 5 )/2. The expressions of the coefficients gj(j = 
1, 2, ...If ) can be found in Appendix. 

Having the amplitude M, it is straightforward to obtain the differential cross section 

by 

dcosO 167r s f4 s ^J 1 y> 

Where N c is the color factor and 9 is the angle between incoming electron e~ and 
outgoing bottom quark b. 

7.3 Numerical results 

In the numerical calculations the following values of the parameters have been used 



ff 



m e = 0, m u = 0.005Get> , m c = lAGev, m t = 175Gev, m s = 0.17Gev, 

m b = 4AGev, m w = 80AlGev, m z = 91 A87Gev, T z = 2.5Gev, a = — 

f 28 

In order to keep the unitary condition of CKM matrix exactly, we employ the 



standard parameterization and took the values [f2 



sia = 0.220, s 23 = 0.039, s 13 = 0.0031, 5 13 = 70° 

Numerical results are shown in Figs. 2, 3, 4. 

In Fig. 2, we show the total cross section o to t of the process e + e~ — > bs as a function 
of the center-of-mass energy ^/s. There are three peaks, corresponding to the pole of 
neutral gauge boson Z°, a pair of charged gauge boson W threshold, and a pair of 
top quark tt threshold respectively. In most of high energy region, total cross section 
is the order of 10~ 3 fb, which is too small to be seen at LEP II or planning NLC 
colliders. Therefore, even a small number of bs events, detected at LEP II or NLC, 
will unambiguously indicate new FCNC couplings beyond SM. Smallness of the total 
cross section can easily be understood. One has 

E l M l 2 = e 8 | E V^V jc f(x j: yj)\ 2 

spins j=u,c,t 

- e \V tb V ts —^— — \ Xj , yj=0 + ...\ , (5) 

due to GIM mechanism, where Xj = rnj/m^yj = rrij/s, and denote the less im- 
portant terms for y's > 200 Gev. Assuming -^:\x-,y=o = 0(1); one obtains from eqs. 
(4), (5) 



a ~ 10~ 3 /& 



7 Bottom-Strange Associated Production at High Energy e + e Colliders in Standard 
Model 91 

at y/s = 200 Gev. 

We fixed the center-of-mass energy y/s at 200 Gev. Differential cross section of the 
process at the energy as a function of cos# is shown in Fig. 3. 
The Fig. 4 is devoted to the backward-forward asymmetry 

A _ Io /2 gg - g/g % d9 m 
Jo + J n/2 M dU 

as a function of y/s. 

To summarize, we have calculated the process e + e~ — > 6s in SM. We found that 
the total cross section is of the order of 10~ 3 fb in the high energy region which is still 
too small to be seen at LEP II or planning NLC. However, it is worth to note that 
the total cross section at Z resonance may reach as large as 10 _1 fb. Therefore, it is 
possible to see the process if a luminosity reaches 100-1000 fb -1 . In addition to that, 
the process is of a good place to search for new physics. 

7.4 Appendix 

91 = m s (B a (m 2 - m 2 )(m 2 - m 2 J(m 2 + 2m 2 J + B b m 2 (mt ~ 2mgm? + mj + m 2 m 2 w + m 2 m 2 w - 
2m 4 /) — Bo m b( m s ~~ 2m 2 m 2 + rrij + m 2 m 2 v + m 2 m 2 w — 2m/)){a\ — 4a2S^) + 
2m b m s {2C^ + Cl x m 2 b + C$m 2 + Cf(m 2 b + m 2 - 2m 2 w ) + Cf 2 s + C|s + 
Ci2( m ft ~ m l + s ))(°3 + 6a 4 s^, - 8a 4 s* ) - QC^m h m s (a 3 + 4a 4 s£,(<4 - s^)) + 
\2m h m s m 2 w {Ci + )(a 3 + 6a 4 c^,4 - 2a 4 4) (7) 

c/2 = -a 3 m| + 8a 4 m24-6m^(C 1 d 1 mg + C2 i 2S + C 1 d 2(mf-m2 + s))(a 3 + 8a 4 c2 4) + 

2/ D a/ 2 2\/ 2 2 \ , / 2 2 22 22, 4, 22 

m b m s {B {m b — m s ){m^ — m w ) + [m b m s — m b m^ — m s m^ + + m,-m w — 
2m%){B b - B c ) + m 2 w (B b (2m 2 - m 2 ) - B c (2m 2 - mg)))(ai + 6a 2 s 2 w - 4a 2 4) " 
6C$m 2 w (a 3 (m 2 b -m 2 + s) + 4a 4 s 2 w {(m 2 b - m 2 s ){c 2 w - s 2 w ) + 2sc 2 w )) - 6C$ (a 3 (m 2 + 6m 2 J 
Aa A s 2 w {c 2 w m 2 + I2c 2 w m 2 w - m 2 s 2 w )) + 2(2C e 00 + C e u m 2 + C e 22 s + C e 12 (m 2 - m 2 + s)) 
(a 3 (m 2 + 2m 2 /) + 4a 4 s£,(3m£, - 2m 2 s 2 w - Am^s 2 /)) + 2Clm 2 (a 3 (m 2 + m 2 - m 2 - 2m 2 /) 
2a 4 s 2 u {3m 2 - 3m 2 - 4m\s 2 w - ^m 2 s s 2 w + km 2 s 2 w + 8m^4)) + 2CI(a 3 s(m 2 + 2m 2 w ) - 
2a i s 2 l] (3m 2 m 2 — 3m 2 m 2 — Gm^s + 4m|ss^ + Sm^ss 2 /)) — 6Cfm^(a 3 (m^ — m 2 + s) + 
Aa A s 2 w {m 2 b c 2 w - m 2 b s 2 w - 2c 2 w m 2 s + 2c 2 w s)) - 6C$m 2 w (a 3 (s - m 2 - m 2 ) + 
Aa 4: s 2 v (2sc 2 u - ml + 2m 2 s 2 w - 2c 2 w m 2 )) + 2Cf (a 3 (2sm^ + m^m^ + m\m 2 - 2m 2 w m 2 /) + 
2a 4 sj|(3mjmj — 3m\m 2 + 6sm^, — Am1m 2 s 2 w — 4m 2 m 2 s w — 8sm^4 — 

+ 8m^m%)) (8) 



+ 



+ 



53 = -a 5 m 6 m s (2L>4 + Z? 3 ) + m s (B%(m 2 b - m 2 )(m 2 - m 2 J(m 2 + 2m 2 J + B b m 2 s (mt - 2m\m 2 + 
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m 4 j + m b m 2 u + m 2 m 2 u — 2m A w ) — B^m^(m\ — 2m 2 m 2 + ml- + m 2 m 2 + m 2 m 2 — 2m* ))(ai + 
2a 2 4 - 4a 2 4) + 2m 6 m s (2C e + C^mg + C e m 2 + Cf (mg + m) - 2m 2 J + Cf 2 s + Cfs + 

c i2( m b - m s + s ))( a 3 - 3a 4 + 10a 4 4 - 8a 4 4) ~ 6C# m 6 m s (a 3 - 2a 4 (l - 2s^) 2 ) + 
12m b m s m 2 w (C$ + C?)(a 3 - a 4 (3 - 4s* )(1 - 24)) (9) 

g 4 = -a 5 (2D f Q0 - (2D f 13 + D{)(m 2 - t) + 2D f 23 t + D f 3 t) - a 3 m 2 - 4a 4 mf4(l - 24) - 
6m* (C&mg + C d 22 s + Cf 2 (m 2 - m 2 + s))(a 3 - 4a 4 4(l - 24)) + 

2/na/ 2 2\/ 2 2 \ , ub/ 22 22 22, 4 , o 2 2 22, 

m b m s {B Q {m b — m s ){m^ — m w ) + B {m b m s — m b m^ — m s m^ + rrij + 2m b m w — m s m w + 

2 2 n 4 \ n c/ 22 22 22, 4 2 2 , n 2 2 22 n 4 \\ 

m^m w — 2m w ) — B {m b m s — m b rrij — m s m,j + rrij — m b m w + 2m s m w + rrijm w — 2m w j) x 

(oi - 3a 2 + 8a 2 4 - 4a 2 4) — QC2 m w( a 3( m b -m 2 + s) + 2a,i{c 2 v m 2 - <? w m\ - 2c 2 w s + 

m 2 b s 2 w + 2c 2 w m 2 b s 2 w - m 2 s s 2 w - 2c 2 w m 2 s s 2 w + 4c 2 w ss 2 w - 2m 2 b s A w + 2m 2 s s A w )) - 6C$ (a 3 (m 2 + 6m 2 ) + 

2a 4 (m 2 s 2 w - c 2 w m 2 - I2c 2 w m 2 w + 2c 2 w m 2 s 2 w + 2Ac 2 w m 2 w s 2 w - 2m 2 s A w )) + 

2(2C e + C^ml + C e 22 s){a 3 {m 2 + 2m 2 J + 2a 4 (2mf4 - 3m 2 w + Wm 2 w s 2 w - Am 2 si - Sm 2 w s 4 J) + 
2CQm 2 (a 3 (ml + rn 2 - rn 2 - 2m 2 w ) + a 4 (3m 2 - 3m 2 + 4m 2 4 + 10m s4 ~ 10m 2 4 ~ 8 "44 ~ 
8m|4 ~~ 8m 2 4 + 8m|4 + ^ m w s w)) + 2Cf (a 3 s(m 2 + 2m 2 u ) + a 4 (3m 2 m 2 — 3rn 2 rn 2 — 
6m 2 s — 6m 2 m 2 4 + 6m 2 m 2 s 2 u + 4m 2 ss 2 u + 20m^s4 — 8m 2 s4 — 16m^ss^,)) + 
2C 1 e 2 (m 2 - m 2 + s)(a 3 (m 2 + 2m 2 ) - 2a 4 (3m 2 - 2mJ4 - Wm 2 w s 2 w + 4m|4 + 8«44)) ~ 
6Cfm^,(a 3 (mg - + s) + 2a 4 (mgs^ - 3c 2 w m b + Qc 2 w m 2 b s 2 w - 2m^4 + 2 4* ~ 4 44* + 
2cljU - 444 n )) ~ 6Com^,(a 3 (s - m 2 - m 2 ) + 2a 4 (mg4 + 4 m fe _ ^ c w m b s w ~ ^fn 2 s 2 w - 
2m b s w + 4m^4 + 2m s4 _ 4m^44 ~ 2s 4 + 4«4 s w)) + 2Cf (a 3 (mfm^ + mfm^ + 
2m 2 m 2 w — 2rr? w t — 2m 2 w u) + a 4 (3mgm| — 3m 2 m 2 + 6m s m^ + 10mgm^4 — ^ m b m 'j s w ~ 
20m 2 m 2 u s 2 u — 8mgm^4 — 8mgm^4 + ISm^m^s^, — 6sm^, + 20sr?4s^ — 16sm^s^)) (10) 

g 5 = -24a 4 m s 4(Cfmf - 2C 2 d m2 ) - Am s (Cf im 2 b + Cgm| + C{{m\ + m^ 2 - 2m 2 ))(a 3 + 6a 4 4 - 
8a 4 4) - QrrisiC^ml + 2C$m 2 w + Cf{m 2 b - m 2 + 2m 2 ))(a 3 + 4a 4 44 - ^st) - 

d 2 2 2 22222 22 22 22 

QC 12 m s {a 3 (m b - rrij - 2m w ) + 4a 4 s w (m b c U) - c w m^ - 4:C w m w - m b s w + mjS w )) - 
ACt 2 m s {a 3 {m 2 b - m) - 2m 2 J + 2a 4 s 2 {3m 2 b - 6m 2 w - Am 2 b s 2 w + \m)s 2 w + 8m 2 4)) (H) 

56 = 24a 4 m b 4 (C|m 2 - 2C 2 d m 2 y ) - 12m b m 2 (C d + 2Cf)(o 3 + 8o 4 44) ~ 

4m fe m 2 (C^ + 2C{ )(a 3 - 8a 4 s^) - 6C^m 6 (a 3 (m 2 + 2m 2 ) + Aa A s 2 w {c 2 w m 2 + Ac 2 w m 2 w - m 2 s 2 w )) + 
6Cf 2 m fe (a 3 (m 2 - m 2 - 2m 2 y ) + 4a 4 s 2 (4 m s ~ 4 m j ~ 4 4 ra S; ~ TO s4 + m |4)) ~ 
4Cf 1 m fe (a 3 (m 2 + 2m 2 ) + 4a 4 4(3"i 2 - 2m 2 s 2 - 4m 2 s 2 )) + 4Cf 2 m b (a 3 (m 2 - m 2 - 2m 2 y ) + 
2a 4 4(3m 2 - 6m 2 - 4m 2 4 + 4mf4 + Sm^ 2 )) (12) 

g 7 = 2a 5 m s (D f 23 + D{) + 12a 4 m s (Clm 2 - 2C 2 V 2 y )(l - 24) - 4m s (C 1 e 1 m 2 + C^m 2 + C\(m 2 b + 
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m 2 — 2m 2 u ))(a 3 — 3a4 + lOa^ 2 — 8043^,) — 6m s (Cf 1 m^ + 2Com 2 + 

Cf(ml - m 2 + 2m 2 ))(a 3 - 2a 4 (4 - c^) 2 ) - 6Cf 2 m s (a 3 (ml - m| - 2m 2 w ) + 2a 4 (c 2 m 2 - 

2 2 , a 2 2 , 2 2 1 2 22 22 o 2 22 o 2 22 n 2 4 , o 2 4\\ 

c w m b + 4c w m w + m b s w + 2c w m b s w - m jSw - 2c w m jSw - 8c w m w s w - 2m b s w + 2m j sJ) - 
ACf 2 m s (a 3 (ml - m 2 - 2r?4) + a 4 (6m 2 - 3m 2 , + 10m|4 - 4m 2 s 2 - 20m 2 u s 2 - 8m^4 + 
8m 2 4 + 16m 2 4)) (13) 



g 8 = -2a 5 m b (D f 22 - d( 3 + D f 2 ) - \2a A m b {C e 2 m 2 - 2C 2 i m 2 w ){l - 2s 2 J - 12m b m 2 w (C$ + 2Cf)(a 3 - 
4a 4 4(l - 24)) - 4m b m 2 (C e + 2C 1 e )(a 3 + 4a 4 4(l " 24)) " 6C 1 d 1 m 6 (a 3 (m 2 + 2m 2 ) - 

o i 2 2 , a 2 2 2 2 o 2 22 o 2 2 2 , o 2 4 \\ , i i 2 2 

2a i {c w m j + 4c w m^ - rn^-s^ - 2c w m^s w - 8c w m w s w + 2m ? s w J) + 6G 12 m fe (a 3 (m s - m^ - 

n 2 \ ry / 2 2 2 2 a 2 2 22 o 2 2 2 i 2 2 i 2 2 2 i o 2 2 2 i 

^^w ) Z(X^.yC w 7Tl s CyjTTlj ^±CyjTflyj ???< s Syj ZC^TTt^S^ ~~\~ 7TtjS w ~~\~ ZC^TTljS^ ~\~ oC u , TH/^y S, u , ~\~ 

2m 2 ,sf u — 2m 2 4)) — 4Cf 1 m;,(a 3 (m 2 + 2m 2 u ) + 2a4(2m 2 4 — 3m^, + lOm^s^ — 4m 2 s^ — 
8m 2 4)) + ^Cf 2 m b (a 3 (m 2 - m 2 - 2m 2 w ) + a 4 (6m 2 - 3m, 2 + 10m 2 4 — 4m 2 s 2 - 
20m 2 w s 2 w - 8m 2 4 + 8mf4 + 16m 2 4)) (14) 
g 9 = 2a 5 (D f 12 -2D f 13 + D f 22 -D f 23 + D f 2 ) (15) 
g 10 = -a 5 m s (2D f 13 + 2D f 23 + d{) (16) 
5ii = a 5 m b (2D{ 3 + D{) (17) 

where is defined as 

1 1 

°1 — TTT^ o 9 — 77 — 5 9\ 9~' a 2 



1927r 2 sm; ) m 2 m 2 u (m 2 . — m^sj, ' 7Q&Tt 2 m b m 2 s m^{m 2 — m 2 ){m 2 — im z Yz — s)c 2 IJ s 

1 



4 



a 3 — ^„ 9 9 9 > °4 — „„„ 9 9/9 • m \ 9 A 1 a 5 



967r 2 ,sm 2 s 2 7§%ii 2 m 2 w {m 2 z — im z Tz — s)c 2 ll s^ u 32ir 2 sf v 

where c w = cos9 w and s w = sin6 w . In the presentation of gj above, we have used the 
definition of scalar integrals Bs, Cs,and Z)g||13||,and these functions, Bs, Cs, and Ds, 
with superscripts a,b,...,f have the arguments 

(0,m 2 ,m 2 ), (m 2 b ,m 2 ,m 2 w ), (m 2 ,m 2 ,m 2 w ), (m 2 b ,m 2 , s,m 2 w ,m 2 ,m 2 J 
(ml, m 2 , s, m 2 , m 2 , m 2 ), (0, m 2 , m 2 , 0, t, s, 0, m 2 w , m 2 , m 2 ) 

respectively. Here rrij denotes the mass of up-type quark u, c, t. 
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Figure 1: Typical Feynman diagram of prosess e + e — > bs 
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Figure 2: Cross section of the process e + e — > 6s as a function of y^. 
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Figure 3: Differential cross section of the process e + e — > bs, where = 200 GeV. 
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Figure 4: Afb of the process e + e — > bs as a function of y/s. 
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8 B — ► X s r + r in a CP softly broken two Higgs dou- 
blet model 

ABSTRACT 

The differential branching ratio, forward-backward asymmetry, CP asymmetry and 
lepton polarization for a B-meson to decay to strange hadronic final states and a 
t + t~ pair in a CP softly broken two Higgs doublet model are computed. It is 
shown that contributions of neutral Higgs bosons to the decay are quite significant 
when tan (3 is large. And it is proposed to measure the direct CP asymmetry in 
back-forward asymmetry. 



8.1 Introduction 

The origin of the CP violation has been one of main issues in high energy physics since 
the discovery of the CP violation in the Kq — K systerm in 1964 [jlj. The measurements 
of electric dipole moments of the neutron and electron and the matter-antimatter 
asymmetry in the universe indicate that one needs new sources of CP violation in 
addition to the CP violation come from CKM matrix, which has been one of motivations 
to search new theoretical models beyond the standard model (SM). 

The minimal extension of the SM is to enlarge the Higgs sectors of the SM . It 
has been shown that if one adheres to the natural flavor conservation (NFC) in the 
Higgs sector, then a minimum of three Higgs doublets are necessary in order to have 
spontaneous CP violations ||. However, the constraint can be evaded if one allows 
the real and image parts of </>i~</>2 have different self-couplings and adds a linear term 
of Re(0^0 2 ) in the Higgs potential (see below Eq. (0) with ml = 0). Then, one can 
construct a CP spontaneously broken two Higgs doublet (2HDM), which is the minimal 
and the most "economical" one f\ among the extensions of the SM that provide new 
source of CP violation. Furthermore, in addition to the above terms, if one adds a 
linear term of lm(0f $2), then one has a CP softly broken 2HDM |4]]. 

Flavor changing neutral current (FCNC) transitions B — > X s ^ and B — > X s l + l~ 
provide testing grounds for the SM at the loop level and sensitivity to new physics. Rare 
decays B — > X s l + l~(l = e,/i) have been extensively investigated in both SM and the 
beyond || [?J. In these processes contributions from exchanging neutral Higgs bosons 
(NHB) can be safely neglected because of smallness of -^-{l = e, //). The inclusive 
decay B — ► X s t + t~ has also been investigated in the SM, the model II 2HDM and 
SUSY models with and without including the contributions of NHB [|], |], [K| |TTJ . In this 



note we investigate the inclusive decay B — > X s t + t~ with emphasis on CP violation 
effect in a CP softly broken 2HDM, which we shall call Model IV hereafter for the 
sake of simplisity. We consider the Model IV in which the up-type quarks get masses 
from Yukawa couplings to the one Higgs doublet H 2 and down-type quarks and leptons 



2 Comparing the Model III 2HDM in which CP is explicitly violated, the CP spontaneously 
broken 2HDM has only two new parameters besides the masses of the Higgs bosons in the large tan (3 
limit (see below) . In this sense it is the most " economical" . 
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get masses from Yukawa couplings to the another Higgs doublet H\. The Higgs boson 
couplings to down-type quarks and leptons depend on only the CP violated phase £ 
which comes from the expectation value of Higgs and the ratio tg/3 = ^ in the large 
tg/3 limit (see next subsection), which are the free parameters in the model. Because 
the couplings of the charged Higgs to fermions in Model IV are the same as those in 
the model II, the constraints on tan f3 due to effects arising from the charged Higgs are 
the same as those in the model II. Constraints on tg(3 from K — K and B — B mixing, 
T(b — > sj),T(b — > ctu t ) and Rb have been given 

°-W<0.52(^|) (1) 

(and the lower limit m#± > 200Gev has also been given in the ref. [0). It is obvious 
that the contributions from exchanging neutral Higgs bosons now is enhanced roughly 
by a factor of tg 2 (5 and can compete with those from exchanging 7, Z when tg/3 is large 
enough. Because the CP violation effects in B — > X s t + t~ come from the couplings of 
NHB to leptons and quarks, we shall be interested in the large tan f3 limit in this note. 
The constraints on £ can be obtained from the electric dipole moments (EDM) of the 
neutron and electron, which will be analysed in the next subsection. 



8.2 Model description 

Consider two complex y — 1, SU{2) W doublet scalar fields, <pi and 2 . The Higgs 
potential which spontaneously breaks SU(2) x U(l) down to U(1)em can be written 
in the following form M: 

i=l,2 

+m\Re{(j) J [ 4> 2 ) + m|Jm(^0 2 ) 

+a 3 [(0M)(^ 2 )] + U(4>M(<Pt<Pi)\ 

+A 5 [Re(0+0 2 )] 2 + A 6 [Im(0+0 2 )] 2 (2) 

Hermiticity requires that all parameters are real. The potential is CP softly broken 
due to the presence of the term m|/m(05 l "0 2 ). It is easy to see that the minimum of 
the potential is at 

thus breaking SU(2) x £7(1) down to U(1)em and simutaneously breaking CP, as 
desired. It should be noticed that only for A5 7^ Ae, the phase £ can't rotated away 
as usual, which breaks the CP-conservation. If 777,4=0 in (2) then the potential is CP 
invarint. It has been shown that the CP spontaneously breaking happens at (3) Hl3| . 
We limit ourself to the case of ml ^ in the paper and shall investigate the 7774=0 



case in a separate paper 14 



8 B — > X s r + r in a CP softly broken two Higgs doublet model 101 



m ..!\-, -^jr 1, <4) 



In the following we will work out the mass spectrum of the Higgs boson. For charged 
components, the mass-squared matrix for negtive states is 

Diagonalizing the mass-squared matrix results in one zero-mass Goldstone state: 

G~ = sin (3 fa + cos (3 fa , (5) 

and one massive charged Higgs boson state: 

H~ = cos/902 — sin (3 fa, (6) 
m H - = \K\{vl + v 2 2 ) 1 (7) 

where tan/3 = v-ijv\. Correspondingly we could also get the positive states G + and 
H + with the same masses zero and | A4 1 (f 1 + respectively. 

For neutral Higgs components, because CP-conservation is breaking, the mass- 
squared matrix is 4 x 4, which could not be simply separated into two 2x2 matrices 
as usual. However, in the case of large tan (3 which is we intrested in, the neutral parts 
can be written as separately two 2x2 matrices and one of them is 

/ A 5 +A 6 +(A6-A 5 )cos(2g) (A B -A 5 )sin(2g) \ 

V 2 I (A 6 -A 2 5 )sin(2g) A 5 +A 6 +(A 5 ~A 6 ) cos(2g) j ■ ( 8 ) 

V 2 2 / 

Diagonalizing the Higgs boson mass-squared matrix results in two eigenstates: 







H° 2 




(9) 



with masses 



m H o = X 5 v 2 

m 2 H o = X 6 vl, (10) 
where = cos£ and = sin^. The diagonalizing of the 4x4 neutral Higgs mass- 



squared matrix has been analytically carried out under some assumptions in Ref. [15] 
and the results reduce to Eq. (|) and (|i~Cf ) in the case of large tan/3. 

The another 2x2 matrix can be similarly deal with. Because the couplings of the 
third physical neutral Higgs boson and neutral Goldstone to down-type quarks and 
leptons are not enhanced for large tan f3 case in which we are interested, we do not 
show the explicit results. 

Now, we turn to the discussion of the Higgs-fermion-fermion couplings. After com- 
pleting the transformation from the weak states to the mass states, the couplings of 
neutral Higgs to fermions which are relevant to our analysis are 
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where / represents down-type quarks and leptons. And the couplings of the charged 
Higgs bosons to fermions are the same as those in the CP-conservative 2HDM (model 
II, for examples see Ref. This is in contrary with the model III in which the 

couplings of the charged Higgs to fermions are quite different from model II. It is easy to 
see from Eq. (|TID that the contributions come from exchanging NHB is proportional to 
^/^Gps^c^rrtf/ cos 2 /?, so that the constaints due to EDM translate into the constraints 
on sin 2£ tan 2 (3 (1/ cos (3 ~ tan (3 in the large tan (3 limit). According to the analysis in 
Ref. []T7|], we have the constraint 



sin 2£ | tan/? < 50 (12) 

from the neutron EDM. And the constraint from the electron EDM is not stronger than 
Eq. ([12]). It is obvious from Eq. that there is a constraint on £ only if tan f3 > 50 
and the stringent constraint on tan/5 comes out and is tan/3 < 50 when £ = 7r/4. 

8.3 Formula for B -> X s t + t~ 

Inclusive decay rates of heavy hadrons can be calculated in heavy quark effective theory 
(HQET) ]n| and it has been shown that the leading terms in I /nig expansion turn out 
to be the decay of a free (heavy) quark and corrections stem from the order 1/toq [19[. 
In what follows we shall calculate the leading term. The transition rate for b — > st + t~ 
can be computed in the framework of the QCD corrected effective weak hamiltonian, 
obtained by integrating out the top quark, Higgs bosons and W ± , Z bosons 

ATI 10 10 

Heff = -^V*V2Ea(A*)0i(A*) + Y, C QMQi{ti) (13) 

i=i j=i 

where Oi(i = 1, • • • , 10) is the same as that given in the ref. |J, Qi's come from exchang- 
ing the neutral Higgs bosons and are defined in Ref. [plj. The explicit expressions of 
the operators governing B — > X s t + t~ are given as follows: 

7 = (e/16n 2 )m b (s La a^b Ra )F^, 

8 = (e/167r 2 )(s La 7 / "6 La )f7 M r, 

9 = (e/167r 2 )(s La 7 / "6L Q )f7 A1 75r, 
Qt = (e 2 /lQ7T 2 )(s La b Ra )(fr), 

Q 2 = (e 2 /l6ir 2 )(s La b Ra )(f l5 r). (14) 

At the renormalization point /i = the coefficients Cj's in the effective hamilto- 
nian have been given in the ref. and Cq/s are (neglecting the 0(tg(3) term) 

c Ql (m w ) = m ; m :^ { e ^(m+m)}, 

m b m T tg 2 (3x t ^ , , 

°^ mw) = 2*in*9 w { J? M rf iflBi + f2Ei)h 
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where 



2 

C Qa (m w ) = —^-(C Ql (m w )+C Q2 (m w )), 
m T gj 

2 

C Q4 (m w ) = —^-(C Ql (m w )-C Q2 (m w )), 
m T gj 

C Qi {m w ) = 0, z = 5,...,10 (15) 



= o, 


i = 5, - ■ • , 10 






= -s^ D Hl = 






= % D H2 = 




Bh 1 


- 2 , BH 2 - 


c ? + zs ? 
2 



/l 

/a 



x t lnx t %H±lnxn± — Xtlnxt 



x t 



x H ± - x t 



x t lnx t 



XH±lnxu± 



(x t - l)(x H ± - 1) (x H ± - x t )(x H ± - I] 



with Xi = mf/m^. In Eq. (|T5|) , i?j are given by 



(16) 



Eh 2 

Cl 

c 2 



-(-S^Ci +c 5 c 2 ), 
\{c i c 1 + s^c 2 ), 

—x H ± + C£X Hl (c£ + is{) + S£Xh 2 (s£ - ic{), 
i {—x H ± + S£X Hl (s£ - ice) + C£Xh 2 (c£ + is^)) . 



(17) 



Neglecting the strange quark mass, the effective hamiltonian ( |T3"D leads to the fol- 
lowing matrix element for b — > st + t~ 



M 



GpOL 

71 



-V tb V: s [Ct fS s Ll ,b L f^r + C 9 s L j„b L f^j 5 r 



+ 2C 7 m b s L ia» v ^b R T^T + C Ql s L b R rr + Cq 2 s l & r t7 5 t 



;is) 



where || ||, [20 



U 8 



+ 



{g( — ,s) 

m b 



c 8 

3 



^ , , „ M$ - q 2 - lM Vi V Vi 



(19) 



with s = q /ml, q = (p T + + p T -) ■ I n (jjj) fl , (^ £ >^) arises from the one-loop matrix 



element of the four-quark operators and can be found in Refs. || ^TJ. The second 
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term in braces in (|T9|) estimates the long-distance contribution from the intermediate, 



ip', ip" ■■■ H In our numerical calculations, we choose fc(3Ci + C 2 ) = —0.875 p2 |. 

The QCD corrections to coefficients Cj and C , Q i can be incooperated in the stan- 
dard way by using the renormalization group equations. Although the Cj at the scale 
/i = 0(rrib) have been given in the next-to-leading order approximation (NLO) and 
without including mixing with Qi, we use the values of Ci only in the leading order 
approximation (LO) since no Cq { have been calculated in NLO. The Ci and Cq. with 
LO QCD corrections have been given in Ref. [|i~0[ |. 



C 7 (m») = rT 16/23 

-0.012C Q3 (m w )], (20) 



C 8 (m b ) = C 8 (m w )+^^[-4(l-^ 11/23 ) + ^(l-^ 29/23 )]C 2 (m w ),(21) 

a s (m w ) 33 87 

C 9 (m b ) = C 9 (m w ), (22) 
C Qi (m 6 ) = !)-^ /A C ft K), z = l,2, (23) 

where 7q = —4 JZ^] is the anomalous dimension of sl&r, /3 = 11 — 2nj/3, and 77 = 
a s (m fe )/a s (m w ). 

After a straightforward calculation, we obtain the invariant dilepton mass distribu- 
tion 

dr(g ^ +T ~> = ^^ n / n (i-.«i-*)« DW 

ds Aix 2 f{m c /m b ) s |l/ cfe | 2 

9+2 oj.2 o 

= |C 8 e// | 2 (l + — )(1 + 2 S )+4|C 7 | 2 (1 + — )(1 + -) 
s s s 

9+2 9+2 

+ |C 9 | 2 [(1 + 2s) + -y (1 - 4s)} + 12Re(C 7 C 8 e// *)(l + =-) 
+^ Ql | 2 (s - 4t 2 ) + ||C Q2 | 2 S + 6Re(C 9 C* 2 )t (24) 

where s=g 2 /m 2 , t=m T /m&, 5(5 — > X C W) is the branching ratio, / is the phase-space 
factor and f(x)=l — 8x 2 + 8x 6 — x 8 — 24x 4 In x. 

The CP asymmetry for the B — > XJ + l~ and £> — > X S Z + Z~ is defined as 

We also give the forward-backward asymmetry 



1 JJ*fi+A*£~J>M 
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where z = cos 6 and 6 is the angle between the momentum of the B-meson and that of 
l + in the center of mass frame of the dileptons t + t~. Here, 

E{s) = Re(C 8 e// C*s + 2C 7 C* 9 + C e 8 ff C* Q1 t + 2C 7 C* Q2 t). (27) 

The CP asymmetry in the forward-backward asymmetry for B — ► X s t + t~ and B — > 
X s t + t~ is defined as 

Acp(S) " A( S ) + A( S )- (28) 

It is easy to see from Eq. (|24|) that the CP asymmetry is very small because the 
weak phase difference in CVC| arises from the small mixing of O7 with Q 3 (see Eq. 
(20|)). In contrast with it, A 2 CP can reach a large value when tan/? is large, as can be 



seen from Eq. ( p?|) and (0). Therefore, we propose to measure A 2 CP in order to search 
for new CP violation sources. 

Let us now discuss the lepton polarization effects. We define three orthogonal unit 
vectors: 



(-L 



e N 



Pi 

\Pi\ ' 

Ps X pi 



\Ps X pi I 

e T = e N x e L , 

where pi and p s are the three momenta of the i~ lepton and the s quark, respectively, 
in the center of mass of the £ + l~ system. The differential decay rate for any given 
spin direction n of the i~ lepton, where n is a unit vector in the £~ lepton rest frame, 
can be written as 

^T = K^) [ 1 + ^ LiL + PNSN + PT ^ ) - fi \ ' (29) 

where the subscript "0" corresponds to the unpolarized case, and Pl, Pt, and P/v, 
which correspond to the longitudinal, transverse and normal projections of the lepton 
spin, respectively, are functions of s. From Eq. fl29|), one has 

dT dY 

_ ( n = a) - — in = -ej 

— [n = a) + — (n = -ej 
as as 

The calculations for the Pj's (i = L, T, iV) lead to the following results: 

Pl = (1 _^^£>, 

s -D(s) 

= 37T (1 4t 2 )1/2 Ay(s) 



4sV2 v s 7 D(s) ' 

3vrt £> T (s) v 

Pt - "2^W (31) 
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2 ' 



where 

D L (s) = Re(2(l + 2s)C e 8 ff C* 9 + 12C 7 C* 9 -QtC Ql C* 9 -3sC Ql C* Q 
D N (s) = Im (2sC Ql C; + sC Ql Ci ff * + sC Q2 C; + AtC 9 C* 7 + 2tsC e 8 ff *C 9 

D T {s) = Re (-2C 7 C* 9 + AC e 8 ff C* 7 + ^\C 7 \ 2 - C e 8 ff C* 9 

-LeirW./'i 2 s ~ ^ n r* s n r* s n e ffn* I (io\ 

Pj (i=L, T, N) have been given in the ref. [9], where there are some errors in P? and 
they gave only two terms in D^, the numerator of P^. We remind that P^ is the 
CP-violating projection of the lepton spin onto the normal of the decay plane. Because 
Pn in B — > X s l + l~ comes from both the quark and lepton sectors, purely hadronic 
and leptonic CP-violating observables, such as d n or d e , do not necessarily strongly 



constrain P N [231. So it is advantageous to use Pn to investigate CP violation effects 



in some extensions of SM [p5|1 . In the model IV 2HDM, as pointed out above, d n and 



d e constrain J\ sin2£| tan/? and consequently Pn through CQ i (i = 1,2) (see Eq. 
8.4 Numerical results 

The following parameters have been used in the numerical calculations: 

m t = 175Gev, nib = 5.0Gev, m c = 1.6Gev, m T = 1.77Gev, r] = 1.724, 

rriH 1 = lOOGev, rriH 2 = tuh± = 200Gei;. 

Numerical results are shown in Figs. 1-9. From Figs. 1 and 2, we can see that 
the contributions of NHB to the differential branching ratio dT/ds are significant when 
tan p is not smaller than 30 and the masses of NHB are in the reasonable region, and 
the forward-backward asymmetry A(s) is more sensitive to tan (3 than dT/ds, which is 
similar to the case of the normal 2HDM without CP violation [IIJ . 



The direct CP violation A l CP (i = 1, 2) and CP-violating polarization P^ of B — * 
X s t + t~ are presented in Figs. 3-7, respectively. As expected, A l CP is about 0.1% and 
hard to be measured. However, A 2 CP can reach about 10%. A 2 CP is strongly dependent 
of the CP violation phase £ and comes mainly from exchanging NHBs as expected. 
From Figs. 6 and 7, one can see that Pn is also strongly dependent of the CP violation 
phase £ and can be as large as 5% for some values of £, which should be within the 
luminosity reach of coming B factories, and comes mainly from NHB contributions in 
the most of range of £. 

Figs. 8 and 9 show the longitudinal and transverse polarizations respectively. It 
is obviously that the contributions of NHB can change the polarization greatly, espe- 
cially when tan/5 is large, and the dependence of Pi on CP violation phase £ is not 
significant in the most of range of £. The longitudinal polarization of B — * X s t + t~ has 
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been calculated in SM and several new physics scenarios ||. Switching off the NHB 
contributions, our results are in agreement with those in Ref. |§. 

In summary, we have calculated the differential braching ratio, back-forward asym- 
metry, lepton polarizations and some CP violated observables for B — ► X s t + t~ in the 
model IV 2HDM. As the main features of the model, NHB play an important role in 
inducing CP violations, in particular, for large tan (3. We propose to measure J\^p 7 the 
direct CP asymmetry in back-forward asymmetry, in stead of A^ P , the usual direct 
CP violation in branching ratio, because the former could be observed if tan (3 is large 
enough (say, > 30) and the latter is too small to be observed. It is possible to discrim- 
inate the model IV from the other 2HDMs by measuring the CP-violated observables 
such as Aq P , Pjv if the nature chooses large tan (5. 
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Figure 1: Differential branching ratio as function of s, where £ = tt/4, solid and dashed 
lines represent tan j3 — 10 and 30, dot-dashed line represents the case of switching off 
CQ i contributions. 
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Figure 2: Backward-forward asymmetry as function of s, where £ = 7r/4, solid and 
dashed lines represent tan /3 — 10 and 30, dot-dashed line represents the case of switch- 
ing off Cq i contributions. 
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Figure 3: A^ p as function of £, where s = 0.8, solid and dashed lines represent 
tan (3 = 10 and 30. 
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Figure 4: A 2 CP as function of £, where s = 0.8, solid and dashed lines represent 
tan/5 = 10 and 30. 
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Figure 5: A 2 CP as function of s, where £ = 7r/4, solid and dashed lines represent 
tan/5 = 10 and 30. 
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Figure 6: P N as function of s, where £ = vr/4, solid and dashed lines represent tan (3 = 
10 and 30, dot-dashed line represents the case of switching off C Qi contributions. 
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Figure 7: P N as function of £, where s = 0.8, solid and dashed lines represent tan/5 = 
10 and 30, dot-dashed line represents the case of switching off C Qi contributions. 
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Figure 8: P L as function of £, where s = 0.8, solid and dashed lines represent tan/5 = 
10 and 30, dot-dashed line represents the case of switching off Cq 1 contributions. 
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Figure 9: P T as function of £, where s = 0.8, solid and dashed lines represent tan/5 = 
10 and 30, dot-dashed line represents the case of switching off Cq 1 contributions. 
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9 Exclusive Semileptonic Rare Decays B — ► (K, K*)£ + £ 
in Supersymmetric Theories 

ABSTRACT 

The invariant mass spectrum, forward-backward asymmetry, and lepton polariza- 
tions of the exclusive processes B — > K(K*)£ + £~ , I = /i, r are analyzed un- 
der supersymmetric context. Special attention is paid to the effects of neutral 
Higgs bosons (NHBs). Our analysis shows that the branching ratio of the pro- 
cess B K/i + fj," can be quite largely modified by the effects of neutral Higgs 
bosons and the forward-backward asymmetry would not vanish. For the process 
B — ► K*^t + ^r 7 the lepton transverse polarization is quite sensitive to the effects 
of NHBs, while the invariant mass spectrum, forward- backward asymmetry, and 
lepton longitudinal polarization are not. For both B — > Kt + t~ and B — > K*t + t~ , 
the effects of NHBs are quite significant. The partial decay widths of these pro- 
cesses are also analyzed, and our analysis manifest that even taking into account 
the theoretical uncertainties in calculating weak form factors, the effects of NHBs 
could make SUSY shown up. 



9.1 Introduction 

The inclusive rare processes b — > X s £ + £~,£ — e, //, r have been intensively studied in 
literatures |TJ, 0, [| [|, |5|, |5], |7[ P, |IU|, |TT| , As one of flavor changing neutral current 
processes, it is sensitive to fine structure of the standard model and to the possible new 
physics as well, and is expected to shed light on the existence of new physics before 
the possible new particles are produced at colliders. 

It is well known that invariant mass spectrum, forward-backward asymmetries, and 
lepton polarizations are important observables to probe new physics, while the first 
two observables are mostly analyzed. About lepton polarizations, it is known that due 
to the smallness of the mass of it, therefore electron polarizations are very difficult 
to be measured experimentally. So only the lepton polarizations of muon and tau 
are considered in literatures [III, [T2|, [T3], |TJj]. The longitudinal polarization of tau in 



B — > X s t + t has been calculated in standard model (SM) and several new physics 
scenarios F° r B — > X s l + l~ (I = /U,r), the polarizations of lepton in SM are 



analyzed in |12| and it is pointed out that for the \i channel, the only significant 
component is Pl, while all three components are sizable in the r channel. The analysis 
has been extended to supersymmetric models (SUSY) and a CP softly broken two 
Higgs doublet model in refs. and |14| respectively. The reference also gives 

a general model-independent analysis of the lepton polarization asymmetries in the 
process B — > X s t + t~ and it is found that the contribution from Clrlr + Clrrl is 
much larger than other scalar-type interactions. 

Compared with the inclusive processes b — > X s £ + £~,£ = e, /i, r, the theoretical 
study of the exclusive processes B —> K(K*)£ + £~ is relatively hard. For inclusive 
semileptonic decays of B, the decay rates can be calculated in heavy quark effective 
theory (HQET) fTS] . However, for exclusive semileptonic decays of B, to make theoret- 
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ical predictions, additional knowledge of decay form factors is needed, which is related 
with the calculation of hadronic transition matrix elements. Hadronic transition ma- 
trix elements depend on the non-perturbative properties of QCD, and can only be 
reliably calculated by using a nonperturbative method. The form factors for B decay 
into have been computed with different methods such as quark models |f[6fl , SVZ 
QCD sum rules |L7]], light cone sum rules (LCSRs) [HJ |2C], [21|, |2"2fl . Compared 



to the lattice approach which mainly deal with the form factors at small recoil, the 
QCD sum rules on the light-cone can complement arily provide information of the form 
factors at smaller values of s. And they are consistent with perturbative QCD and the 
heavy quark limit. In this work, we will use the weak decay form factors calculated by 
using the technique of the light cone QCD sum rules and given in p3| . 

A upper limit on the branching ratio of B° — > K°*fi + fi~ has been recently given by 
CLEO || : 

BR(B° -»■ K°*fM + fj,-) < 4.0 x KT 6 , (1.1) 
and they will be precisely measured at B factories, these exclusive processes are quite 



worthy of intensive study and have attracted many attentions | 23j , |25 , 2q , 27 , |28 



2], [33], [34], In reference p3| , by using improved theoretical calculations of the 
decay dorm factors in the light cone QCD sum rule approach, dilepton invariant mass 
spectra and the forward-backward asymmetry of these exclusive decays are analyzed in 
the standard model and a number of popular variants of the supersymmetric models. 
However, as the author claimed, the effects of neutral Higgs exchanges are neglected. 



For exclusive processes, as pointed out in [35], the polarization asymmetries of \x and 



t for B — > K^*)fj, + fx and B — > K^*)t + t are also accessible at the B-factories under 
construction. In reference [Q, the lepton polarizations and CP violating effects in 



B — > K*t + t~ are analyzed in SM and two Higgs doublet models. 

As pointed in refs. [|, f4|, in two-Higgs-doublet models and SUSY models, neutral 
Higgs boson could contribute largely to the inclusive processes b — > X s £ + £~,£ — fi, r 
and greatly modify the branching ratio and forward-backward asymmetry in the large 
tan/3 case. The effects of neutral Higgs in the 2HDM to polarizations of r in B — ► 
Kt + t~ are analyzed in [ p3| , and it was found that polarizations of the charged final 
lepton are very sensitive to the tan/?. 

In this paper, we will investigate the exclusive decay B — > K(K*)£ + £~~ , £ = /i, r in 
SUSY models. We shall evaluate branching ratios and forward-backward asymmetries 
(FBA) with emphasis on the effects of neutral Higgs and analyze lepton polarizations 
in MSSM. According to the analysis of [ESJ, different sources of the vector current could 
manifest themselves in different regions of phase space, for the very low s the photonic 
penguin dominates, while the Z penguin and W box becomes important towards high s. 
In order to search the regions of s where neutral Higgs bosons could contribute large, we 
analyze the partial decay widths of these two processes. Beside that they are accessible 
to B factories, our motivation also bases on the fact that to the inclusive processes 
B — > X s £ + £~ , £ = /i, r, neutral Higgs could make quite a large contributions at certain 
large tan/? regions of parameter space in SUSY models, since part of supersymmetric 
contributions is proportional to tan 3 /? [§|. Such regions considerably exist in SUGRA 
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and M-theory inspired models |36| . We also analyze the effects of neutral Higgs to the 
position of the zero value of the forward-backward asymmetry. Our results show that 
the branching ratio of the process B —>■ Kfi + fi~ can be quite largely modified by the 
effects of neutral Higgs bosons and the forward-backward asymmetry would not vanish. 
Because the FBA for B —>■ K £ + 1~ {£ = /i, r) vanishes if there are no the contributions 
of NHBs and the contributionsof NHBs can be large enough to be observed only in 
SUSY and/or 2HDM with large tan/3, a non-zero FBA for B — > K £ + £~ would signal 
the existance of new physics. For the process B — > K*fi + {i~, the lepton transverse 
polarization is quite sensitive to the effects of NHBs, while the invariant mass spectrum, 
forward- backward asymmetry, and lepton longitudinal polarization are not. For both 
B — > Kt + t~ and B —>■ K*t + t~ , the effects of NHBs are quite significant. Our analysis 
manifest that even taking into account the theoretical uncertainties in calculating weak 
form factors, the effects of NHBs could show SUSY up. In a word, our analysis manifest 
that effects of NHBs is quite remarkable in some regions of parameter space of SUSY, 
even for the process B — > K\i" 'f/T . 

The paper is organized as follows. In the subsection 2, the effective Hamiltonian is 
presented and the form factors given by using light cone sum rule method are briefly 
discussed. Basic formula of observables are introduced in subsection 3. Section 4 is 
devoted to the numerical analysis. In subsection 5 we make discussions and conclusions. 



9.2 Effective Hamiltonian and Form Factors 

By integrating out the degrees of heavy freedom from the full theory, MSSM, at elec- 
troweak(EW) scale, we can get the effective Hamiltonian describing the rare semilep- 
tonic decay b — > s£ + £~ : 



Ti 



eil' 



4G F 

'71" 



10 



10 



tW(ECi(A*)O<(A0 + Y,CQ t (n)Qi<Ji)) 



(2.1; 



i=l i=l 

where the first ten operators and Wilson coefficients (WC) at EW scale can be 
found in || [37|[], and last ten operators and WC which represent the contributions of 
neutral Higgs can be found in [§]]. 

With the renormalization group equations to resum the QCD corrections, WCs at 
energy scale \x = m& are evaluated. Theoretical uncertainties related to renormalizat ion- 
scale can be substantially reduced when the next-leading-logarithm corrections are 



included [38 



The above Hamiltonian leads to the following free quark decay amplitude: 



M(b^ s£ + £-) 



V*V tb [C 9 ^ [s^Lb] 



hH +C 10 [s^Lb] £^£ 



-2m b G 



off 



sia 



-Rb 



+ Cqi [sRb] 



C Q2 [sRb] 




1 In our previous papers, e.g., ||, ||, we follow the convention of ref. [|lj for the indices of operators 
as well as Wilson coefificients. In this paper, we use more popular conventions (see, e.g., @). That 
is, 8 — * Og and Og — > O w . 
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where Cg is defined as |39|, HO 



3vr 



C 9 eff (/i,S) = C 9 ( yU ) + F(/i,S) + -C( /U ) £ .., , . . 



\"\. (2.3) 



where s = s/m 2 b ,s=q 2 , C(p) = (3 d + C 2 + 3 C 3 + C 4 + 3 C 5 + C 6 ), and 



-~<?(1, s) (4 C 3 + 4 C 4 + 3 C 5 + C 6 ) - ~s(0, 3) (C 3 + 3 C 4 ) 
-^(3C 3 + C 4 + 3C 5 + C 6 ) . 



(2.4) 



where the function g(rh c , s) comes from one loop contributions of four-quark operators 
and is defined by 



4, 2 8 16 z 2 

--lnz 2 -\ 1 

9 27 9 a 



4 /M _ 1(2 + ^Wctan 



Az 2 < s 

Az 2 > s 
(2.5) 

The last terms in (|2.3|) are nonperturbative effects from (cc) resonance contributions, 
while the phenomeno logical factors Ki can be fixed from the processes [23] B — > 
ifWy. _^ K { *k + £- and as given in the Table. |. 

Exclusive decays -B — > (If, K*)£ + £~ are described in terms of matrix elements of the 
quark operators in Eq. ( |2.2j ) over meson states, which can be parametrized in terms of 
form factors. 

For the process B — > K£ + £~ , the non- vanishing matrix elements are (q = p B — p) 



m B — m K 



% + 



(K(p)\s^b\B(p B )) = f + (s)Up B +p), ^ 
and 

{K(p)\sa^q v {l+ lh )b\B(p B )) = (K{p)\so^b\B{p B )) 



m B -m* K 



fo(s)q„ (2.6) 



{(Pb +p), 



s ~ q»(m B - m 2 K 



)} 



777. _B + 



(2.7) 



K 


J/tf 


*' 


K 

K* 


2.70 
1.65 


3.51 
2.36 



Table 1: ifytfce factors in B -> ifMj/*,*' -> K^£ + t decays calculated using the 
LCSR form factors. 
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While for B — > K*£ + £ , related transition matrix elements are 

(K*(p)\(V-A) ll \B(p B )) = -ze;(m B + m^A^s) + i(p B + p),(e*p B ) ^ 

^ r m B + m K , 

+iq»(e* PB ) — - A (s)) + ^/ 2 ^ (s) . (2.8) 

and 

(^1*7^(1 + 75 )6|5(p B )> = k^e^p^^T^s) 

+ T 2 (s) {e;(m| - m^.) - (e*p B ) (p B + p),} 

+ T 3 (s)(e*p B ) L - -5-^-5- (Pfl +P) ,} (2-9) 

Where e M is polarization vector of the vector meson K*. By means of the equation 
of motion, one obtains several relations between form factors 

a < \ m B + m K* . , s m B - m K * . , , 

M s ) = Ms) = A 2 {s), 

zttik* ZrriK* 
A o (0) = A 3 (0), 

{K*\d^\B) = 2m K *{e*p B )Ms), 

71(0) = T 2 (0). (2.10) 

All signs are defined in such a way as to render the form factors real and positive. The 
physical range in s extends from s min = imj to s max = (1 — rriK,K*) 2 - 

The calculation of the form factors given above is a real task, and one has to rely 
on certain approximate methods. We use the results calculated by using technique of 
LCSRs and given in |23| . And form factors can be parametrized as 



F(s) = F(0) exp( Cl s + c 2 s 2 + c 3 s 3 ). (2.11) 

The parameterization formula works within 1% accuracy for s < 15 GeV 2 and can 
avoid the spurious singularities at s — m 2 B . Related parameters is given in the Table. 
4ofO 



9.3 Formula of Observables 

In this subsection we provide formula of experimental observables, which include dilep- 
ton invariant mass spectrum, forward-backward asymmetry, and lepton polarizations. 

^From eqs. (2.2 - 2.8), it is straightforward to obtain the matrix element of B — > 
K(K*)l+l~ as follows. 

M = ^ V*V tb m B [r; (ir t) + T 2 (£Y 75 1) + S (ft)] , (3.1) 
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where for B -> 

3? = ^(S)^, (3.2) 

*rf = C\s)p, + D\s)q^ (3.3) 

5 = S 1 (s) (3.4) 

and for B -> K*£ + £~, 

r; = A(s)e^e*^p^-^(s)eV + iC(s)(e*-p B )^, (3.5) 

T; 2 = £(s) e^Oplfx* - iF(s) e% + (e* • p B )p M + itf(s) (e* • p B )g M (3.6) 

5 = i2m K *(e* ■ p B )S 2 (s) (3.7) 

with p = p B +Pk,k*- Note that, using the equation of motion for lepton fields, the 
terms in q^ in vanish. 

The auxiliary functions above are defined as 



A'(s 
C'(s 
D'(s 

A(s 

B(s 

C{s 

E(s 
F(s 
G(s 

H(s 



C 9 cS (s) f + (s) + T^C 7 ° s f T (s) , (3.8) 
i + rriK 

Cio /+(«)> (3-9) 

3.10) 



Cio/-(S) " . , C Q2 f (s) , 

2m e (m b - m s ) 



1 - 



Cqi/o(s) , 



(rh b - m s ) 

— — C^{s)V{s) + ^C^S) , 
1 + m K * s 

(1 +m K *) 
1 



C^s)^) + ^(1 - m K »)C 7 cS T 2 (s) 
(1 - m x *)^9 Cff (S)A 2 (s) + 2m b C 7 cS ( T 3 (s) + — ^ T 2 (s 



1 — m 2 K * 



2 -C 10 V(s) , 



1 + m K * 

— -\ — C w A 2 {s) , 
1 + m K * 

Cio 



s 



[(1 + rhK^A^s) - (1 - m^)^2(s) - 2m^. A (s)\ 



where 



3.11) 
3.12) 
3.13) 

14) 

3.15) 
3.16) 
3.17) 

3.18) 
3.19) 

3.20) 
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and to get the auxiliary functions given above, we have used equations of motion 



<f (^17/^75^2 ) 



(mi - m 2 )ipiip2, 
-(nix + m 2 )^i75^2- 



(3.21) 
(3.22) 



The contributions of NHBs have been incorporated in the terms of S\(s), D'(s), H(s) 
and £2(5) • K is remarkable that the contributions of NHBs in D'(s) and H(s) are 
proportional to the inverse mass of the lepton, and for the case I = fi, the effects of 
NHBs can be manifested through these terms. 

A phenomeno logical effective Hamiltonian is recently given in [29]. If ignoring 
tensor type interactions in the phenomenological Hamiltonian (it is shown that physical 
observables are not sensitive to the presence of tensor type interactions ||), it is easy 
to verify that the matrix element of B — > K^*H + £~ can always be expressed as the 
form of the equation (|3.1|) with the auxiliary functions defined as 



where 



A' (a) 
C(a) 
D'(s) 

S 1 (s) 

A(s) 
B(s) 
C(s) 



E(s) 
F(a) 
G(s) 

H(s) 



1 + m K 
1 



m 



K 



2m e (m b - m s ) 



w c Jo{s) 



w C2 f+{§) , 
w C2 f-(s) - 

1 — m 2 K 

(rh b - rh s ) 

2 2) 
— - - — w Cl V(s) - —{w C9 + w Cl0 )Ti(s) 
1 + mx* s 



-(1 + rh K 
1 



w^A^s) + -{I - m K *)(w C9 + w cl0 )T 2 (s) 



A — [(1 -m K *)w C3 (s)A 2 (s) 
1 — m K * 



+(w C9 - w cw ) 1^(1 + m K *)T 3 (s) + 

o 

-w C2 V(s) , 



1 — m 2 K » 



■T 2 (s) 



1 + m K 
-(1 + rh K *)w CA Ai(s) , 

-7—^ — w Ci A 2 (s) , 
1 + m K * 

^ 4 (A 3 (s) - A (s)) + 



w C7 A (s) 



rhf\m b + fh s 



w Cs A (s) 



{fn b + m s 



w, 



ci 



3.23) 
3.24) 
3.25) 

3.26) 

3.27) 
3.28) 



3.29) 

3.30) 
3.31) 
3.32) 

3.33) 

3.34) 

(3.35) 
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w c 2 


- \ Crr 


T t-ifi — LyRL + L<RR) , 


{6.6\J ) 


w C3 


l ( r 


- <-^LR + <-^RL + <~<RR) , 


(O r > r 7\ 

[6.6 1 ) 




= l( C LL- 


Cl_r — Ckl + Crr) , 


(3.38) 


w C5 


= -^{Clrlr 


+ CrllR + ClRRL + Crlrl) , 


(3.39) 




I 

= -^{Clrlr 


+ Crllr — Clrrl — Crlrl) , 


(3.40) 


w C7 


= -^{Clrlr 


— Crllr + Clrrl — Crlrl) , 


(3.41) 




= j{Clrlr 


— Crllr, — Clrrl + Crlrl) , 


(3.42) 


W C10 


= TYl^C BR , 

= m s C SL , 




(3.43) 
(3.44) 
(3.45) 



In the above equations Cll,Clr etc. are defined in ref. ||. Therefore our formula 
given below can also be used to make model independent phenomenological analysis, 
if using Eqs. ((fTEJ-flPD) in stead of Eqs. ((|3T8|)-(|3TT9D). 

Keeping the lepton mass, we find the double differential decay widths T K and T K 
for the decays B — > K£ + £~ and B — > K*£ + £~ , respectively, as 



2-rK 



d 2 r 



2 ™2™5 
B 



Gia 2 m 



dsd-u 



2rK* 



2 n 7r 5 



v; s v tb \ 



Vrti(U 



x {(|A'| 2 + |C"| 2 )(A - u 2 ) + |5i| 2 (s - 4m 2 ) + Re^A'*)^ 
+ |C"| 2 4m £ 2 (2 + 2m 2 K - s) + J Re(C" J D'*)8m/(l - m^) + |Zy| 2 4m/s} , (3.46) 



d 2 r 



dsd-u 



G\a 2 m\ 2 

I Via Vtfcj 



2 n 7T 5 

fUl 2 



x • 



+ 



4 
1 



(s(A + m 2 ) + 4m 2 A) + i—L (s(A + w 2 ) - 4m 2 A) + |<S 2 | 



4m 2 ) A 



|B| 2 (A 



u 2 + 8m 2 K *(s + 2m?' 



+ \F\ 2 ( A — w + 877^. (s - 4m 



c 2rti(U 



2su [Re(BE*) + Re(AF*)] + —^[Re(S 2 B*)(s + rh 2 K * - 1) + Re(S 2 C*)X] 



m K * 



+ 



A 



4m 2 K , 



Cf(A - m 2 ) + |Gf ( A - u 2 + 4m 2 (2 + 2m 



2m^* 



Re(5C*)(l - m 2 K * - s)(A - w 2 ) 
+Re(FG*) ((1 - m 2 ^ - I) (A - u 2 ) + 4m 2 A) 
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2^A 



rn 



K" 



Re{FH*) -Re(Gif*)(l - m 2 K .) + \H\ 



mi 



-sX 



rn 



K* 



Here the kinematic variables (s, u) are defined as 



which are bounded as 



with rhi = mi/m B and 



u(s 



A 
V 



§ = q 2 = (p + +p-) 2 , 

u = (Pb- p~y - (Pb - P+f 



(2m £ ) 2 < s < (1 - %,jf. 
— u(s) < u < ii(s) , 



1 + m^K* + s 2 ~ 2s - 2m\ K ,{1 + s) 



Ami 



(3.47) 



(3.48) 
(3.49) 



(3.50) 
(3.51) 



(3.52) 
(3.53) 

(3.54) 



Note that the variable u corresponds to 8, the angle between the momentum of the 
-B-meson and the positively charged lepton £ + in the dilepton CMS frame, through the 
relation u = —u(s) cos 6* |]4l[ . 

Integrating over u in the kinematic region given in Eq. ( |3.51| ) we get the formula 
of dilepton invariant mass spectra (IMS) 



dr 



K 



ds 

dT K * 
ds 

D K * 



G\a 2 m\ 

2 10 7T 5 



\V*V tb fu(s)D 



K 



u(s) 



+ \C'\ 2 4m e 2 (2 + 2m 2 K - s) + Re(C'D'*)8m e 2 (l - m 2 K ) + \D'\ 2 4m e 2 s 



G\ a 2 m% 



\V*V tb \ 2 u(s)D 



2 10 7T 5 

^5A(14 2—) 4 \hf, , 
3 s 3 



(3.55) 

(3.56) 
(3.57) 



m 2 , ,„, 2 „u{s^ 2 



+ |S 2 | 2 (s-4m 2 )A 



Am 2 K * 
A 



1 

2m^» 



|Cf(A- 



3 +8m^(s + 2m, 2 )) + |F| 2 (A- 



8m^,(s - 4m 2 )) 



?2(s)' 



) + \G\ 2 A 



+ 4m^(2 + 2m 2 K , - s) 



u\s 



Re(BC*)(X - ^V~)(l - - s) 
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+ Re(FG*)((A 



U(S 



-)(1 -m|* - s) +4m^A) 



-2- 



■A 



m 



K* 



Re(FH*) - Re(GiT)(l - m 2 K . 



+ -r^sA Lff 



(3.58) 



A"* 



Both distributions agree with the ones obtained in |23|, [35|], if Cq 1j 2 are set to zero. 
The differential forward-backward-asymmetry (FBA) is defined as 



A FB (s) 



U(S) dz dT 
dsdu 



, dT 
du- 



u(s) dsdu 



U(S) dz dF 
dsdu 



, dT 
du- 



u(s) dsdu 



For B — > K£ + £ decays it reads as follows 

-D' x = -2m^(s)Re(5iA / 



H A K 



ds 

For B — > K*£ + £~ decays it reads as follows 

u{s)\s [Re(BE*) + Re(AF*)] 



(3.59) 



A AK* 



ds 



+ 



rriK* 



Re(S 2 B*)(l - s - Th 2 K .) - Re{S 2 C*)X 



(3.60) 



We can read from ( 3.59 ), the FB asymmetry of the process B — > K£ + £~ does not 
vanish when the contributions of NHB are taken into account. With it, our analysis 
below also show the contributions of NHBs can even be accessible in B factories. 

The lepton polarization can be defined as follows 



dT(n) 
ds 



1 dT, 

2 W 



1 + (P L e L + P N e N + P T &r) ■ n 



(3.61) 



where the subscript "0" corresponds to the unpolarized amplitude, and Pl, Pt, and Pn, 
correspond to the longitudinal, transverse and normal components of the polarization 
vector, respectively. 

For the process B — > K£~£ + , the P^ , P? , an d P$ , are derived respectively as 



pK D K 



pKjjK 



P?D K 



H 

7T\/lu(s) 



ARe(v4 / C"*) - 3m/(l - rh 2 K )Re{C* Sx) - 3m e sRe(D'* S ± ) H3.62) 



2 



- Im(A'Sl) + 2m e lm(C'D' 



(3.63) 



mi 



[l-m 2 K )Re(A'C'*) 



+sRe{A'D' 



(s - Amj] 



■Re(C'Sl 



(3.64) 
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D K is defined in Eq. ([336]). For the process B -> K*£~£ + , the Pf *, Pf *, and 



P/v*, are derived respectively as 



PfD K * 



v\ 2 -^Re(AE*) + {X + 12 f K ' ) Re(BF*) 
I 3 3m A K , 

A< ' ~ ™K' ~ —Re(BG* + CF*) + -^-R e (CG*] 



3m 2 K * 



3m 2 K , 



2rho\ r 



r N u 



+ 

m K 

—irV§u(s) \ the 



4m 



K 



Re(FS;) - sRe(HS;) - (1 - m 2 Kt )Re(GS* 2 \ 
lm(FG*)(l + 3m 2 K , - s) 



(3.65) 



m K * 



+ lm(FH*)(l - rh 2 K , - s) - Im(GP*)A 

+ 2m K *m e [Im(PP*) + Im(AF*)] 

- (1 - m 2 K * - s)Im(PS*) + AIm(C<S 2 *) 1 



(3.66) 



PfD K * 



71 



rn 



K 



+ 



A 



m 



K* 



-Re(BF*) + (1 - rh 2 K *)Re(BG*) + sRe(BH*) 
Re(CF*) - (1 - m 2 K ,)Re{CG*) - §Re(CH* 



+ 



(s — Am, 



m K *m t 

D K * is defined by Eq.(ggg). 



[1 — m 



K* 



s)Re(FS;) - ARe(G<S*) 



(3.67) 



9.4 Numerical analysis 

Parameters used in our analysis are list in Table |2|. Considering that the branching 
ratios of B —>■ K £ + 1~ and B — > K* £ + £~ are not very sensitive to the mass of mj,, we 
neglect the difference between the pole mass and running mass of b quark. 

The Wilson coefficients in the SM used in the numerical analysis is given in the 
Table ||. C7 is defined as 



C e 7 ff 



C-j — C5/3 — Cq 



(4.68) 



Cq 1 2 come from exchanging NHBs and are proportional to tan 3 /? in some regions 
of the parameter space in SUSY models. According to the analysis in 0, 36], the nec- 
essary conditions for the large contributions of NHBs include: (i) the ratio of vacuum 
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m b 


4.8 GeV 


m c 


1.4 GeV 


m s 


0.2 GeV 




0.11 GeV 


mtau 


1.78 GeV 


M B 


5.28 GeV 


M k 


0.49 GeV 


M K * 


0.89 GeV 


M J/psi {M r ) 


3.10(3.69) GeV 


r B 


4.22 x 10~ 13 GeV 




8.70(27.70) x 10- 5 GeV 


T{J/psi -> £+£-) 


5.26 x 10~ 6 GeV 


T(4>' -> 


2.14 x 10~ 6 GeV 




1.17 x 10~ 5 GeV" 2 


a- 1 


129 


\v t *v tb \ 


0.0385 



Table 2: Values of the input parameters used in our numerical analysis. 



Ci 


c 2 


c 3 


c 4 




c 6 


cf 


c 9 


Cio 


C 


-0.248 


+1.107 


+0.011 


-0.026 


+0.007 


-0.031 


-0.313 


+4.344 


-4.669 


+0.362 



Table 3: Wilson coefficients of the SM used in the numerical analysis. 
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SUSY models 


R7 


i? 9 


-Rio 


Cqi 


Cq 2 


SUSY I 


-1.2 


1.1 


0.8 


0.0 


0.0 


SUSY II 


-1.2 


1.1 


0.8 


6.5(16.5) 


-6.5(-16.5) 


SUSY III 


1.2 


1.1 


0.8 


1.2(4.5) 


-1.2(-4.5) 



Table 4: Wilson coefficients of the SUSY used in our numerical analysis. Ri means 
Ci/Cf M . SUSY I corresponds to the regions where SUSY can destructively contribute 
and can change the sign of C7, but the contributions of NHBs are neglected. SUSY 
II corresponds to the regions where tan/3 is large and the masses of superpartners are 
relatively small. SUSY III corresponds to the regions where tan/3 is large but the masses 
of superpartners are relatively large. In the last two cases the effects of NHBs are taken 
into account. The contributions of NHBs are settled to be different for both the case 
£ — [A and £ = r, since Cq 12 are proportional to the mass of lepton. The values in 
bracket are for the case £ = r 

expectation value, tan/3, should be large, (m) the mass values of the lighter chargino 
and the lighter stop should not be too large (say less than 120 GeV), (m)mass splitting 
of charginos and stops should be large, which also indicate large mixing between stop 
sector and chargino sector. As the conditions are satisfied, the process B — > X s 7 will 
impose a constraint on C7. It is well known that this process puts a very stringent con- 
straint on the possible new physics and that SUSY can contribute destructively when 
the signature of the Higgs mass term \i is minus. There exist considerable regions of 
SUSY parameter space in which NHBs can largely contribute to the process b — > s£ + £~ 
while the constraint of b — > S7 is respected (i.e., the signature of the Wilson coefficient 
C7 is changed from positive to negative). When the masses of SUSY particle are rel- 
atively heavy (say, 450 Gev), there are still significant regions in the parameter space 
of SUSY models in which NHBs could contribute largely. However, at these cases C7 
does not change its sign, because contributions of charged Higgs and charginos cancel 
with each other. We will see it is hopeful to distinguish these two kinds of regions of 
SUSY parameter space through observing B — > K^*H + £~. 

As pointed out in HQ, the contribution of NHBs is proportional to the lepton 
mass, therefore for £ = e, contributions of NHBs can be safely neglected. While for 
cases £ = \x and I — r, the contributions of NHBs can be considerably large. To 
investigate the effects of NHBs in SUSY models, we take typical values of 67,9,10 an d 
Cq x 2 as given in Table £|. The SUSY model without considering the effects of NHBs 
(SUSY I in Table 4) is given as a reference frame so that could the effects of NHBs be 
shown in high relief. 



Numerical results are shown in Figs. 1-4. In Fig. |l|(a), the IMS of B — > K/j, + n~ 
is depicted. We see that at the high s regions, NHBs greatly modify the spectrum. 
While at the low s region, the effects of NHBs become weak. In Fig. [I](b), the FB 
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(c) (d) 

Figure 1: The IMS(a), FBA(b), LP(c), and TP(d) of the process B -> Kfi + fi~. The 
solid line, dashed line, dot line and dashed-dot line represent the SM, SUSY I, SUSY 
II, SUSY III respectively. Both the total (SD+LD) and the pure SD contributions are 
shown in order to compare. 



9 Exclusive Semileptonic Rare Decays B 



(K,K*)£ + £ in Supersymmetric Theoii&S 



asymmetry of the B — > Kfi + fi~ is presented. Fig. 0(b) shows that the average FB 
asymmetry in B — > Kfi + fi~ 0.02. To measure an asymmetry A of a decay with the 
branching ratio Br at the na level, the required number of events is N = n 2 /{Br A 2 ). 
For B — ► Kfi + fi~, the average FB asymmetry is 0.02 or so, the required number of 
events is 10~ 12 or so. Therefore it is hard to observe the derivation of FB asymmetry 
from the SM. In Fig. 0(c) and Fig. 0(d), the longitudinal and transverse polarizations 
are given. The effect of NHBs to the longitudinal polarization is weak but the effect 
to the transverse is remarkable. 



In Fig. 0(a) and Fig. 0(b) the IMS and FB asymmetry of B — > Kt + t~ are 
presented respectively. For SUSY II, the effects of NHBs to IMS is quite manifest, and 
the average FB asymmetry can reach 0.1. For SUSY III, the average FB asymmetry 
can reach 0.3. Therefore, in order to observe FBA, the required number of events 
should be 10~ 9 or so and 10~ 8 , respectively, so that in B factories, say LHCB, these 
two cases are accessible. In Fig. 0(c) and Fig. 0(d), the longitudinal and transverse 
polarizations are drawn respectively. The effects of NHBs are also very obvious. 

Figs. 3 and 4 are devoted to the decay B -> K*l + l~ . In Fig. §, the IMS, FB 
asymmetry, and polarizations of B — > K*fi + {i~ are given. We see that this process 
is not as much as sensitive to the effect of NHB as B — > Kfi + fi~. However, the 
contribution of NHBs will increase the part with positive FB asymmetry and will be 
helpful to determine the zero point of FB asymmetry. Fig. 0(d) depicts the transverse 
polarization of the B —>■ K*fi + fi~ , and the effect of NHBs is quite obvious. The zero 
point of the FB asymmetry can be slightly modified as shown in Figure 0(b) due to 
the contributions of NHBs. 



In Fig. 0, the IMS, FB asymmetry, longitudinal and transverse polarizations of 
the B — > K*t + t~ are depicted. The effect of NHBs does show in great relief. It is 
worth to note that IMS, FBA, and lepton polarizations for B -> K* £ + £~ in MSSM 
without including the contributions of NHBs are also significantly diffident from those 
in SM, while for B — > K £ + £~ they have little differences from those in SM. Therefore, 
compared to the process B — ► K £ + £~ , more precise measurements for B — ► K* £ + £~ 
are needed in order to single out the contributions of NHBs. 

Normal polarizations for both B — > K £ + £~ and B — > K* £ + £~ are small and can be 
neglected because the imaginary parts of Wilson coefficients are small in SUSY models 
without CP violating phases which are implicitly assumed in the paper. 

The behavior of IMS(a), FBA(b), LP(c), and TP(d) shown Figs 1-4 can be under- 
stood with the formula given in the Section 3. With Eqs. ( |3.56| ), (3.10) and (3.11), 
we see that the contributions of NHBs are contained in the terms of Si and D' . At 
the high s regions, it is these two terms which are important. This explained the 
behavior of IMS given in (a) of Fig. 1 and Fig. 2. The Eq. ( |3.59| ) shows that the 
FBA is proportional to the mass of the lepton. For the case B — ► Kfi + fi~, due to 
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Figure 2: The IMS(a), FBA(b), LP(c), and TP(d) of the process B -> Kt + t- 
line conventions are the same as given in the legend of Fig 1. 



The 
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Figure 3: The IMS(a), FBA(b), LP(c), and TP(d) of the process B -> K*fx + fi~. The 
line conventions are the same as given in the legend of Fig 1. 
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Figure 4: The IMS(a), FBA(b), LP(c), and TP(d) of the process B -> K*t + t- . The 
line conventions are the same as given in the legend of Fig 1. 
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model 


A 


B 


c 


D 


E 


tot(SD) 


tot(SD+LD) 


SM 


LCSR 


0.353 


54.707 


0.032 


4.566 


0.076 


0.573 


59.736 


SVZ 


0.215 


22.918 


0.015 


1.593 


0.026 


0.299 


24.767 


SUSY I 


LCSR 


0.425 


54.723 


0.037 


4.576 


0.086 


0.675 


59.847 


SVZ 


0.179 


22.910 


0.011 


1.586 


0.019 


0.236 


24.704 


SUSY II 


LCSR 


0.556 


54.865 


0.131 


4.833 


0.849 


2.067 


61.233 


SVZ 


0.348 


23.009 


0.068 


1.726 


0.321 


1.002 


25.473 


SUSY III 


LCSR 


0.429 


54.727 


0.040 


4.584 


0.109 


0.717 


59.889 


SVZ 


0.181 


22.912 


0.012 


1.590 


0.028 


0.255 


24.723 



Table 5: Partial decay widths for B — > K/j, + jj, . LCSR means the approach light- 
cone QCD sum rules, SVZ means the SVZ QCD sum rule [jT/j] . Character A 

~ 2 ~ 2 ~ 2 "2 "2 

means the region (s , (m^ — 5) ), B ((m^, — 5) , (m^ + 5) ), C ((m^ + 5) , (m^> — 5) ), 

" 2 " 2 * 2 _fi 

D((m^,/ — 5) , {rh^i + 5) ) and E ((m^/ + 5) , s max ). The unit is Tb X 10 , which is 

4.22 x 10~ 19 GeV. 5 is selected to be 0.2 GeV. 6 is normalized with M B 



smallness of the mass /i, the FBA does not vanish but is hard to be measured. While 
for the case B — > Kt + t~, the mass r is quite large and observing FBA is relatively 
easy. For SUSY II, though the numerator of FBA is comparatively large, the large 
IMS suppresses the value of FBA; for SUSY III, the numerator is relatively small, but 
the FBA do demonstrate the effects of NHBs more manifestly, as shown in Fig. 2(b) 
due to smallness of IMS. The Eqs. (3.63) and (3.64) show that for the case £ — fi, the 
contributions of NHBs to P^,Pt are suppressed by the mass of /i. But for the case 
I = t, the contributions of NHBs become quite manifest both for SUSY II and SUSY 
III. The term with D' in Eq. (|3.64|) will change its sign when there exists relatively 
not too small contributions of NHBs, the fact deduced from Eq. (3.10), that explains 
why the sign of TP is changed. The difference between the case SUSY II and SUSY 
III is small, the reason is just the same as stated in the analysis of FBA. 

Since the terms incorporating the contributions of NHBs is proportional to A as 
shown in Eq. (|3.58|) , which approaches zero at high s regions; while at small s regions, 
the effects of NHBs are dwarfed by the other contributions. Therefore, only when Cq 1 
are quite large could effects of NHBs be manifest, as shown in Fig. 3(a) and Fig. 4(a). 
According to the Eq. (|3.60f) , at high s regions, the effects of NHBs would be suppressed 
by A and 1 — s — rh\*. The same suppression mechanism exists for LP. This suppression 
mechanism explains the fact that the processes B — ► K* £ + £~ are not sensitive to the 
effects of NHBs. However, when there exist large contributions of NHBs, the sign of 
TP will be changed, as indicated in both Fig. 3(d) and Fig. 4(d). 

The partial decay widths (PDW) are listed in Tables. 5, 6, 7, and 8. We see that at 
the high s region, for the process B —>■ Kl + l~, l—fi,r, the contributions of NHBs do 
show up, as expected. For B K*l + l~ , the effects of NHBs in the high s region is 
signifiacnt when l=r while they are small for l=/i. It can be read out from these four 
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model 


A 


B 


c 


D 


E 


tot(SD) 


tot(SD+LD) 


SM 


LCSR 


0.930 


83.257 


0.141 


9.976 


0.258 


1.882 


94.562 


SVZ 


2.943 


111.278 


0.147 


7.504 


0.137 


3.639 


122.008 


SUSY I 


LCSR 


1.627 


83.402 


0.198 


10.085 


0.330 


2.915 


95.64 


SVZ 


4.517 


111.423 


0.183 


7.552 


0.149 


5.291 


123.825 


SUSY II 


LCSR 


1.178 


83.431 


0.234 


10.164 


0.352 


2.677 


95.360 


SVZ 


2.801 


111.292 


0.156 


7.525 


0.145 


3.522 


121.918 


SUSY III 


LCSR 


1.631 


83.407 


0.201 


10.092 


0.334 


2.938 


95.664 


SVZ 


4.518 


111.425 


0.184 


7.553 


0.150 


5.296 


123.830 



Table 6: Partial decay widths for B — > K*jj, + jj, . Other conventions can be found in 
Table 5. 



model 




A' 


B' 


tot(SD) 


tot(SD+LD) 


SM 


LCSR 


1.884 


0.094 


0.132 


1.978 


SVZ 


0.659 


0.036 


0.054 


0.695 


SUSY I 


LCSR 


1.884 


0.086 


0.131 


1.970 


SVZ 


0.655 


0.025 


0.038 


0.680 


SUSY II 


LCSR 


2.022 


1.496 


1.674 


3.519 


SVZ 


0.726 


0.552 


0.637 


1.278 


SUSY III 


LCSR 


1.874 


0.094 


0.129 


1.968 


SVZ 


0.651 


0.026 


0.035 


0.677 



-s 2 

Table 7: Partial decay widths of B — > Kt + t~ . Character A' means (s , (rh^ — 5) ), B' 

* 2 

means ((m^ + 5) , s max ). The unit is T B x 10~ 6 , which is 4.22 x 10~ 19 GeV. 



model 




A' 


B' 


tot(SD) 


tot(SD+LD) 


SM 


LCSR 


4.045 


0.096 


0.183 


4.141 


SVZ 


3.029 


0.048 


0.102 


3.076 


SUSY I 


LCSR 


4.088 


0.173 


0.327 


4.261 


SVZ 


3.052 


0.072 


0.159 


3.124 


SUSY II 


LCSR 


4.148 


0.266 


0.460 


4.413 


SVZ 


3.054 


0.084 


0.167 


3.138 


SUSY III 


LCSR 


4.078 


0.168 


0.312 


4.246 


SVZ 


3.050 


0.071 


0.156 


3.121 



Table 8: Partial decay widths of B — > K*t + t . Other conventions can be found at 
Table 7. 
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table that the results are consistent with the Fig. |I](a),|2](a), 0(a), and£|(a). In order to 
estimate the theoretical uncertainty brought by the methods calculating the weak form 
factors, we use the form factors calculated with LCSR and SVZ QCD sum rules (SVZ) 
method [17|. For B — > K £ + i~ , PDWs calculated with form factors obtained by SVZ 
method is 50% of those by LCSR approach; while for B — ► K* £ + i~ , PDWs increase 
100% or so. We see that at low s regions the theoretical uncertainty can reach from 
100% to 200%. Another point worthy of mention is that the contribution of resonences 



domainate the integerated decay width, as had been pointed out in [EE 



9.5 Conclusion 

We have calculated invariant mass spectrum, back-forward asymmetry, and lepton po- 
larizations for B — > K £ + i~ and B — > K* £ + i~ l=/i, r in SUSY theories. In particular, 
we have analyzed the effects of NHBs to these processes. It is shown that the effects 
of the NHBs to B — > Kt + t~ and B — > K*t + t~ in some regions of parameter space 
of SUSY models are considerable and remarkable. The reason lies in the mass of the 
r, which can magnify the effects of NHBs and can be see through from the related 
formula. The numerical results imply that there still exist possiblities to observe the 
effects of NHB in B — > Kfi + and B — > K*fi + {i~ through IMS, FB asymmetry and 
lepton polarizations of these processes. In particular, for B — > Kji + \i~ in the case 
of SUSY II, the partial width in the high s where short distance physics dominants 
can be enhanced by a factor of 12 compared to SM. Our analysis also show that the 
theoretical uncertainties brought in calculating of weak form factors are quite large. 
But the effects of NHBs will not be washed out and can stand out in some regions of 
the parameter space in MSSM. If only partial widths are measured, it is difficult to 
observe the effects of NHBs except for the decay B — > Kt + t~ . However, the conbined 
analysis of IBS, FBA, and lepton polarizations can provide usefull knowledge to look 
for SUSY. Finally, we would like to point out that FBA for B — > Kl + l~ vanishes (or, 
more precisely, is neglegiblly small) in SM and it does not vanish in 2HDM and SUSY 
models with large tan/3 due to the contributions of NHBs. However, only in SUSY 
models and for \—t it is large enough to be observed in B factories in the near future. 
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